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ABSTRACT 


i 

During July 1969 a reconnaissance magnetometer survey was conducted 
in the Valley of Ten Thousand Smokes with both total- and vert i cal -fie Id 
magnetometers . 

The large, sharp, narrow total magnetic anomalies observed over a 
zone of relict fumaroles in Broken Mountain Valley showed spectacular 
agreement with the surficial geology. Such a correlation is a strong 
indication that accumulations of magnetic minerals have been preserved 
along these fissure vents at shallow depths. Since large magnetic ano- 
malies were measured near fumarolic markings along all of the traverses, 
it is proposed that the retention of sublimates along fumarolic vents 
is common throughout the Valley. 

The generally concentric contours of the vertical magnetic anomaly 
at the head of the Valley suggest that the dome of Novarupta is merely 
the surficial expression of a very massive con i ca 1 -shaped intrusive 
centered just northeast of the dome. Corresponding offsets in the ano- 
malies along adjacent radial traverses, however, Imply the presence of 
concentric faulting around the dome, probably in the bedrock as well as 
in the pyroclastic flow. 

Profiles across the various branches of the Valley indicate that the 
flow is very heterogeneous. The variations in thickness and susceptibility 
implied by the total magnetic anomalies are consistent with the hypothesis 
of fissure-feeders for the flow. 

The magnetometer survey indicates that the pyroclastics in the Valley 
may be over 150 meters thick. Such an estimate is compatible with the 
volume of eruptive material needed to compensate for the subsidence 
surrounding Novarupta as well as a sizable amount of other regional sub- 
s I dence. 
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CHAPTER 1 


THE VALLEY OF TEN THOUSAND SMOKES 

1 . 1 Foreward 

On June I, 1912, the pyroclastic flow of the Valley of Ten Thousand 
Smokes and its many fumaroles, the collapse caldera of Katmai Volcano, 
and the plug dome of Novarupta volcano were formed by one of the world's 
largest recorded volcanic eruptions. 

After sixty years, many of the facts of this eruption are subject 
to controversy; chief among these are: the origin and mode of emplace- 

ment of the flow, as well as its thickness and internal composition and 
structure; the source of the fumarolic emanations and the chance that 
some of the fumarolic sublimates have been preserved; the cause of the 
banded ejecta which characte r i zes Novarupta and some portions of the flow 
the events leading to the collapse of Katmai caldera; and the cause of 
the marginal terrace throughout the Valley. 

1.2 Description of the Valley 

The Valley of Ten Thousand Smokes is situated about half-way down 
the Alaska Peninsula, within Katmai National Monument. it is bounded 
on the south by the volcanic peaks of Mounts Katmai, Trident and Mageik; 
on the west by the sedimentary Buttress Range; on the east by the vol- 
canics of Mount Griggs (formerly Knife Peak) and sedimentary mountains; 
and on the south by the Ukak River which runs along the base of sedi- 
mentary Mount Katolinat. 

Three format ional units are exposed in the vicinity of the Valley 
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f~y 1912 Pyroclastic flow 
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Figure- 1.1. Generalized geologic map of the Katmai area 
based on Keller and Reiser (1953, Plate 29) 
and Curtis ( 1968 , Figure l). Three major 
format i ona 1 units ore exposed in the Volley 
of Ten Thousand Smokes: the Naknek sedimentary 

formation, Quaternary Vo Iconics, and the 1912 
Pyroclastic Flow. 
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average thickness of at least 45 meters. The volume of these deposits 
has been estimated to be no less than 3-8 cubic kilometers (Kienle, 1970, 
p . 664 1 ) . 

The Valley is Y-shaped, being 17 kin long and from 4 km to 12 km 
wide. Figure 1.2 presents much of the nomenclature for the Valley. For 
ease in reference, the various branches of the Valley have been distin- 
guished as the lower valley, middle valley, southern branch and south- 
eastern branch. Between the southern and southeast branches are Novarupta 
Volcano and sedimentary Baked and Broken Mountains. Broken Mountain Val- 
ley lies between the latter two mountains, and heads at Greasy Pass 
(Griggs called this '‘Greased Hill", 1922, p. 24l) which connects these 
mountains just north of Novarupta. The Valley extending west of Nova- 
rupta is known as Novarupta Basin (Allen and Zies, 1923). Mageik Basin 
(Fenner, 1925) and Trident Basin refer to the minor depressions just north 
of each peak, respectively. The ridge separating Knife Creek Valley from 
Trident Basin was first referred to by Fenner (1923, pp . 3^“35) > who 
considered it to be the remains of a terminal moraine; it will be called 
Fenner Ridge. The southern peak of Broken Mountain which rises adjacent 
to Novarupta has been designated as Stumbling Mountain. 

The Katmai Trail traverses the Valley from the village of Savonoski, 
to Katmai Village via Katmai Pass. The pass lies between two old vol - 
can i c domes; Falling Mountain and Mount Cerberus. Both of these moun- 
tains are surrounded by the Valley of Ten Thousand Smokes pyroclastic 
flow which extends southward through the Pass into the valley of Mageik 
Creek. 

in the southeastern branch, Knife Creek heads at the base of Mount 
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Katmai where it issues from the Knife Creek Glaciers. The glaciers of 
Mount Mageik are the main source of the River Lethe in the southern 
branch of the Valley. The River Lethe has cut a deep gorge in the Upper 
Valley. Belov; Baked and Broken Mountains, both rivers have cut impres- 
sive narrow courses over 30 meters in depth. Except for the stream cuts 
and a marginal terrace called either the strand line or "high water mark", 
the general profile of the Valley floor is flat. 

The present volcanic activity in the immediate vicinity of the Val- 
ley is restricted to a few fumaroles. In 1969 , steam was issuing from 
near the summits of Mounts Griggs, Trident, Martin and Mageik, two pits 
at the margin of the Valley at the south-western edge of Baked Mountain, 
along the crater rim of Novarupta, and the fractures on the southern 
slopes of Broken Mountain. The ground was perceptibly warm at the vent 
of a small aromatic fumarole on the ridge line of Baked Mountain as well 
as at some of the fumaroles at the terminus of Broken Mountain Valley. 

1.3 History of Investigations within the Valley 

In October 1 83 8 , Spurr visited this region during a reconnaissance 
study for the U. S. Geological Survey. The account (Spurr, 1900) of his 
party's journey from Savonoski on the Savonoski River to Katmai Village 
on the Shelikof Strait via Katmai Pass yields the only geologic commen- 
tary on this area prior to the 1912 volcanic holocaust, it is most 
difficult to reconcile the present topography with the Katmai region 
prior to 1912 as shown on Spurr's reconnaissance map (Spurr, 1900, Map 
No. 11). The agreement of the map south of the Aleutian range is ex- 
cellent, but just to the north of the range correlation Is practically 
impossible until one reaches Naknek Lake. There has been some speculation 
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that Spurr did not traverse the Valley of Ten Thousand Smokes, but rather 
the valley of Windy Creek. This supposition is based on the angle at 
which the Katmai Trail continues north of the Pass on his map. Of parti- 
cular interest is the fact that the angle between the true strike of the 
Aleutian Mountains and that depicted on Spurr 1 s map is approximately the 
magnetic declination. It is known that Spurr's journey through the pass 
was hurried and such an error is unders tandab le . His party was delayed 
in Katmai Village, thus allowing time for better mapping south of the 
range. After studying Spurr's original field notes, Forbes concluded 
that Spurr did travel up the valley which was to become the Valley of 
Ten Thousand Smokes, though none of the evidence either way is conclusive 
(Forbes, personal commun i cat i on , 1971). The following conclusions are 
based on the supposition that Spurr did indeed map the Valley of Ten Thou- 
sand Smokes. 

The topography of the head of the Valley in 1898 was much different 
than today. Most obvious differences were the presence of a half-mile 
long lake on the northwestern side of the pass near the summit, which 
was dammed by the debris of three volcanoes rising above it, and a cone- 
shaped mountain adjacent to the lake on the west. Both of these features 
would have been located at the head of the present valley. One can assume 
that this mountain, which Spurr describes as "having a cone scarcely 
modified by erosion" (Spurr, 1900, p. Ik6) , was the ancestral Novarupta. 

It stood at least 1070 meters high (Spurr, I960, Map. No. 11), Origi- 
nally it would have encompassed, or at least overshadowed, Baked and 
Broken Mountains. 

Spurr remarks that the only instance in which the otherwise 
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horizontal green sedimentary strata forming the mountains bordering the 
valley showed any folding was adjacent to one of the volcanoes. In the 
Valley of Ten Thousand Smokes, the Naknek strata are folded where the 
southern end of the Buttress Range abuts Mount Mageik. 

Spurr describes the valley floor as relatively flat except for 
marginal terraces at 305 meters and 30 meters deep river gorges. The 
valleys were filled with glacial drift composed of stratified sands, 
gravels and even boulders near the pass, as well as some sedimentary rock 
fragments containing Jurassic fossils. Nearer the pass the surface was 
strewn with boulders, some forming sharp hillocks. The only unmodified 
drift was that in recently abandoned moraines. Several such moraines 
were damming mountain gorges. 

The National Geographic Society was responsible for the earliest 
investigations to the Katmai region following the 1912 eruption. In the 
summer of 1912, Martin (1913) visited Kodiak and Katmai Village. in 
1916, Griggs (1922) led a scientific party to the area, but not until the 
end of the field season did they venture north of Katmai Pass and dis- 
cover the Valley of Ten Thousand Smokes, Griggs returned in the follow- 
ing years to head the investigations in the Valley. In 1917, Shipley 
( 1920 ) took gas samples and measured the temperature of several fumaroles, 
studied the encrustations near some of the vents, and made a ground tempera 
ture profile across the terminus of Novarupta Basin. In 1918, Sayre and 
Hagelbarger (1919) continued the study of the temperatures of the fuma- 
rolic emanations. In 1919, Zies and Allen conducted even more extensive 
investigations of the fumaroles and their gasses and sublimates. The i r 
reports have become classics (Allen and Zies, 1923 ; Zies, 192Aa; Zies, 
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1 92 *ib ; Zies, 1929). In 1919, Fenner accompanied Griggs 1 expedition. 
Fenner returned to the Valley In 1923 to complete his studies. His geo- 
logic research resulted in a series of papers published over a span of 
30 years from 1920 to 1950. (Host of these are listed in the bibliogra- 
phy.) 

For many years no scientific parties visited the area; then a new 
round of studies began in the 1950's and is still continuing. In 1952, 
Wilcox of the U. 5. G. S. collected a suite of specimens from alteration 
zones adjacent to a fumarole at the terminus of the Valley. Lovering 
analyzed these samples for major and minor constituents (Lovering, 1957) • 
In 1953, the National Park Service promoted geological surveys within 
the KatmaS region. Curtis, Juhle and Williams took part in this re- 
examination of the Valley (Williams, 195*0. Curtis conducted a detailed 
study of the distribution of several distinct layers of tephra (Curtis, 
1968). Then, in the early 1960’s, a group of European investigators 
visited the Valley. Their views are presented in an article by Bordet 
et al . ( 1963) . 

With the construction of the Baked Mountain vol canolog i cal research 
station (BHM) a new series of investigations began. From this base, 
Kubota and Berg (1967) used seismic techniques to locate magma chambers 
in the area. Ward and Matumoto (I 967 ) studied the seismicity of the 
region and conducted a limited seismic refraction profile near the 
terminus of the Valley to test the applicability of hammer seismology 
in determination of the thickness of the py roc) as t i c depos i ts . Sbar 
and Matumoto (1971) conducted several such seismic refraction profiles 
across major branches of the Valley. Kienle (I 969 , 1970) surveyed four 
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gravity traverses across the main arms of the Valley. In 1969, Gedney 
et al, (1970) made two explosion seismic refraction profiles near Baked 
Mountain, This seismic work was continued in 1970 by Kienle and Bing- 
ham, who conducted two more profiles and in 1971 > Kienle added 5 wore 
short profiles. The ground magnetometer survey treated in this report 
was conducted in 3969. l n 1970 and 1971 > Anma and Stone continued the 
magnetometer survey, primarily from helicopter. 

1.4 Areas of Controversy 

a) The Source of the 1912 Pyroclastic Flow 

Immediately following the ash-falls of 1912 and the concurrent dis- 
appearance of the top of Mount Katmai , it was mistakenly assumed that 

/* 

Katmai was the sole source of the eruption. Later it became apparent 
that Novarupta had also contributed to the pyroclastic deposits in the 
Val ley . 

At first, the early travelers to the Valley of Ten Thousand Smokes 
thought the pyroclastic fill had been emplaced by a mudflow. Shipley 
{1920) proposed this sequence of events to account for the Valley fill: 
first, Novarupta exploded throwing out vast amounts of material; this 
ejecta fell on the snow-covered northern slopes causing much melt; 
heavy rain accompanied the eruption and aided this hot slush to slide 
into the valley forming a gigantic mudflow. The eruption of Mount Kat- 
mai followed the emplacement of the mudflow (Shipley, 1920, p. 141). 

The mudflow hypothesis was quickly abandoned once its obvious inade- 
quacies were exposed. The arguments against a mudflow are presented by 
Fenner (1920, pp . S77“578) . Basically, the existence of liquid water 
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at eruptive temperatures near incandescence would be incredible, 

Griggs (1922) postulated that a granitic bathol i th was approaching 
the surface over a wide area encompassing the volcanic range and the 
Valley, He believed that this mechanism could account for eruptions 
from Katmai and Novarupta, as well as from fissures throughout the Valley. 
A bathol i th would also res conti nued fumarolic activity. Griggs 

interpreted the eruptive events as beginning with the opening of many 
vents in the V alley floor and the release of lava through these. He 
supposed that the lava gave off so much gas following its release, that 
it became a fiery suspension of incandescent fragments buoyed up by the 
gasses they themselves were evolving. The masses accumulating about 
individual vents ran together until they covered the entire valley floor 
and then, under the influence of gravity, the entire mass poured down the 
Valley much like a flooding river. Next, Novarupta went into typical 
explosive eruption, followed closely by violent explosions from Mount 
Katmai . Meanwhile, the craters of Mounts Mageik and Martin supposedly 
opened too. 

The relatively undisturbed, horizontal Naknek sedimentary strata 
of the Valley walls led Fenner (1925b) to propose that the source of 
the pyroclastic flow was a sill intruded at shallow depths. The d i ml — 
nuation of furmaroHc activity by 1923 (estimated to be 1/3 of that 
in 1912 ) further indicated that the mass of hot material was not great. 
Fenner deduced that the source of the sill was a magma reservoir situated 
beneath the volcanic zone of Mounts Katmai, Trident, Mageik and Martin. 

As the magma rose from this reservoir, a portion found release northward 
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shattering the lower slopes of Mount Trident and fracturing the valley 
areas as they were raised relative to the surrounding more massive 
mountains. These fractures served as vents for the magma. Fenner (1920, 
p, 589; 1925b, p. 202) felt that Novarupta was similar to the other 
feeders for the tuff flow, although it was unique in that it broke out on 
a slope rather than along the Valley floor. Chance conditions were such 
that later this vent became enlarged and erupted great quantities of 
pumice and ash and extruded a lava dome. 

According to Fenner, much ejecta was also thrown out from Mount Kat- 
mai following the flow. At first Fenner (1925b, p, 201) believed the 
channel supplying magma to the Valley region from the chamber underlying 
the volcanic chain was not necessarily the same conduit supplying material 
to the Katmai crater; but in his last paper Fenner (1950b, pp. 707-708) 
acknowledged the interconnection of the vents and assumed that the Valley 
sill escaped from the Katmai conduit before the magma reached an explosive 
stage in the crater. Fenner (1920, p. 606) attributed the formation of 
the Katmai crater to collapse of the crater walls and incorporation of 
this material in the new magma. 

Fenner was the first to point out that the topographic regime of 
the Valley region would have prevented the present distribution of pyro- 
clastics had Katmai been the primary source. Indeed, Fenner proved that 
“from no single source .. .coul d the material well have reached all the 
areas where it is to be found" (Fenner, 1923, P* 17). Fenner (1950b, 
pp. 707-708) recognized Novarupta as a major source, but he believed 
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that much of the tuff flow erupted from fissures throughout the Valley. 

He did not appreciate the amount of subsidence around Novarupta al- 
though he did realize that the settling around the dome did seem refer- 
able to a collapse of the roof over the body of magma (Fenner, 1925b, 
p. 219). 

Williams (195^, pp. 58-59) was convinced that no sill was injected 
beneath the Valley and that the fumarolic gasses and associated subli- 
mates were derived from the fragmental ejecta itself. He attributes 
the pyroclastic eruptions to "glowing avalanches" issuing from swarms of 
fissures at the head of the Valley. These fissures were supposedly align 
ed along a zone essentially paralleling the volcanic axis, although off- 
set to the north. He also concluded that volumnous amounts of the two 
magmas were erupted from Mount Katmai , leading to the wholesale collapse 
of the summit to form the huge caldera. According to Smith (I960, pp . 
809-81 0) the pyroclastic fill of the Valley was erupted from fissures at 
the head of the Valley and emplaced by flowage of fragmental material 
which was itself continuously emitting hot gasses. Similarly, Bordet, 
Marinelli, Mi ttenpergher and Tazieff believe that the ignimbritic de- 
posits had been poured out through a swarm of fissures as an 1 'ove rf lowi ng 
glowing cloud" consisting of an emulsion of gas, glass splinters, hard 
particles of pumice, and intratelluric phenocrists. They were impressed 
by the apparent interconnect! ons underlying the Katmai district volcanoes 
They assume that these interconnections probably correspond to regional 
tectonic faults (Bordet et al . , 1963, pp. 7-8) . Ward and Matumoto ( 1 9 67) 
also attributed the main ash flow to Novarupta or fissures near the head 
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of the Valley. 

Finally, Curtis's ( 1 968 ) teph ra studies prove that Novarupta was the 
main source of eruptive material in 1912. Curtis distinguished 9 layers 
of tephra overlying the tuff flow. For A of these he was able to measure 
sufficient sections to prepare isopachous maps (a composite of these is 
presented in Figure 5.2). The contours of these maps close about Nova- 
rupta, and the general trend of mos t of the other layers also attests 
that the dome was their source. Apparently, only the last thin layers 
were erupted from Mount Katmal . 

Curtis thought that the layering of the tephra reflected individual 
eruptive events. He recognized the possibility that orientation of the 
conduit may have been partially responsible for the differing distribution 
patterns exhibited by successive tephra layers. However, he thought that 
the low correlative value of the stratification among various outcrops 
was due to variations in wind direction, speed and turbulence, and the 
eruptive stage. Since one of the first layers shows evidence of having 
been deposited by running water, Curtis concludes that temperatures 
during this part of the eruptions were sufficiently high to promote rapid 
melting of glaciers and snowfields adjacent to the Valley. 

Sbar and Matumoto's (1971) seismic refraction profile in Novarupta 
Basin shows a general thickening and greater complexity for the flow 
here than in the other branches of the Valley. This evidence also sup- 
ports the assumption that Novarupta was a major source of the flow. 

Curtis (1968) deduced that two magma chambers took part in the 
1912 eruptions. He believes that the rhyolitic chamber underlies 
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Novarupta, while the andesitic chamber underlies the arcuate chain of 
four volcanoes: Katmai through Martin. It appears then that the conduits 
of the two chambers were connected throughout most of the 1912 eruptions. 
Supposedly a column of magma supported the summit of Katmai during most 
of the eruption, allowing the summit to slowly collapse as the magma found 
release at Novarupta. When all activity had ceased at Novarupta and the 
two chambers were again distinct, a small amount of ash was erupted from 
Mount Katmai. In the waning eruptive stages of Novarupta, subsidence 
of the surrounding area began. The total subsidence seems to have amount- 
ed to over 250 meters. At the time of eruption therefore, Novarupta's 
vent would have been at a sufficient elevation to supply tuff to all areas 
where it is found. Curtis (1968, p, 192) concluded that “Novarupta and 
its radial and, possibly also, concentric fissure systems were the source 
vents for the great tuff flow . 11 Curtis ( 1968 , p. 19*0 believes that the 
continued activity of fissure fumaroles within a mile radius of the dome 
strongly attests to the existence of conduits leading to a magnetic 
source at depth, although no feeder dikes of any kind have been discovered 
in this area of disturbance. 

If it had had an elevation some 250 meters (800 feet) higher than 
today, Novarupta is the one location in all of the Valley of Ten Thousand 
Smokes from which the tuff flow could have reached all the places where 
it is found (Curtis, 1968, p. 192), Eruptions from this higher vent 
could easily have moved into all branches of the Valley and down the 
valley of Mageik Creek, too. Spurr's map of the Katmai region in 1 89 8 
substantiates the notion of a higher ancestral Novarupta. 
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b) The Banded Ejecta of Novarupta 

The banded structure of sane Novarupta lava and pumice has been the 
subject of niuch debate. According to Griggs (1922, p. 297) the new 
eruptive magma (rhyolite) dissolved the old rock (andesite) from the sum- 
mit of Mount Katmai by the process of “overhead-stopping". Fragments of 
the old rock in the process of sinking contaminated the newly extruded 
lava, resulting in the andesitic bands in the rhyolitic extrusives. Fen- 
ner (1950b, pp. 708-710) also believed that the primary eruptive magma 
was the rhyolite, but that the andesitic streaks are fragments of volcanic 
glacial debris and Naknek sediments which were attacked by the rhyolites 
as it was erupted from fissures in the valley floor, and through the de- 
trital material covering the valley floor. While the magma was at first 
quiescent in the vents, it would have assimilated the wall rocks of Nak- 
nek sediments and the overlying volcanic glacial drift. This contaminated 
lava would next have mingled with the eruptive rhyolite. Fenner there- 
fore considered that the dark scoria and inclusions represented partially 
digested rocks from the conduit walls. He was able to correlate spatial 
differences in the pyroclastic flow with the probable composition of 
the glacial drift at each location (Fenner, 1950b, pp . 707'710), 

Others attributed the andesitic inclusions to processes of magmatic 
differentiation. Fenner (1926, p. 772) thought that the crystallization 
process was incompetent to explain the form of variation of Katmai rocks. 
But, Forbes* (personal communication, 19 68) studies suggest that the 1912 
pyroclastics are differentiates of an andesitic magma which was subjected 
to a period of stillstand and subsequent fractionation by gravity settling 
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of mafics. The basic bands may be due to flow differentiation of mafic 
cumulates dragged up along the vent contacts during the terminal move- 
ment of magma along the conduit system (see also sections 2.3b and 5-3). 

Another view was proposed by Williams (195^, p. 58), who considered 
that the intermingled rhyolitic and andesitic ejecta were formed by simul- 
taneous discharge of the two magmas from the same or closely adjacent 
fissures at the head of the Valley. On the other hand, Curtis (19 68 , pp. 

and 207) proposed that a rhyolitic magma chamber underlying Nova- 
rupta was contaminated by andesitic lava from beneath Mount Katmai and 
its probable interconnections with other recently active volcanoes in the 
area. Under this regime, the andesitic lava would have reached Nova- 
rupta via a conduit which became closed in the later stages of eruption. 
The distribution of dark, bands in Novarupta dome suggests that the con- 
duit supplying andesitic magma to the rhyolite was sheet-like in cross 
sect i on. 

c) Nature of the Fissure Fumaroles 

The presence of numerous high temperature fumaroles for several years 
after the eruption was interpreted by Griggs (1922) as support for his 
hypothesis that a bathol i th was approaching the surface beneath the 
valley region. Similarly, the distribution of fumarolic activity led 
Shipley (1920, p. 1 A9) to conclude that "the mudflow is either in inti- 
mate contact with a heated mass of the earth's crust or the outlets for 
the gaseous emanations from the magma are well distributed beneath the. 
flow". Shipley (1920, p. 1A2) proposed that fumaroles are principally 
located along cracks in the mudflow which had been formed as the mudflow 



dried and contracted over an uneven bedrock topography. 
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Allen and Zies also believed that certain features of the fumaroles 


indicated the presence of a lava body beneath the old valley floor. They 
attributed the rectilinear alignment of the fumarolic cracks to their 
deep-seated fissure origin. They could not believe that a reasonable 
(about 60 meters) thickness of pumiceous material would be able to 
account for the high temperatures persisting in many fumaroles through- 
out the Valley in 1919, and since the structure of the pumice showed that 
it had lost most of its gasses during eruption, it could not be the 
source, of the amounts of gasses still being exhaled in 1919 (Allen and 
Zies, 1923, p. 95). Furthermore, analysis of the metallic sublimates 
from some of the fumaroles as contrasted to the unaffected ash evidenced 


that the prevalent magnetite could not have been derived from the sur- 
rounding pumice (Zies, 192ka, p. 166). 

Fenner (1925b) interpreted the dlminuation of fumarolic activity by 
] 323 as support for his hypothesis that a sill, not a batholith, had 
been intruded beneath the Valley in 1912. According to Fenner, bedrock 
fissures served as feeders for the tuff flow. The distribution of the 
fumaroles was similarly controlled. Such a deep-seated fissure origin 
of fumaroles was challenged by Williams (195*, P- 58) who observed that 


the fumarol e cracks gradually disappear downward within the avalanche 
deposits. In concordance with Williams, Curtis (1968, p. 186) reports 
that his examination of all the bedrock exposed below the ash revealed 
no feeders for either the tuff flow or the fumaroles. The old fumarole 
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were observed to extend more than 12 meters into the tuff flow. 

In a recent seismic survey by Sbar and Matumoto (1969, p. 3^0) , 
however, it was noted that in a number of places where the refraction 
profiles crossed fumarole lines, a discontinuity in the deeper horizons 
was observed. To them, this correspondence suggested that the vents 
penetrated to the old valley floor. Gedney et al. (1970, P. 2623) sug- 
gest that the correspondence of fissure fumarole lines with a fault 
in the lower seismic horizon shows that some fumarolic activity did ori- 
ginate through bedrock faults, probably from the degassing of residual 
magma in a subsurface reservoir. 

Correlation or the seismic, gravimetric and magnetic data across 
the mouth of the southern branch of the Valley (see Figure ^.10) suggests 
an undulating, if not faulted, surface for the bedrock. Since the ex- 
posed Naknek strata are nearly horizontal in this region, and Spurr (1900) 
reported the ,1898 valley as generally level, one must conclude that the 
floor was subsequently faulted; perhaps in' connection with the 1912 erup- 
tion as suggested by Fenner (1925b, pp . 20^-206). The positioning of the 
fumaroles above steep bedrock slopes in no means proves that the vents 
have their "roots" in the bedrock, or below. Rather, it may only suggest 
that the bedrock relief controlled the establishment of fumaroles by in- 
fluencing the location of faulting within the flow, by guiding the ris- 
ing vapors, and by concentrating the ground water. 

The general consensus (lovering, 1957; Smith,. I 96 O; Ward and Matu- 
moto, I 967 ) is that the majority of the valley fumaroles were generated 
by degassing of the pyroclastics and the vaporizing of ground water. 
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This process could account for the waning of fumarolic activity as the 
deposits cooled. 

A deep-seated origin, however, is suggested for some of the fuma- 
roles at the head of the Valley. Griggs (1.922, pp . 238-239) describes 
several locations at the head of the valley where the fumaroles were 
observed to extend into undisturbed sandstone strata. Ward and Matu- 
moto (I 967 ) contend that the prominence and persistence of fumaroles 
along the bench (strand line) on the western side of Baked Mountain re- 
quire more explanation than merely sweating out of the flow, which would 
necessarily thin near the edge of the Valley. According to Forbes (per- 
sonal communication, 1 3 68 ) , the fumaroles which are now active along 
linear trends on Broken Mountain, the summit ridge of Baked Mountain 
and the solfatara fi el <1 adjacent to Novarupta are degassing through fis- 
sures which cut the underlying sediments. The persistence of the fuma- 
roles at the head of the Valley have been attributed to the presence of 
magma in the subsurface in this area (Lovering, 1957; Smith, I960; Forbes, 
personal communication, 1988). 

Allen and Z i es (1923, p. 152) report that the highest fumarolic 
temperatures in 1919 followed a discontinuous zone about Baked and Broken 
Mountains with an extension out into the middle valley (refer to Figure 
5.4). This zone is also marked by faulting, perhaps related to the 
subsidence of Novarupta in the waning eruptive stage. Lovering (1357, 
pp. 1586, 1588) 'also remarks that the areas richest CO^ and the sulfur 
acid gasses in 1 9 1 9 surround Baked Mountain and Broken Mountain near 
Novarupta. Sulfur enrichment is characteristic of active volcanism. 
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not later degassing of effuslves. Drainage of a magma chamber under- 
lying the head of the Valley through the major vent at Novarupta could 
have led to settling of the Baked-Broken Mountain complex, as well as to 
the obvious subsidence in the immediate vicinity of Novarupta. Such 
collapse has been postulated by Fenner (1925b, p. 210). Escape of mag- 
matic gasses along such an encircling zone of bedrock faulting could be 
responsible for the elevated temperatures and sulfur content peculiar 
to this area in the early years after the 1912 eruption. 

The seismic recordings of Kubota and Berg (19&7) suggest the pre- 
sence of several magma chambers in the neighborhood of the Valley. 

Magma chambers were located at intersections of ray paths (epicenter to 
recording station) which shew no S-phase, indicating transmission through 
a medium of low rigidity. Magma chambers were found beneath Mounts Kat- 
mai and Trident., Mount Griggs, and Mounts Mageik and Martin. None of the 
data collected indicate the presence of magma beneath the Novarupta 
area; however, one of the inadequacies of the mode of location is that it 
cannot detect small pockets near or beneath the recording stations and 
one of their stations was operated from the Baked Mountain Research Sta- 
tion. 

d) Possible Concentrat i ons of Magnetic Minerals 

It has been much debated whether or not any of the abundant accumu- 
lations of magnetite observed throughout the Valley in the first few 
years of intense fumarolic activity could have been preserved. For 
example, only one magnetite accumulation has been described in the 

in 1913* 2ies (1924a, p. 166) observed "loosely coherent 


1 i terature . 
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octahedra of almost pure magnetite" forming the lining of the roof of a 
series of fissure vents occupying an area approximately 100 ft by 100 ft 
(30.5 m x 30.5 m) at fumarole 148 in the middle Valley. The deposit was 
concentrated along ten parallel fissures. "The depth of the deposit 
varied f rom four to nine inches.. .the width varied ft om twe 1 ve to f i f teen 
inches" (zies, 1929 , p. 16), and "in places the deposit was four inches 
thick" (Zies, 1924a, p. 1 66 ) . Fenner revisited this site In 1923 and 
reported the absence of visible magnetite (Zies, 1929? P* 1 6 ) . 

Both Shipley (1920) and Zies (1924a) recognized that although ex- 
posure to high temperature emanations could increase the Iron content 
of the ash nearby, prolonged exposure tended to decrease the i ron-conten t . 
Zies (1929) offered this explanation: rising hot acid gasses may have ex- 
tracted metallic constituents from their source areas, and also along 
their routes. As these emanations cooled, metallic sublimates would form 
a lining in the conduits, many encrustations would be built up around 
the vents, and the nearby ash would be altered. Fumarolic magnetite and 
hematite, formed by the hydrolysis of iron transported as a halide in 
the vapor phase, would thus be concentrated in the upper portions of 
high temperature vents and under proper pressure conditions, as the 
temperature of the emanations dropped, magnetite might even have been 
deposited along the vents at considerable depths. When the temperature 
of these fumaroles fell to a point where active condensation of the acid 
steam could take place in or near the vent, however, the acid gases would 
go into solution and leach the surrounding area, decomposing the magne- 
tite and releasing the iron unless the conduit had developed some form of 
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protective cover. Certainly many sublimates were leached by the acid 
solutions or removed by runoff over the years. 

If the vents had become blocked off early in their history, and thus 
protected from the ready access of surface waters, a deposit resembling 
a mineral vein might result (Zies, 1929, p. 60; Fenner, 1923 , p. 51). 
Early investigators in the Valley report instances. in which vents were 
blanketed by deposits from heavily laden waters (Griggs, 1922, p. 245), 
or filled in by the constantly shifting surficial ash (Fenner, 1925b, 
pp. 206-207). In other instances, surficial accumulations of fumarolic 
clays could certainly have formed an effective seal over some of the 
vents . 

In a recent study of the alteration zones around fumarole No. 1, 
Lovering (1957, p. 1596) discovered that the iron content was high just 
outside the inner lining of the vent. His findings are presented in 
Table 1.1 and Figure 1.4. Apparently here the magnetite and hematite, 
which had been precipitated by hydrolysis early in the life of the 
fumarole, were protected from the leaching acids as the temperature of 
the emanations dropped. 

Accumulations of iron have been observed in other ash flow deposits. 
Mackin (1952, p. 1333) di scovered veinlets of crystalline hematite occur- 
ring in joints and minor breccia zones in the upper lithoidal unit in the 
ignimbritic deposits of the Iron Springs District of Utah. Gilbert (1938 
p» 1851) identified the grains of magnetite and hematite at the top of 
a welded tuff in eastern California as products of sublimation of gasses 
In extensive studies of the Bishop Tuff, a region apparently similar to 
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Table 1.1 

Variation (Weight Percent) of Ferrous- and Ferric-oxide 
in Samples from Alteration Zones about 
Fumarole No. 1, from Lovering (1957, Table 1, p. 1593). 
Refer to Figure 1.4 for explanation of zones. 


Zone 

1 

2 

3 

4 

5 

6 

7 

£ Fe 2 0 3 

1 .22 

4.59 

3.11 

3.20 

1.15 

l .99 

1.33 

% FeO 

.58 

-74 

1.41 

• 99 

.59 

.90 

1.38 


Table 1 .2 

Variation of Ferrous- and Ferric-oxide in the 
Bishop Tuff, after Sheridan (1970, Table 1, p. 864) . 


Average for 
una 1 tered tuff 

% Fe 2 0 3 0.43 

% FeO 0.19 


Average for 
unaltered fuma roles 

0.93 

0.24 


Average for rare 
intensely altered inner 
zone of some fuma roles 

4.6 

1 .53 



Orifice 



Figure 1.4. Variation of ferrous, ferric and total Iron in alteration zones around 

fumarole No. 1 (Zies, 1929) in the Valley of Ten Thousand Smokes in 1952 
(After Lovering, 1957, Figures 3 and 5). Refer to Table 1.1. 
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the Valley of Ten Thousand Smokes, Sheridan (1970, pp . 860 and 86^) dis- 
covered that the great majority of fumarolic mounds have no significant 
overall chemical difference from the surrounding tuff, although their 
iron content is increased. A few fumaroles showed an inner zone of ex- 
treme black and red discoloration about their central joints. The average 
iron content was found to vary among unaltered tuff, fumarolic tuff, 
and the rare inner zone altered tuff. These variations are given in Table 
1.2. Sheridan also noted that although some fumarolic fractures are 
coated with hematite and opal, minerals that might have formed during the 
early fumarolic stage are now missing from most fumaroles. He attributes 
this absence to probable attack by acidic fluids in the latter part of 
fumarolic activity. 

There may also be some secondary concentrations of magnetite. Kienle 
discovered several small accumulations of magnetite grains in depressions 
near streams where they we re apparently deposited during periods of large 
runoff (Kienle, personal communi cat i on , 1970). 

e) The Cause of the Marginal Terrace 

Almost throughout the Valley is a marginal terrace about 100 meters 
above the valley floor. The cause of this terrace has been much debated. 
According to Shipley (1920), the terrace marks the highest level of the 
flow; thus the synonym "high water mark". Others thought that the 
terrace was formed as the center of the flow compacted after coming to 
rest. Fenner (1929b, pp . 20^-206) did not believe there was a sufficient 
volume of pyroclastic fill to account for the strand line by mere com- 
paction, besides, he believed there was evidence that the flow was indurated 
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prior to any downward movement of the valley floors. Fenner proposed 
that as the magma was drained through the sill into the valley regions, 
the area readjusted. The terrace, or strand line, was formed along the 
valley walls as the bedrock settled over the body of the intruded magma. 

Curtis (1968, p. 186) attributed the strand line or "high water 
mark" around much of the margin of the Valley to gradual compaction of 
the tuff flow after the emplacement. He concluded that probably the 
compaction was proportional to thickness and thus would have been great- 
est over the old stream channels; consequently, the location of present 
streams would be over the old channels. Partial welding within the flow 
(Smith, i960), subsequent melting of an ice lens buried during the erup- 
tion (Hamilton, personal communication, 1970 » or the existence of 
buried glacial terrace (as described by Spurr, 1900) are other explana- 
tions for the configuration of the cross-valley profile. 

f) The Age of the Lower Tuff Unit 

There has been some speculation that some of the tuff in the Valley 
of Ten Thousand Smokes originated prior to the 1912 pyroclastic flow. 

Fenner (1923) concluded that some of the tuffs which had been thrown 
out by explosion craters at the base of Mt . Griggs represented a pre-1912 
tuff. These particular tuffs, although similar in appearance to the re- 
cent tuffs, are much more indurated and their inclusions are much more 
decomposed. Fenner regarded these as representative of "a hardened tuff 
of considerable age that formed a layer at or near the old surface along 
a stream or wet piece of ground at the time when the sandflow spread 
over it" and he supposed that "the ensuing explosions broke it up and 



threw it broadcast 11 (Fenner, 1923, p. 23). 

Curtis (1968, p. 204) describes two sections of tuff within Katmai 
caldera which are evidently of pre-1912 origin. Forbes (personal commu- 
nication, 1971) also suggests that the lower tuff unit exposed at several 
locations in the upper Valley, in particular at Fissure Lake, may well 
be of pre-1912 origin, 

g) Thickness of the Flow 

The objective of most of the invest igat ions in the Valley of Ten 
Thousand Smokes has been estimating the total volume of the pyroclastics 
erupted in 1912, The primary unknown factor in such a determination is 
the thickness of the flow. The several estimates for the thickness of 
the tuff are given in Table 1.3. As more geophysical evidence is gathered, 
the more conservative become the initial predictions of Fenner (1923, p. 
67 ); but besides the geophysical evidence there is another, to date 
apparently ignored, indication of the thickness of the deposits. Spurr 
(1900, p. 1 A6) describes the pre-1912 Valley as generally level except 
for river gorges and a distinct horizontal terrace about 305 meters above 
the floor. If one assumes that the so-called "strand line" is indeed an 
expression of this terrace since buried beneath the pyroclastic material, 
it is possible to deduce a minimum thickness for the flow. Today the 
terrace along the western flank of Baked Mountain is approximately 
100 meters above the level of the present floor. It is clear therefore 
that the pyroclastic deposits in this branch of the Valley are at least 
200 meters thick. Since the terrace also is covered by a mantle of ash, 
the flow in the valley must be thicker by this amount of ash. Also, any 



Table 1.3 


Previously Suggested Maximum Tuff Thickness (meters) 


! nves 1 1 gator 

Method 

Lowe r 
Valley 

Middle 
Va 1 1 ey 

Southern 

Branch 

Southeast 

Branch 

Broken Mtn . 
Valley 

Novarupta 
Bas i n 

Fenner 

i nduct i on 

— 

30 

60 

60 

— 

— 

Curtis 

geomorphology 

120 

200 

225 

150 

50 

— 

Ki en ie 

gray i metrics 
(Mode 1 P) 

65 

— 

70 

69 

25 

25* 

Ki en le 

gravi met r i cs 
(Model F) 

140 

— 

150 

170 

6o 

50* 

Matumcto 
and Ward 

seismic 
refract 1 on 

46 

— 

— 

— 

— 

— 

Gedney 
et a 1 . 

seismic 
refract i on 

— 

— 

CM 

— 

36 

— 

Sbar and 
Matumoto 

seismic 

refraction 

— 

— 

92 

22 

27 

80 

Ki en 1 e 

seismic 
refract t on 

>46 

>63 

... 

>100 

104 

>34 


-value queried since no bedrock reference data for these calculations. 
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subsidence of the valley floor during the eruption would have the effect 
of increasing the thickness of the flow above the 200 meter estimate. 

By imposing the gradients of adjacent valleys to the Valley of Ten 
Thousand Smokes according to the elevation of the base of the flow at the 
end of the valley, Curtis (1968, p. I87) recons tructed probable pre -1912 
Valley profiles. Based on this reconstructed profile, the pyroclastic 
fill is probably between 700 and 900 feet thick and would amount to -some 
2.63 cubic miles (11 km ) . 

In 1966 , Kienle ( 1 369) completed four gravity traverses across the 

main brarjches of the Valley. He determined that the minimum average 

thickness of the pyroclastic deposits is fairly uniform and ranges from 

35 to 40 meters. His estimates show minimum ash thickness of 6 to 8 

meters over the ridges of Baked and Broken Mountains, respectively (Kienle, 

1989, p. 138), and maximum thickness of 70 meters over the buried river 

channels (Kienle, 1970 , p. 66 - 47 ) - Kienle ( 1970 , p. 6659 ) estimates the 

3 

volume of the flow to be between 3.8 and b ,~j knr , depending upon the 
assumed density contrast. 

Various thicknesses have been suggested from seismic refraction sur- 
veys. Sbar and Matumoto ( 1971 ) determined that the main body of the tuff 
varies in thickness from 20 to over 70 meters. Gedney et a). (1970) 
found a maximum thickness of 50 meters for the pyroclastics in the 
southern arm of the Valley, and of 25 meters in Broken Mountain Valley. 
Kienle's (personal communication, 1971 ) profiles imply tuff thicknesses 
varying from 50 meters in the lower valley to over 100 meters in Broken 


Mountain Valley. 



The reduced data from all available seismic refraction profiles is 
presented in Appendix D. The implications of these data are further 
discussed in Chapters IV and V. 

1.5 field Procedure for the Ground Magnetometer Survey 

In the hope of answering some of the many questions about the pyro- 
clastic flow in the Valley of Ten Thousand Smokes, a ground magnetometer 
survey of the region was undertaken in 19&9. 

The purpose of magnetometer surveying is to ascertain local spatial 
variations in the magnetic field which can be related to local geologic 
structure. It Is required, then, to remove the gross effects of the 
geomagnetic field and its diurnal variations from all survey data. Since 
the temporal changes in the Earth's field can be considered constant with 
in an area the size of the Valley of Ten Thousand Smokes (Grant and West, 
1965, p, 207 ), it Is possible to use the record of the dally fluctuations 
in the field taken at a local base station to apply a diurnal correction 
to all magnetometer survey data. Such a magnetometer base station was 
established at the Katmai Vol canol og i ca 1 Research Station on Baked Moun- 
tain. It is marked by a wooden post ups lope and to the west of the 
Baked Mountain Hut. Diurnal fluctuations in the geomagnetic field 
were monitored at this station with a Varian Model V-493S Rubidium Vapor 
Magnetometer. The sensing head of the magnetometer was buried adjacent 
to the post designating the magnetometer base. The Lamor frequency out- 
put of this instrument was monitored from within the camp. Due to mal- 
function of the recorder, the Lamor frequency display had to be manually 
recorded, usually at half-hour intervals during field survey. The Lamor 
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frequency can be converted to magnetic intensity to a precision of 0.01 
gammas . 

The magnetometer base station also served as a check point for all 
the field survey instruments. In order to be able to compensate for 
possible instrumental drift, the various instruments were read at the 
base prior to and following each day of surveying. 

Relative spatial variations in the vertical magnetic field were 
measured with a Clsco-Sharpe Model MF1-100 fluxgate magnetometer, having 
a precision of + 2 gammas. The instrument displays the vertical field 
In gammas relative to an arbitrary zero. This arbitrary zero position- 
ing permits direct comparison with any chosen datum; in this case the 
fluxgate was set at 150 gammas at the Baked Mountain Hut Magnetometer 
Base Stat ion. 

Since the fluxgate magnetometer is a compact, light-weight instrument, 
easily carried and operated by a single person, it was ideal for less de- 
tailed work and for “spot readings" on a reconnaissance basis. Due to 
the high magnetic latitude of Katmai National Monument (magnetic inclina- 
tion of 71-2 degrees) anomalies in the vertical field are comparable to 
anomalies in the total field. 

An El sec proton precession magnetometer type 592/132 FS was used 
for survey measurements of the total geomagnetic field and its vertical 
gradient. The recording unit displays the number of proton precessions 
per three second interval. This value converts to gammas of field strength 
within 0.1 gammas. The adjustable sensor probe support facilitated 


readings at two heights for determination of the vertical gradient of 



the total field. 


The proton precession magnetometer was read at least twice at each 
station to ensure that correct values were recorded and to guard against 
errors caused by short term fluctuations in the geomagnetic field. If 
the two readings were not In close agreement, the station was occupied 
until a steady reading was obtained. This procedure was repeated for 
both upper and lower probe positions if the vertical gradient was being 
measured. 

A Brunton and a Meridian magnetic compass were used to take magnetic 
bearings along the traverses, and to shoot angles to prominent features 
for location of specific stations. Isolated "spot reading" stations, and 
most initial and final traverse stations, were located by t r i angul at 1 on 
on at least three we 1.1 -separated landmarks, weather permitting. In the 
cases where low cloud cover precluded such t rl angu 1 at i on , the location 
of the station was found only by reference to the U. S. G. S. map A--4, 
Katmai Quadrangle. 

During the surveys, straight line traverses were maintained by 
visual sighting back along flagged stakes marking previous stations, as 
well as by compass bearing and fore-sighting to landmarks. 

Total field stations were spaced along the surface at intervals 
of small integer (usually one or two) multiples of 28.05 + 1.00 meter, 
which corresponds to the length of the electric cable between the sens- 
ing head and recorder of the Elsec magnetometer. When adverse topography, 
such as stream beds, dictated adjustment of this standard station spacing 
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scheme, the offset was paced or estimated. 

The vertical field traverses, except for OC 1 below, were paced. 

In the cases where the fluxgate traverses follow previous total field 
traverses, (e.g., B-B", B-B tM , b-b") the distance between nearby paced 
vertical field stations and the flagged stakes marking total field sta- 
tions was recorded. This procedure permitted the proper coincidence 
of the two data sets. Stations along traverse C-C 1 were spaced with 
a fifty-foot (15.2 m) length of nylon cord. This traverse, which coin- 
cides with one of the seismic survey lines (Gedney et al., 1970; seismic 
profile A-A‘), was preceded by a transit theodolite survey. The esti- 
mated or paced distance between surveyed positions (marked with small 
stakes) and nearby vertical magnetometer stations was recorded, thus 
allowing verification of location and spacing of the magneotmeter sta 
tions. It was also beneficial to correlate the elevation differences 
derived from the level survey with the elevation recorded by the alti- 


meters during the magnetometer traverse. 

Usually, altimeter survey accompanied the magnetometer traverses, 
with altimeter elevations recorded at selected stations. 

The altimeter (one of two Paulin altimeters having precisions of 
+ 2.5 and + 5 feet) was set to 2550 feet at the magnetometer post before 
each survey. Marked breaks in slope, stream beds, etc., were noted in 
the field notes to supplement the altimeter data in the production of 
elevation profiles. For logistic reasons, several of the traverses near 
Novarupta were conducted without an altimeter. In the cases of traverses 
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elevations were taken from the U. S. G. S. map A-4, Katmai Quadrangle. 
Clearly, elevations obtained in this manner lack precision. 

The time of occupation for selected stations was recorded, together 
with brief descriptions of the surface topography and geology, including 
notes on the crossings of old fumarole lines, nearby inactive fumaroles, 
position of snow fields* etc. 

Several photos were taken in conjunction with the conduct of the 
magnetic survey (and collection of susceptibility samples). These were 
an invaluable aid in the months of analysis following the field operations 

The reduced magnetometer data and other various details of the ground 


survey are given in Appendix A. 



CHAPTER I I 


MAGNETIC SUSCEPTIBILITIES FOR ROCKS 
FROM THE KATMAI AREA 

2.1 Purpose of Magnetic Susceptibility Determination 

The parameter distinguishing rock type in magnetometer surveys 
is the magnetization. In recent deposits such as the Valley of Ten 
Thousand Smokes pyroclastic flow, the magnetization can be assumed to 
be in the direction of the present geomagnetic field, thus the defini- 
tive property of the rock types in the study area is their magnetic 
suscept I b i 1 i ty . 

An invaluable aid in deducing the structure and composition of 
geologic bodies responsible for observed magnetic anomalies Is computer 
modeling. The magnetic susceptibility assigned to the model body pro- 
ducing an anomaly like the observed anomaly implies a composition. Ob- 
viously, for the interpretive method used in this study, it would be most 
convenient were each rock type in the region to possess a unique suscepti- 
bility. Model studies could then distinguish the rock type responsible 
for the observed anomalies. 

Hand samples were collected in the Valley region of i960. Among 
these are samples of banded Novarupta rhyodacite, Lethe River vitroclastic 
tuff, Nakenk sediments and surficial Valley deposits. Unconsolidated 
ash was collected in polyethylene bags from a cut bank along Knife Creek, 
the eastern stream cut in Broken Mountain Valley, and a pit fumarole in 
the fumarole area in Broken Mountain Valley. Locations of these sampling 
sites are shown In Figure 2.1. Subsequent sampling in 1970 and 1971 has 
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Figure 2. 


Location of magnetic susceptibility sampling 
sites 1 through X. For data from individual 
sites, refer to Table 2.1. 
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made available more types of rocks from additional sites for suscepti- 
bility determination. The susceptibility data are given in Tables 2.1 
and 2.2. 

2.2 Determinat ion of Susceptibility 

Magnetic susceptibilities were measured on a susceptibility bridge 

developed from that described by Collinson et al . (1963). Bulk suscepti- 

~8 

bilities down to about 5 x 10 emu/cc can not be measured with this system. 
The field In which the rocks are measured is in the range of 1 to 25 
oersteds, so that in general only the initial susceptibility is measured 
(Stone, personal communication and Collinson et al., 1963). 

The susceptibility meter used for the above samples was calibrated 
for samples of 2.54 cm diameter and 1 cm height. The rock samples were 
cored and sliced into discs of these dimensions. The number of discs 
obtained from each sample varied according to the size, shape, and character 
of the original sample. 

Some of the hand specimens were too friable to withstand coring. 

Among these are: |V~2> V- 1 , V-2, V-3, and VfM. These samples were 
impregnated by immersion in plexiglas dissolved in ethylene dichloride, 
methylene dichloride, or d ] ch 1 oro-ethane . Additional solution was poured 
over these rocks as they dried. This method, described by Noltimeir 
(I 967 ), has the advantage of being an endothermic reaction thus eliminating 
the chance of introducing a thermal component of magnetization. After 
impregnation these samples could be cored and cut into discs of appropriate 
size. 

The unconsolidated samples had settled severely during transportation 
from Katmai to College. Therefore, to insure that any given small portion 
of the sample was representative, of the whole, a splitting technique 
was required. This was accomplished by applying the method discussed by 
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Muller (1967, pp. k2-bk) . The final splits of the unconsolidated samples 
were placed in plexiglas holders of inner diameter 2,5b cm and depth 1 cm 
and weighed using an Ainwortb analytic balance; accuracy - 1.0 mg. Be- 
cause It was not possible to measure densities while in the field, an 
average density has been assumed for conversion of measured bulk suscepti- 
bilities to representat i ve volume susceptibilities. Both Griggs (1922, 
p. 293) and Kienle (1969, p. 128) determined an average density for Valley 
ash of 1.03 grn/cm. Using this density and letting K represent the measured 
bulk susceptibility of a sample of mass M, the standarized volume suscepti- 
bility, K', can be found by: 

K 1 - 5-22 (K/M). 

A representative volume susceptibility for each original field specimen 
was calculated by simply averaging all 2,5b cm x 1.0 cm samples from each 
field specimen. This is a satisfactory technique, except in the cases 
of the samples from sites II and X, which are bimodai. Refer to Tables 
2.1, 2.2, 2.3, and 2 ,b and Figure 2.3. 

2.3 Discussion of Results 

a) Naknek Sediments 

Beyond the northern terminus of the Valley, the Naknek section has 
been sampled by several i nves t i gators (refer to Table 2.2). Although 
the susceptibilities of these samples range considerably, of the series 
collected by Packer (personal communication, 1971) on Mount Katolinat, 
about one-fourth clustered near 20 x 10 ^ emu/cc. The present study has 
also determined an average susceptibility of 20 x 10 ^ emu/cc for altered 
Nakek sediments in the Valley. 

A small aromatic fumarole located on the first Nakek outcrop above 
the Baked Mountain camp is still active; immediately adjacent to the vent 
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Table 2.1 

Susceptibility Studies of Rocks' from the 
Valley of Ten Thousand Smokes 


Susceptibility'' Approx. 

( x 10~6 emu/cc) Total Depth Rock 


Sample 

High 

Low 

Average 

Discs 

(feet) 

Type 

Color 

C onime n ts 

1-1 

406.9 

322.8 

364.5 

8 

0 

P 

G 

surface 
rubb 1 e 

1-2 

8.6 

0.9 

4.9 

4 

0 

1 A 

W 

fumarole cap 

! I-! 

1743. 8 

980.3 

1348.2 

9 

0 

1A 

G** 

Lethe indu- 
rated ash 

1 1-2 

1755.5 

834.6 

1310.8 

7 

0 

1 A 

G** 

same 

1 1 1-1 

1376.0 

1255.0 

1323.8 

2 

2 

UA* 

G 

air-fall from 
Knife Creek 
cut bank 

MI-2 

1177.0 

1107.4 

1142.2 

2 

3 

UA 1 

G 

same 

1 11-3 

1120.1 

1061 .6 

1090.3 

2 

4 

UA* 

G 

same 

111-4 

1461.2 

1303.4 

1382.3 

2 

* 4 

UA 

G 

same 

111-5 

1141.4 

976.7 

1053.3 

2 

6 

UA 

G 

same 

! I 1-6 

1370.6 

1021.3 

1196.0 

2 

6 

UA 

P 

s ame 

111-7 

985.6 

980.6 

983.1 

2 

8 

UA* 

Y 

same 

1 1 1-8 

12 84.2 

1193.8 

1246.8 

2 

8 

UA 

P 

same 

1 1 1-9 

72.7 

65.0- 

63 .O 

2 

9 

UA' 

R 

same 

1 1 1-10 

60.5 

49.5 

55.9 

2 

1 1 

UA 

G 

same 

1 1 Ml 

228.0 

222.4 

225.7 

2 

1 1 

UA' 

G 

same 

MI-12 

1175.7 

1044.4 

1110.1 

2 

1 1 

UA' 

Y 

same 

IV-1 

1.7 

0.0 

0.5 

9 

O’ 

i A 

W 

fumarole cap 

IV-2 

83.9 

37-3 

51.1 

5 

0 

FA 

B 

same 

V-i 

10.5 

7.0 

9.1 

6 

10 

FA 

W 

s t rearn- 
d i ssected 
fumarol e 

V-2 

465- 1 

42.0 

197-2 

7 

10 

FA 

B 

same 

V-3 

92. 1 

0.0 

18.8 

8 

10 

FA 

0 

s ame 

v-4 

1318. 1 

1040.5 

1 172.2 

Z 

3 

UA 1 

0 

nearby 
a i r- f a 1 1 

V-5 

1128.0 

913.9 

1010.2 

4 

6 

UA 

G 

same 

VI-1 

314.6 

305.9 

310.2 

2 

15 

UA 

B 

stream- 
d i ssected 
fumarol e 

VI-2 

624.1 

534.6 

577.1 

2 

15 

UA 

B 

same 

VI 1-1 

23.3 

15.2 

19.5 

12 

1 

FA 

R 

pit fumarole 
in fumarole 
grid 

VI 1-2 

50.2 

39.1 

43-9 

2 

1 

UA' 

Y 

same 

VI 1-3 

37.4 

26.7 

31.7 

2 

2 

UA' 

R 

same 

VII -4 

69.2 

57-6 

64.3 

2 

3 

UA* ' 

B 

same 
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Table 2.1 (continued) 



Suscepti bi 1 i ty* 
( x 10"^ emu/cc) 

Total 

App rox . 
Depth 

Rock 



Sample 

H i gh 

Low 

Average 

Discs 

( feet) 

Type 

Color 

Comments 

V! 1-5 

38.4 

22.2 

31.1 

2 

4 

UA‘ 1 1 

W 

p i t f umarol e 
in fumarole 
grid 

VI 1 1-1 

552 .0 

533 - ^ 

539.7 

4 

0 

IS 

G 

Naknek, 1st 
outcrop 
above BMH 

VI 1 1-2 

22.5 

18.0 

20.1 

5 

0 

! S 

G 

Naknek, 2nd 
outcrop 
above BMH 

Vi 1 1-3 

22.0 

19.3 

20.7 

4 

0 

IS 

G 

Naknek, 3rd 
outcrop 
above BMH 

VI ( 1-4 

31.0 

24.2 

28.5 

3 

0 

IS 

G 

Naknek, 4th 
outcrop 
above BMH 

IX-1 

— 


5058.3 

1 

0 

IS 

R 

h ighly- 
al te red 
Naknek 

X-l 

1650.6 

199.3 

995.6 

17 

0 

L 

G*** 

: Novarupta 
banded 
vol can i cs 

X-2 

793.8 

179.9 

356.it 

9 

0 

L 

G 

Novarupta 
rhyol i te 


Symbols denoting rock type: Symbols denoting color: 


A 

ash 

W 

whi te 

A‘ 

ash with some 1 ap i 1 1 i 

G 

gray 

A 1 ’ 

ash with some c 1 ay 

B 

b 1 ack 

P 

pumi ce 

Y 

ye 1 1 ow 

S 

sed i ments 

0 

orange 

L 

1 ava 

R 

red 

U 

unconsol 1 dated 

P 

pink 

1 

i ndurated (f i rmly) 



F 

friable 



*$us 

ceptibility is in units of 10 ^ emu/cc. 




'-’'Refer to Table 3.5 for details. 
''-*Refer to Table 3.4 for details. 
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Additional Susceptibility Data for 
Rocks in the Katmai Region 

Sample Suscept > b» 1 1 ty* 


Reference 

1 dent i f i cat ion 

High 

Low 

Average 

Rock Type 

Location 

Chantry-Pr i ce 

1—3 


3 

0 

3 

sed iment 

near Overlook 

B i ngham 

K 1—K 6 


330 

190 

171 

sed iment 

along Margot Creek 

Packer 

KTM 1 1 2 — KTM 

145 

340 

17 

94 

sed iment 

Mt . Katol i nat 

Stone 

GR 1--GR 8 


2400 

280 

1200** 

a ndes i te 

Mt. Griggs 

S tone' 

KR 1 — KR 23 


3200 

1900 

2474 

andes i te 

Mt. Katmai crater rim 

Trible 

KT 1 1 


— 

— 

3500 

andes i te 

Mt. Trident 1953 flow 

Tr i b 1 e 

Trident Cone 


— 

— 

1933 

andes i te 

T r 1 dent Cone 

Trible 

KFM 1 — KFM 2 


2500 

2250 

2375 

andes i te 

Fa f 1 i ng M tn . 

5 tone 

S 15—S 18 


3000 

1 400 

2300 

andes i te 

Cerberus Dome 

Stone 

M 1--M 6 


3800 

3400 

3567 

andes i te 

Mt. Mageik 

Trible 

KN 3-1-65 
KM 3-2-65 
KN 1-3-65 


330 

230 

287 

rhyol i te 

Nova rupta centra 1 dome 

Trible 

KN 4-2-65 


— 

— 

1180 

banded lava 

Nova rup ta east margin 

Trible 

NRC 1 — N PX 3 


1377 

704 

1084** 

banded pumi ce 

Novarupta crater 

Trible 

NOVA- 4 


— 

— 

2200 

tuff 

Novarupta Basin fumarole crater 

S tone 

S 10--S 14 


1600 

1200 

1400 

tuff 

Fissure Lake 

Stone 

S 1— S 9 


1200 

570 

826 

tuff 

Corner Lake, upper unit 

Trible 

S 0 



— 

1230 

tuff 

Corner Lake, lower unit 

Trible 

LR 1 


— 

— 

720 

tuff 

Lethe River 

S tone 

KA--KI 


1320 

720 

367 

tuff 

Middle Valley 

Stone 

J 1— J 3 


753 

500 

658 

tuff 

Juhle Creek 

Trible 

KPFP 2-70 


— 

— 

1025 

tuff 

Peckish Springs 

Packer 

PFK 1--PFK 6 


2600 

740 

1583 

tuff 

Upper Kni fe Creek 

Trible 

BMV 1 1--8MVI 

2 

1420 

1007 

1214 

intrusive 

Broken Mtn. Valley 


''Susceptibility is in units of 10 ^ emu/cc. "^Apparently bi modal . 


\_n 

vn 
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the ground Is perceptibly warm and mosses are present, although no 

emanation is visible. Sample VIM taken from this outcrop exhibits a 

-6 

susceptibility in excess of 500 x 10 emu/cc. It is assumed that ex- 
posure of elevated temperatures, as well as to fumarolic emanations, led 
to chemical alternation of this sedimentary section resulting in an enhance- 
ment of its magnetic susceptibility. 

There have been several reports of Naknek sediments which were obviously 
altered as a result of the 1912 volcanic eruptions; Griggs (1922, pp. 243“ 
244) described dull brick-red Naknek exposures in the Baked-Broken Mountain 
area. He attributed their "baked" look to chemical change by excessive 
heat and perhaps other agents. In 1923 when Fenner (1950a, p. 607) descended 
into Katmai crater (the lake was reduced to a small lagoon within the horse- 
shoe island), among the rock fragments on the floor he found "much shale, 
some of it having a ’.baked 1 look." Fenner (1923> P- 35) also reports the 
presence of many flat bits of reddish shale and sandstone throughout the 
tuff-flow; he describes them as follows; 

"The bits of shale and sandstone in the tuff bear 
evidence of having been exposed to heat. The color of the 
Naknek sediments which underlie the Valley and form the hills 
and mountains outside of the volcanic range, is generally of 
a somber green and green i sh-gray . The fragments included in 
the tuff are almost always reddened or blackened if of shale, 
and reddened or whitened if of sandstone, 11 

Only one sample of such highly altered Naknek sediments was collected. 

This sample is streaked with red and yellow throughout, although the 

stains are darkest on the outside. This sample was part of the surface 

rubble of the Broken Mountain Valley slope of Greasy Pass. It possesses 

a magnetic susceptibility of 5058.3 x 10 ^ emu/cc. 

Zies (1929) points out that the volcanoes of the Katmai area have 

their orifices in sedimentary rocks, so that it is reasonable to expect 
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that metamorph i c changes will have occurred at the contacts of the two 
types of rock. Zies also considers that the heat from the large igneous 
mass represented by Novarupta would be quite effective in removing the 
volatiles from adjacent sedimentary rocks. One could therefore expect to 
find extensive zones of altered Naknek rocks (which could posses high sus- 
ceptibilities) adjacent the eruptive conduits throughout the valley region. 

The preceding descriptions suggest that alternation is commonplace 
among Naknek rocks In the Valley, where surfaces of fragments have been 
subjected to elevated temperatures and vapor-phase contamination. For 
the model studies, a magnetic susceptibility of 5000 x 10. emu/cc (sample 
1 X - ( ) will be associated with sedimentary rocks of this history, 
b) Lavas 

The volcanics bordering the head of the Valley are predominately 

lavas of intermediate composition. The similarity of these rocks is 

apparent from the andesitic compositions of several samples from the Valley 

area (Fenner, 1926, pp. 676 - 679 .; Ray, 1967» pp. 1**1, !**5)- 

Forbes et al. ( 1 969 , p. 118) found that all the andesities erupted 

by Mt. Trident throughout the recent years are similar, and furthermore 

that "all andesites erupted by Alaskan orogenic volcanoes in continental 

settings from 1912 to the present are highly siliceous and remarkably 

similar in bulk composition." They propose that this material Is generated 

by the anatexis of 1 owe r crustal material. 

Table 2.2 contains measured susceptibilities for several samples 

of andesite from the volcanoes bordering the head of the Valley. Their 

average susceptibility is 2*0! x 10 ^ emu/cc, although ‘one sample ranged 

-6 

as high as 3^00 x 10 11 emu/cc. For modeling purposes, a value of 2500 x 
“6 

10 emu/cc can be used. 
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Acidic volcanics were erupted in 1912, primarily from Novarupta; 
however, Fenner (1926, p. 676 ) sampled a rhyolitic boulder on the southern 
rim of Katmai crater. The origin of the rhyolite lava has been much de- 
bated. A complicating factor is that a considerable amount of the ejecta 
attributed to Novarupta, as well as the margin of the dome itself, exhibits 
a banded structure (refer to section 1.4b). According to Fenner (1923, 
p. 56 ), a great part of the banding is due to streaks of dark-brown or 
nearly black scoria which often contain great quantities of phenocrysts 
within a mass of light-gray glass almost without phenocrysts. Several 
samples of the Novarupta volcanics have been analyzed (Fenner, 1923* P- 57; 
Zies, 1929, p. 56; and Forbes, personal communication, 1971) • Host of the 
investigators have identified the dark bands as andesitic (Griggs, 1922, 
p. 297; Fenner, 1950b, pp. 707 - 7 1 0 ; Williams, 1954, pp. 58-59; and Curtis, 
1968, p. 194). The analyses of Forbes et al., (1969, p. 118), however, 
suggest that the basic inclusions in the mixed i ava are basaltic. 

The central dome of Novarupta is almost pure rhyolite. The average 
susceptibility of the rhyolitic lava is 304.1 x 10 ^ emu/cc. The banded 
lava of the margin of the dome and the banded pumice found In the moat 
(or crater) is of higher susceptibility. This is to be expected since 
the dark bands contain a greater percentage of ferr i -magnet i c minerals 
than the light mass. An analysis by Fenner (1923, P- 57) shows that the 
light bands consist of 0.82% Fe^ and 1.43% FeO; whereas the dark bands 
are 3.4% Fe^O^ an< ^ ^*53% FeO, This dual nature of the banded rocks results 
in a bi modal distribution of susceptibilities for Novarupta and its banded 
ejecta; this distribution is seen in sample X~1 (see Table 2.3). The aver- 
age susceptibility of the banded lava and pumice is IO 85 x 10 ^ emu/cc. 

,p '6 

For modeling purposes, a value of 1250 x 10 emu/cc may be used for this 



Table 2.3 


Bimodal Susceptibility of Novarupta Volcanics 


i b i 1 i ty 


X-1 


1 e Disc 

(x 10'° emu/cc)* 

Dark bands 

1 

1283 . AO 

yes 

2 

230.80 

no 

3 

199.33 

no 

A 

1300.88 

yes 

5 

1230.34 

yes 

6 

1650.58 

yes 

7 

1265.91 

yes 

8 

1528.19 

yes 

9 

283.26 

no 

10 

290.25 

no 

11 

272.26 

no 

12 

1592.30 

yes 

13 

226. )A 

no 

1A 

1 A 5 S .25 

yes 

15 

1353.79 

yes 

16 

1 405 - 79 

yes 

17 

1 353-34 

yes 

usceptibi 1 i ty i ; 

in units of 10 ^ emu/ cc. 
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mixed material; whereas a value of 250 x 10 emu/cc wi 1 1 be used for the 
rhyo 1 i t i c mote rial. 

Although none of the banded pumice on the Valley floor was sampled, 
a piece of light gray pumice about 5 inches in diameter was obtained 
(sample l-l). This gray pumice is representative of much of the surface 
rubble in the middle and lower Valley. Its susceptibility of 364. 5 x 10 
is close to the value for the Novarupta rhyolite from which this pumice 
was no doubt derived; apparently it was erupted from Novarupta while the 
conduit was supplying only rhyolite. 

c) Glacial Deposits 

In traveling the Katmai Trail during a geological reconna t ssance of 
southwestern Alaska in 1893, Spurr (1900) traversed the valley which was 
to become the Valley of Ten Thousand Smokes. He describes the valley. as 

\ 

filled with about 30 m of glacial drift consisting of stratified gravels, 
sand, and near Katmai Pass, boulders, too. He also reports that fragments 
of the Naknek sediments were included in the valley drift. It is difficult 
to arrive at a susceptibility representat i ve of this glacial material. It 
•is reasonable to assume that this detritus is predominately composed of 
fragments of the basic lavas from the volcanic peaks at the head of the 
Valley which carry many glaciers. The portion of fragments of Naknek sedi- 
ments from the Valley walls included in the drift probably varies locally 
but is always small. If the drift were composed entirely of volcanic 

material, the highest susceptibility it would possess would be that of the 

-6 r 

parent lava in situ; 3800 x 10 emu/cc is the largest measured susceptibility 
for lava. At the other extreme. If the bulk of the drift were from Naknek 
sediments, its susceptibility could be neglible. Obviously, it is impossible 
to predict the suscep t i b i 1 i ty of the glacial drift lying between the pyroclastic 
flow and the Naknek bedrock. A further complication is that the constituents 
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of the drift would be expected to be in random orientation, thus remanence 
would alter the apparent suscept i b l 1 i ty of the material were it analyzed 
under the assumption of magnetization under the present field. 

The situation is further complicated by the lack of knowledge of 
the temperature regime during the 1912 eruption. A maximum fumarolic 
temperature of 6A5 degrees C. was measured in 1919; this occurred at 
fumarole No. 153 which is located in the middle Valley (Allen and Zies, 
1923) p. 10b). Lovering deduced that fumarole No. 1, located at the 
terminus of the ash-flow* had an initial temperature between 800° C and 
900° C (Lovering, 1957, P» 1590). These data imply that the Implacement 
temperature of the pyroclastics was well over 600°C. 

According to Nagata, the Curie temperature of igneous rocks can be 
estimated If the ratio Fe^O^/ (Fe^O^ - FeO) is known. An average ratio 
of 0.k5 is representative of the igneous rocks at the head of the Valley 
(samples 575, 568, 526, $83, and 1A7 of Fenner, 1926, pp. 676 and 682; 
sample 153b of Fenner, 1923, p. 57) . This ratio corresponds to Curie 
temperatures in the realm of k00°C to 600°C (Nagata, 1961, p. 138). It 
therefore seems probable that some realignment of magnet ) zat ion occurred 
within the fragments in the upper portion of the glacial fill of the 
Valley upon emplacement of the ash. Thus, the apparent susceptibility 
of this limited zone of reheated drift will increase slightly. Also, 
such a high thermal regime could result in compositional alterations 
within the drift. It was found that after exposure to extreme temper- 
atures and/or volcanic fumes that Naknek sediments could attain suscepti- 
bilities as high as 5000 x 10 ^ emu/cc. However, Curtis (1968, pp. 1 8 -4 ~ 
185) reports that of the several examples throughout the Valley, of 
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morainal boulders of andesite in direct contact with the tuff flow, none 
show the slightest sign of alteration. 

In the foregoing discussion, it is clear that a singular suscepti- 
bility cannot be assigned to the glacial drift lying between the 1912 
pyroclastics and the Naknek basement. Depending on the composition, 
proportion of components, and history of this material, its susceptibil- 
i ty could conceivably range from 10 to 10 " cgs units. This indeter- 
minate suscept i b i 1 i ty must be considered in any geological i nterpretat i on 
of model studies. 

d) Pyroclastic Flow 

The top of an indurated deposit Is exposed where traverse C-C' 
crosses the River Lethe. In the field the rock Is deep gray, but upon 
drying it assumes a lighter color with prominent bands of brown and white. 
Preliminary petrographic studies show this rock to be vi troc] ast i c , and 
not welded as was first assumed. A few flow structures are also present 
(Forbes, personal communication, ■ 1971) « 

The average susceptibility for the Lethe indurated ash (samples If"? 
and 11-2) is 1330 x 10 ^ emu/cc. The susceptibilities of Individual 
samples vary from 834.6 x 10 ^ emu/cc to 1755.5 x 10 ^ emu/cc. Examinat- 
ion of the character of the tuff comprising the individual discs disclosed 
the correlation of colored streaks with measured susceptibility; white 
and brown streaks are prominent in discs of lower suscept I b i 1 i ty , black 
bands are associated with the higher vaiues (see Table 2.4). 

An exposure of similar tuff was observed along the River Lethe above 
the falls which, are in the area of sample 1-1. Other outcrops 



Table 2 .4 


Susceptibility Studies of Lethe Indurated Ash 


Sample 


Measuri ng 
di sc 


Susceptibility General 
(x 10' 6 emu/cc)* description 


1 

1318.37 

black 

2 

980.33 

b rown 

3 

1318.37 

black, brown 

4 

1085.23 

b 1 ack 

5 

1 108.55 

b rown 

6 

I545.b7 

brown, black 

7 

1510.70 

b 1 ack 

8 

1522.36 

b ] ack 

9 

1743.84 

black 

1 

1493.22 

black, wh 1 te 

2 

1458.25 

black 

3 

1195.97 

black, brown, white 

4 

1755.49 

black 

5 

1242.60 

black and wh 1 te 

6 

834.62 

black, white prominent 

7 

1195.97 

b 1 ack , wh 1 te 


Suscept i bi ]i ty is 


i n 


units of 10 


-6 


emu/cc. 
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resembling the Lethe indurated tuff are to be found in connenction with 
explosion craters in Novarupta basin and at the base of Mt. Griggs; at 
the base of Mt . Katmai where the indurated exposure extends beneath the 

Knife Creek Glaciers; along faults at the base of Mt. Mageik by Fissure 

\ 

Lake where the indurated section is overlain by glacial fill through 

which fumarole vents can be traced, and by Corner Lake where two tuff 

sections are exposed; at the southwest base of Baked Mountain active 

fumaroles expose an indurated section beneath the ash (Forbes, personal 

commun icat ion , 1971)* Subsequent to the present study, several of these 

other tuffaceous outcrops have been sampled. Although Individual samples 

-6 , -6 

possess susceptibilities ranging from 500 x 10 emu/cc to 2600 x 10 

emu/cc, in general, the susceptibilities of the tuff throughout the 

-6 

Valley appear fairly uniform at .its average of 1187 x 10 emu/cc (refer 
to Tables 2.1 and 2.2). In conjunction with the strong visual similarities 
of these rocks, their common susceptibilities lend support to a hypothesis 
of singular origin for this indurated materia) throughout the Valley. 

In the lower Valley the deposits possess a definite pink cast and 
appear to be homogeneous and indurated. The composition of Lovering's 
(1957, p. 1593 ) sample ( 7TL53 ) of normal ash near the terminus of the 
flow is remarkably similar to Forbes’ (personal communication, 1971) sample 
(KN 2-3-65) of Novarupta rhyodacite. Probably, the entire flow here Is 
more rhyolitic than farther up the Valley, where the samples for this study 
were obtained. Unfortunately, no susceptibilities have been determined 
for the flow in the lower Valley. It can only be assumed that the suscep- 
tibility of this portion of the ash flow is similar to that of the rhyolitic 
lava of Novarupta and the uniform pumice (sample 1-1); i.e., approximately 
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-6 

300 x 10 emu/cc. 

Probably the tuffs are the product of partial welding associated 
with the cooling, degassing, and compaction of the pyroclastic flow of 
1912. The variables controlling the presence, extent, 'and degree of 
welding within a given ash flow are discussed in detail by Smith (i960) . 
Briefly, an ash flow of sufficient emplacement temperature and/or thick- 
ness will form a central lens of indurated material of thickness and quality 
dependent upon the rates of welding, cooling, and crysta 1 11 zat i on pecul far 
to that particular deposit. 

e) Air-fall Pyroclastics 

Air-fall pyroclastics were collected 'from approximate depths of 
1 and 2 meters in a stream cut in the eastern drainage of Broken Mountain 
Valley. The top 0.3 meter of the ash in this area appears to be reworked 
and was not sampled. - 

The re-deposited nature of the upper ash is implied by the manner 
in which it is draped over the gentler sloping banks of the stream. 

This mantle of grayish ash was probably originally deposited upstream and 
on the adjacent mountain slopes. Transportation of air-fall pyroclastics 
from these slopes was recognized by Forbes who reported that over the 
years great quantities of ash and pumice have been transported down the 
valley sides by spring melt waters, in the form of a slurry-like mass of 
suspended ash and ice (Forbes, personal communication, 1 968) . 

The next layer consists of about 1 meter of orange-stained ash with 
a few lapilli. This is underlain by at least 2 meters of gray ash. 

These two layers were sampled, and probably represent air-fall. No attempt 
has been made to correlate these layers with the tephra section of Curtis. 
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More air-fall pyroclastics were collected In the southeastern branch 
of the Valley. Knife Creek Is braided near the crossing of traverse 
B-B ] and the valley is quite flat to the east. Not far above the tra- 
verse crossing. Knife Creek narrows, and a k meter cut bank rises abruptly 
on the western edge of the stream; it extends for about a kilometer down- 
stream. Plate ! shows the multi-colored layers of air-fall exposed by the 
stream cut. The layering of this cut bank should represent the air-fall 
History of this part of the Valley; however, it is difficult to correlate 
this sequence with the published sections of Fenner (1923, Table 1; and 
1950, Table 2 and pp. 712-71 and Curtis (1.968, pp. 162-183 and pp. 196- 

201). Tentative assignment of Curtis* system to the Knife Creek section 
is shown in Figure 2.2. 

Discrepancies can be attributed to the fact that very few beds at any 
given outcrop have correlative value for more than a few hundred meters. 

The many changing conditions during an eruption such as fluctuation in 
intensity of eruption, the wind direction and speed, an d the amount of 
turbulence within the rising clouds of tephra, can account for such 
spatial variation (Curtis, 1 968 , p. 167 ). 

Samples 111-9, -10, and -11 exhibit relatively low susceptibilities 
which have similar values to those assigned to fumaroiic surface deposits. 
Inclusion of these values among the fumaroiic susceptibilities only raises 
the average to 106 . 9 x 10 ^ emu/cc, or omitting VI-1, VI-2, and 111-11, 

45.9 X 10 6 emu/cc (see section 2.3f). One is ied to the conclusion that 
following deposition of this part of the Knife Creek section, there was a 
period of little or no air-fall during which fumaroiic venting occurred 
from the ash en masse. The occasional overlapping of the red and white 
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Figure 2.2. Comparative tephra sections based on Curtis ( 1 968 ) and the section exposed 
at sample site III. Plate I is a photograph of this same site. Magnetic 
suscept 1 b i 1 i t fes for these samples are given in Table 2.1. 
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layers at this depth could be indicative of fumarolic activity or of 
resorting by wind. The red zone might well denote fumarolic staining. 
Fenner (1950b, p. 712) concludes that the yellow and crimson colors at 

the top of certain ash layers is probably caused by gas emanations. it 

is also of interest to note that the zone of low susceptibility is centered 

about a thin white layer which has been tentatively correlated with 
Curtis' layer E. From his tephra study, Curtis concluded that "following 
the deposition of layer D was a period of relative quiescence. A thin 
bed of ash... termed layer E, was deposited." This period of quiescence 
could relate to a time of fumarolic venting as indicated by the low 
susceptibility ash at this depth. Curtis also deduced that at the time of 
the eruption of layer F, "confl icting winds must have been blowing in the 
vicinity of the headwaters of Knife Creek" (Curtis, 1968, p. 198). This 
condition could account for the overlapping at the top of the low suscepti- 
bility tephra. 

The tephra section at site 111 contains several zones and splotches 
of red and yellow. The manner of such staining Is discussed by Fenner 
( 1926 , p. 741): 

"In thick beds of ejected pumice certain strata are deeply 
colored -- crimson, yellow, or yellow-brown. This is due to 
mere staining of the pumice with iron oxides, but the lumps are 
deeply impregnated. Presumably the iron was deposited by the 
little residual gas that oozed out subsequent to the first violent 
inflation of the pumice. It is remarkable that some strata show 
this staining and others, above and below, are free from it." 

Since the anomalous low susceptibilities probably represent altered 

ash, they are omitted in arriving at an average susceptibility for air- 

fall of 1155.5 x 10 emu/cc. 
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Table 2.5 


Magnetic Susceptibilities of Some Common Constituents 
of Fumarolic Alterations Reported by 
Shipley (1920), Zies (1929), and Lovering (1957) 


Material 

(Susceptibility x 10“^ 

sulfur 

-15.5 to 700.0* 

corandum 

-37-0* 

silica 

-29-6* 

opal 

diamagnetic 

kaol in i te 

diamagnetic 

mon tmori 1 1 on i te 

d i amagnet i c 

apat i te 

d i amagnet i c 

halite 

-30.3* 

ammonium chloride 

-36.7* 

lead chloride 

-73.8* 

arsenic sulfide 

-70.0* 

spha ler? te 

-25.0* ‘ 

galena 

-84.0* 

gypsum 

-74.0* 

a 1 um 

d i amagnet i c 

lead sulfate 

-697-* 

ba r ? um su 1 fate 

-71 .3* 

Mate rial 

Susceptibility emu/cc 

py r i te 

.000005 to .0002** 

i ron chi or i des 

.00998 to .01475* 

i ron sul fates 

.0102* 

FeO 

.0072* 

hemat i te 

.00004 to .0001** 

magnet E te 

.04 to 2.0** 

magnetite crystals 

6.3 to 24.0** 

i Irnen i te^ 

.03 to .14** 

f rankl i n i te ^ 

.036** 

py r rhot i te?- 

.007 to .028** 

specu 1 a r i te" 

.003 to .004** 


^ 2 i es reported traces of Zn , Ti, & Mn in fumarolic magnetite. 

^-These minerals were probably present, 

*CRC Handbook of Chemistry and Physics, 1965" 1966, P* £"95 to E-100. 

Jakosky , 1950 , p. 164 and 165- f 
***Suscept i b I 1 i ty is in units of 10“° emu/cc. 
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indicated by sharp, high, narrow magnetic anomalies which show excellent 
correlation with surface markings and fumarolic lineations in the fuma- 
role grid in Broken Mountain Valley. This is discussed fuithei in Chapter 
III, Sharp anomalies were observed near fumaroles along magnetometer 
traverses across the Valley floors as well. Unfortunately no samples 
were taken from below the surface in connection with these anomalous 
highs in the magnetic field. 

Since only exposed sections of fumaroles were sampled, the study 
was regrettably limited to those areas which would have experienced the 
most leaching and weathering. The exposure of these fumaroles, as well 
as the diamagnetic nature of many of the i'ncrust ions , account for the low 
susceptibilities encountered in this study. The presence or dominance of 
hematite in other fumarolic samples is probably responsible for the relat- 
ively higher susceptibilities found. Those values greater than 100 x 10 
emu/cc probably indicate the presence of some magnetic minerals. 

2.k Conclusions 

Several of the rock types have indistinguishable susceptibilities 
(refer to Figure 2.3). Both Naknek sediments and exposed fumarol i c ash 
can be represented by a susceptibility of 50 x 10 emu/cc. There is 
considerable overlapping of susceptibilities for air-fall pyrocl ast ics , 
tuff and the banded volcanics. All three have susceptibilities near 
1250 x 10 emu/cc. Glacial drift may well have an apparent suscepti- 
bility in this range, also. The range of susceptibility of many of the 
other rock types overlap,, too. 

The clustering of susceptibilities near 1250 x 10 ^ emu/cc could 
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Table 2 .6 


Summary of Susceptibility Data 



Suscept i b i 1 i ty 
Measured Measured 

(x 10' 6 

emu/ cc) * 

Material 

High 

Low 

Average 

Mode 1 i ng 

Naknek sediments 

340 

0.3 

94 

50 • 

altered Naknek 

— 

— 

5058 

5000 

andesitic lava 

3800 

1 400 

2431 

2500 

mixed lava 

1377 

704 

1085 

1250 

rhyol i t i c 1 ava 

356 

230 

304 

250 

glacial d r « ft 

— 

— 

— 

1250? 

tuff 

2600 

500 

1187 

1250 

tephra 

•1461 

913 

1156 

1250 

furnarolic alterations 

624 

0 

107 

50 


*Suscept i b i I i ty i s in units of 10 ^ emu/ c c. 
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be Indicative of genetic similarity i n the materials. Novarupta has been 
suggested as the major eruptive vent for the tuff and most of the air-fall. 
The similar susceptibilities of tuff, air-fall and Novarupta's banded 
ejecta certainly bear out this hypothesis. 

The susceptibility data is summarized in Table 2.6; also given are 
the modeling susceptibilities which will represent each rock type in the 
model studies in Chapter IV. Further discussion of susceptibility as 
related to geological interpretation in the Valley can be found in that 
chapter. it appears that there is sufficient suscept ib i 1 i ty data to test 
some of the previous hypotheses concerning the geologic composition of 
the Valley fill. 



CHAPTER HI 


SMALL SCALE MAGNETOMETER SURVEY OVER 
A ZONE OF RELICT FUMAROLES 

3.1 Small-scale Magnetometer Survey 

A distinctive feature of Broken Mountain Valley is the colorful 
plaid of orange and red over tan, produced by intersecting lineat 
zones of fumarolic markings. These bands of remnant fumaroles trend 
approximately parallel and perpendicular to the axis of the valley. 

Near the terminus of the valley, an intersection of two of these 
fumarole lines was chosen for detailed magnetomet ry . A five by ten 
foot ( 1.52 x 3.04 meter) rectangular grid was laid out using a 100 
foot (30.5 m) cloth tape, and covered an area extending 250 feet 
( 76.3 m) along magnetic north and 200 feet (61 m) wide. The location 
of this fumarole grid is shown in Figure 4.1. 

A detailed sketch of the surface features of the fumarole grid 
is shown in Figure 3.1a, and Plate II is a view from magnetic south- 
east. The fumarole remnants delineated by bright splotches of fumarol ic 
clay and discolored ash downslope, stand out sharply against the plain 
light ash. A three (3) meter deep crater vent with steep sides of 
loose ly-coherent , red-stained ash is located at R-32. There is also 
an elongated trough in the unaltered ash; this is probably the result 
of collapse associated with the. fumarole line which parallels it to 
magnetic west. At each end of this trough are minor pits, one of 
which Is almost two (2) meters deep. Excluding these depressions, 
the area of the fumarole grid is relatively flat. 
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Figure 3- la. Surficial geology of the fumarole study area. Contour interval 
is 5 feet. Map units designate the color and composition of the 
surficial deposits. Plate II is a photograph of this same area. 

Figure 3- lb. Total field magnetic anomaly map over the fumarole study area. 

Stations 0“50 were positioned at 3 foot intervals along traverses 
A-V, which were 10 feet apart.- Solid contour interval is 100 y 
with dashed contours at 50 y intervals. The cor respondence of the • 
large magnetic anomalies with the surficial geology (Fig. 3.1a) 
suggests that fumarolic magnetic minerals may be preserved at 
shallow depths. 
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The total magnetic field was measured at five foot (1.52 m) 
intervals along traverses A through V, which are separated by ten 
feet (3.04 m) . The probe height was 2.07 meters. Since readings 
were taken in rapid succession, repeated readings generally were 
not made. The magnetometer readings were reduced to anomalies in 
the total field according to the procedure outlined in Appendix A. 
Figure 3.2 is a computer drawn perspective view of these anomalies 
from magnetic northeast. Contours of these anomalies are presented 
in Figure 3. lb. 

Only two attempts were made to measure vertical gradients with- 
in the grid. At T~35 which is situated just outside the pit fumarole, 
the vertical gradient of the total field is approximately 200 gammas 
per meter. At M- 36 , which is located near a 700 gamma anomaly in a 
zone of orange clay and altered ash forming part of a down-valley 
striking fumarole line, the magnetometer failed to record a steady 
reading in the lower .405 meter probe posi tion. Since the proton 
precession signal decays too rapidly for the magnetometer to function 
in fields with gradients greater than about 650 gammas per meter 
(Hood, 1965, p. 404) , the close proximity of a disturbing body is 
I ndi cated . 

3.2 Analysis of the Magnetometer Data 

Steep horizontal gradients of 75 to 200 gammas per rneter indicate 
that the magnetic bodies responsible for the magnetic anomalies 
observed in the fumarole grid are located at shallow depths. The 
high vertical gradients reported in the previous section also indicate 
near-surface magnetic sources. 
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Figure 3.2. Computer drawn perspective magnetic anomaly map of the 
fumarole study area. Refer to contoured anomaly in 
Figure 3.1b. 


Gammas 
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As a rule of thumb, '.'e width of a magnetic body cannot exceed 
the width of its anomaly. Thus, the bodies causing the anomalies 
in the grid are necessarily narrow. 

The remarkable agreement between the trends of the magnetic 
anomalies and the linear surface expression of former fumaroltc 
activity is apparent upon comparison of Figures 3»ia and 3.1b. The 
possibility of deposition and retention of magnetic minerals along 
fissure vents was presented in section l.Ad. The probability of 
occurrence of fumarolic magnetite in the area of the survey grid 
will be discussed in the next section 

Estimates of the depth to the top of the magnetic bodies caus- 
ing anomalies in the grid were made by applying Peter's "slope 11 
method to profiles of the total residual field (Dohr in, I960, pp . 
312-313). This technique assumes that the source is a thin vertically- 
magnetized, vertical dike of infinite length. Also, the width should 
be of the same order of magnitude as the depth and the thickness 
should be much greater than the depth of burial. Although the formulae 
are for vertical anomalies, they should be equally applicable for total 
anomalies at such a high magnetic latitude as the Valley. The high 
gradients and narrowness of the anomalies are clear evidence that the 
source of the anomalies is indeed thin, shallow, and linear, but the 
thickness of these accumulations has not been resolved. Application 
of this technique to individual profiles within the grid indicates 
that the depth to the magnetic bodies varies from one to three meters; 
deepest below the crater vent and shallowest near the 700 gamma 
anomaly of L - 3 7 . 
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3.3 Implied Source for the Anomalies 

Although It was never the site of extensive studies, Broken 
Mountain Valley did receive some attention by the early investigators. 
In describing Broken Mountain Valley ,• Gri ggs (1922, p. 233) reports: 
"Its mouth is the seat of some of the most vigorous vents in the 
region, and smaller fumaroles occur farther up." Fenner {1923, p. 15) 
was also impressed by the amount of fumarolic activity in this valley. 
In 1919, Allen and Zies (1923, pp. 104-105) examined the emanations 
of ten fumaroles in Broken Mountain Valley, primarily at its terminus. 
Although they found the valley to be one of the hottest areas in the 
Valley of Ten Thousand Smokes, no exceptionally high temperatures 
were encountered in 1919 ; measured temperatures for the. fumarol i c 
gasses ranged from 94°C to 353°C. Examination of Allen and Zies' 
map (Allen and Zies, 1923, p. 80 ) indicates that the grid area was 
not the site of any of their fumarole studies. However, their studies 
in other regions of the Valley of Ten Thousand Smokes imply that such 
an area could have been the site of the subl imination of magnetic 
minerals. They report that the hottest fumaroles were small circular 
crater vents, "say fifteen feet in diameter and ten feet deep" (Allen 
and Zies, 1923, p. 91). They also found the throats of the hotter 
fumaroles lined with a crust of "metamorphosed pumice conspicuously 
colored by oxide of iron. Some of the oxide is bright red and of 
loose texture" (Allen and Zies, 1923, p. 97)- Similar features still 
possessed high temperatures in 1 9 19 ; temperatures from 400°C to 645°C 
were measured for a group of crater vents in the main Valley (Allen 
and Zies, 1923, p. 91). The crater at R~32 appears similar to those 
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studied above, so it can be assumed that this vent was probably a site 
for magnetic mineral accumulation. 

Although crater vents often served as the principal outlets along 
a fissure line, sometimes the fissure itself was a line vent of considei 
able length. Even in 1 9 19 , temperatures above ^00°C were recorded along 
some fissure vents (Allen and Zies, 1923, p. 91). it is thus possible 
that the rectilinear features of Broken Mountain Valley are filled 
fissures which could have been the site of sublimation for magnetic 
minerals. 

In the early years after the eruption, Griggs reported that the 
fissures criss-crossing the center of the valleys were mostly filled 
and "usually marked only by the lines of encrustations or of small 
fumaroles stretched along their courses" (Griggs, 1922, p. 235). This 
early description is applicable to the present situation in Broken 
Mountain Valley. No open fissures were observed at this location; 
primarily the fissure pattern was revealed by lines of fumarolic clay 
and discolored ash. It seems reasonable that these centYal fissure 
fumaroles were closed even early in their history. In this case, 
metallic sublimates along their vents may have been protected from 
leaching by acid condensates in the latter stages of fumarolic 
activity and from subsequent weathering. It is proposed that the 
magnetic sublimates thus preserved along the fissure vents at shallow 
depths are responsible for the narrow, sharp, high magnetic anomalies 
encountered over the linear surface markings indicating earlier 


f uma ro 1 i c a 1 te ra t i on . 
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Broken Mountain Valley is much higher than the other branches of 
the Valley. The pyroclastic flow in this valley is also the thinnest 
in the region (Sbar and Matumoto, 1971)- One would thus conclude 
that the top of the Naknek in this valley is also higher than in the 
other valleys. As was proposed for some of the fissures in the 
southern branch of the valley, settling over bedrock topographic 
features can result in fracturing of the pyroclastic flow. In this 
manner, parallel fractures could have been produced in the pyroclastics 
near the terminus of Broken Mountain Valley where the deposit settled 
over a sharp break in elevation of the underlying Naknek. Contracting 
and settling within the valley confines would then explain the inter- 
secting system of fissures. 

The fissure fumaroles within the grid area intersect at an angle 
of about seventy degrees. Similar alignment of fumarolic fissures 
Is reported by Sheridan for the Bishop Tuff. He found that the fuma- 
rolic fractures have orthogonal intersections; for most, the dihedral 
angle is about sixty degrees. These joints penetrate up to 120 meters 
into the sheet, well within the densely welded, devitrlfied zone. 
Sheridan's studies indicate that the fracturing mechanism was complex, 
related to welding deformation as well as to thermal stress release, 
and somewhat influenced by the underlying topography (Sheridan, 1970, 

pp . 86 0 , 86 1 ) . 

The magnetometer survey within the grid offers no proof of the 
extent of the fumaroles in this area but is a very strong indication. 
All that is demonstrated is near-surface , highly-magnet i c accumulations 
along the fumaroles. Nothing from deeper horizons can be inferred. 
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3.4 Conclusions 

This study has shown that narrow magnetic anomalies encountered 
in the vicinity of fumaroles can be directly related and probably 
are genetically related to the fumaroles. In the fumarole grid, the 
sharp anomalies are restricted along surficial zones of fumarolic 
alteration. Analysis of profiles over these fumarolic markings 
indicates the sources are necessarily shallow, narrow, and elongate. 
Previous studies Imply that this area was one of extreme temperatures 
following the emplacement of the 1912 pyrocl ast i cs , thus, it was 
probably the site of sublimation of magnetic minerals. The existence 
of the anomalies implies that these fissure fumaroles were covered 
early in their history, thereby protecting and preserving the magnetic 
accumulations along their vents. 

It Is regretable that samples were not collected at various 
depths within the fumarole grid. Samples from a few meters depth 
below magnetic highs would be most informative. As for now we can 
only speculate as to the composition and extent of the accumulations 
responsible for the sharp anomalies coincident with the surface 


fumarolic alterations. 



CHAPTER IV 


CROSS -VALLEY MAGNETIC PROFILES IN THE VALLEY 
OF TEN THOUSAND SMOKES 

4,1 Magnetometer Survey in the Valley Region 

Magnetometer traverses were made across each branch of the Valley 
of Ten Thousand Smokes (refer to Figure 4.1). it was hoped that this 
restricted survey could clarify some of the controversy about the 
structure and composition of the ash flow and its confining pre-1912 
valley. Appendix B includes the details of these magneotme ter tra- 
verses and the listings of the reduced data. Profiles of the various 
valley crossings are given in Figures 4.2 through 4,5. 

The jagged nature of the valley profiles makes analysis difficult. 

As discussed in Chapter III, it is probable that these narrow anomalies 
are caused by small, near-surface pockets of magnetic minerals asso- 
ciated with fumarollc alterations. It is unfortunate that the present 
survey was a ground survey and thus particularly sensitive to the effects 
of these small, close features. Figure 4.6 gives two of Anma's (1972) 
aeromagnetic profiles which are near ground magnetometer traverses. The 
smoothed nature of the aeromagnetic profiles clearly demonstrate that 
the short period spatial variations are indeed of shallow origin. 

For comparison with model anomalies in an attempt to determine the 
possible composition of the Valley, it is necessary to ignore these 
narrow anomalies. Several qualities of the present survey disqualify 
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Figure A.t. Magnetometer traverses in the Valley of Ten Thousand 
Smokes, Katmai National Monument, Alaska. Detailed 
maps of dashed areas are presented separately in 
Figures A. 9a and 5-3a. Elevation contour interval 
is 500 feet. 
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Figure 4.3a. Reduced magnetic profile data for Traverse B—B ' 1 . 

Figure 4.3b. Reduced magnetic profile data for Traverse b-b 1 . 

Figure 4.3c. Reduced magnetic profile, data for Traverse B 1 1 -b 1 1 . 

Figure 4.3d. Reduced magnetic profile data for Traverse B-B 111 . 
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Figure 4.5* Reduced magnetic profile data for Traverse E-E* 
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the data from being “smoothed 11 by most mathematical techniques. First, 
the readings are not evenly-spaced; although along the surface the sta- 
tion spacing was fairly regular, it was irregular horizontally. This 
variance violates a prerequisite for most smoothing methods. The second 
failure of the survey lies in the probable structure of the narrow ano- 
malies which are to be removed. These anomal ies are likely caused by 
fuma nolle enc rus tat i ons along compaction fractures in the py.rocl as t i cs , 
thus they can justifiably be termed neither “random" nor “periodic". 
Furthermore, the sparseness of the data was a handicap which prevented 
removal of erroneous readings and regional trends by association with 
adjacent profiles. Since the data cannot be assumed to be of any speci- 
fic mathematical form, no method of smoothing the data by fitting to a 
polynomial is justified. The only analysis technique which is justified 
in this study is the gross visual comparison of the observed data with 
profiles of made! anomalies. 

Nowhere, with the exception of the fumarole grid, Is the data 
sufficiently concentrated to justify an attempt at three-dimensional 
analysis. However, a specific type of two-di mens i onal analysis is 
suggested by the character of the observed anomalies. A preliminary 
examination of the valley profiles indicates that the causative bodies 
tend to be linear, i.e., of length s » gn i f i cantly greater than width. 

This linear trend is evident throughout the study area; within the 
individual branches of the Valley the anomalies are easily correlated 
among adjacent traverses, e.g., parallel profiles C-C' and D-D 1 in the 
southern branch of the Valley; E-E' and N l -N“' in Novarupta Basin; and 



the several profiles in Broken Mountain Valley (Figures 4 . 7 » ^*8 and ^. 9 * 
respectively). The continuance of individual anomalies from one traverse 
to the next also demonstrates that the bodies causing these anomalies 
parallel the axis of the valley. 

The mathematical formulation of a computer program designed to cal- 
culate anomalies for linear magnetic bodies for comparison with observed 
profiles, based on the methods of Sharma (1966) and Heirtzler et al . (1962) 
is presented in Appendix E. In the model study, susceptibilities have 
been assigned under the assumption of uniform magnetization in the pre- 
sent field. 

It should be noted that an erroneous interpretation could result 
if the reman ence, which can be sizable in rocks of volcanic origin, is 
significant in any of the Katmai sequences. 

In the present study, remanence is undoubtedly insignificant in the 
recent volcanic deposits, but it could be a factor in the apparent 
susceptibility of glacial drift and Naknek sediments. Since the sus- 
ceptibility of the Naknek rocks is so small, alteration of the orienta- 
tion of i ts magnet i zat i on woul d be relatively inconsequential. The 
glacial drift, which is surely primarily composed of randomly oriented 
volcanic material, could be expected to display an apparent susceptibility 
much lower than that of its parent igneous mass. Also, the volcanics 
probably lost much of their heavier magnetic minerals as they were eroded 
and deposited by glacial action. Therefore, if during the pyroclastic 
eruptions of 1912 this mass had been heated above the Curie point of 
Its components, its apparent susceptibility would increase, but not to 
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Figure 4.9a. Location of magnetometer traverses in Broken Mountain 
Valley (refer to Figure 4.1), Elevations contoured at 
200 ft. intervals. 

Figure 4.9b. Magnetic anomaly map of a part of Broken Mountain Valley 
Contour interval is 200 y. Mote the continuity of the 
anomaly pattern throughout the Valley. 



the value of the undisturbed parent volcanics, This possibility was 
dealt with previously and was taken into account in Table 4.3. 

4.2 The Character of the Flow 

An important property of the flow which has been detected by several 
geophysical techniques is its heterogeneity. Not only is the flow zoned 
vertically, as would be expected had it undergone welding, but it is 
variable laterally as well. And greater proximity to the eruptive vent(s) 
apparently even further increases the cross-sectional heterogeneity. It 
is very difficult to analyze the geophysical profiles, they are so very 
complicated. Elsewhere, the trends of adjacent geophysical profiles 
suggest that the relief of the old Valley floor is in large part responsi- 
ble for the development of the lateral and vertical variations in the flow. 

In order to propose plausible magnetic models of the flow, it was 
necessary to examine all of the geophysical data taken in the Valley. 

The magnetic susceptibility measurements (Tables 2.1, 2.2 and 4.3) and 
density determinations (Appendix C and Table 4.1) given in this paper 
serve as the basis for transforming the gravity, seismic and magnetic 
data into probable geological relationships. 

There have been relatively few density determinations made on rocks 
from the Katmat region. Appendix C lists the accumulated data to date. 
Although most of the rock types have a fairly small range of densities, 
the density of the tuff is found to vary from values comparable with 
tephra to values comparable with the banded lava of Novarupta. Also, 
the densities measured for the tuff throughout the valley show a strong 
tendency to decrease with increasing distance from Novarupta. The 
density data are summarized In Table 4.1. 



Table ^ . ] 


Summary of Density Data 
(refer to Appendix C for details) 

density (grn/cc) 


material 

low 

high 

average 

teph ra 

0.98 

1.09 

1 .03 

tuff 

1,15 

}.Sk 

va r ies* 

banded lava 

1.70 

2.22 

1.91 

rhyol 1 te g 1 ass 

2.25 

2.30 

2.28** 

andes 1 ti c lava 

2.kh 

2.55 

2 . A3 

Naknek sediments 

2A8 

2.75 

2.62 

> ,; Density of tuff 

increases as 

one a 

pproaches 


Novarupta . 

**Mean, not average. 
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One of the first tasks in analyzing the previous geophysical 
information was to re-examine the available seismic data. If travel - 
time curves were presented in the literature, these data were indi- 
vidually reduced to comp ress i onal velocities and thicknesses for each 
seismic layer {Dobrin, I960, pp . 70"S3) ; otherwise, the smoothed data 
as presented by the original investigator was used. A composite of 
this seismic information Is presented in Appendix D and summarized in 
Tab le 4 . 2 . 

It should be understood that the reduction of seismic refraction 
data assumes that each successively deeper layer possesses a greater 
compress i ona 1 velocity than the overlying layer. If this condition is 
not met in nature, then the depths and thickness for any layers below 
such a low velocity layer will be in error since the thickness of the 
low velocity layer. has not been taken into account. 

Since it is suspected that at least a portion of the flow is welded 
the general nature of welding within an ash flow should be understood. 
According to Smith (i960, p. 831 ), a single ash flow (one cooling unit) 
may display three basic zones. These three types are: no welding, 

partial welding, and dense welding. Emplacement temperature and flow 
thickness are critical factors in determining the character of the zon- 
ing. Usually, when dense welding occurs in a unit, it is enveloped in 
a zone of partial welding, which in turn is surrounded by unwelded 
material. The lower zones of no welding and partial welding are thinner 
than their upper counterparts. Special or extreme conditions of tempera 
ture and pressure can cause some of the zones to be local ly absent. 

Therefore, since the degree of welding is directly proportional 



Table 4.2 


Summary of Seismic Parameters 


Seismic 
Laye r 
Numbe r 

P robab 1 e 
Compos i t i on 

H i ghes t 
Compress ional 
Ve 1 oci ty 
(km/sec) 

Lowes t 

Compress i onal 
Ve 1 oci ty 
(km/sec) 

Number 

of 

Ve 1 oci ty 
Read i ngs 

Ave rage 
Ve 1 oci ty 
( km/sec) 

Number 

of 

Th t ckness 
Read 1 ngs 

Ave rage 
Th i ckness 
(m) 

1 

vjea the red layer 

0.29 

0.1 

13 

0.16 

13 

1 .2 

2 

ai r-f a 1 1 

0.59 

0.25 

30 

0.38 

30 

5.7 

3 

. zone of no we 1 d f ng 

0.74 

0.43 

28 

0.61 

28 

26.6 

4 

zone of partial welding 

1.2 

0.8 

18 ■ 

0.98 

13 

41 .4 

5 

zone of dense welding 

2.44 

1.5 ' 

30 

1 .94 

13 

48.4 

6 

g 1 ac I a 1 drift 

3.36 

2.5 

17 

2.88 

-- 

— 

7a 

Naknek bedrock 

3.8 

3.68 

6 

3 . 73 

-- 

— 


4.4 


4.5 


7b 


igneous bedrock 


4.6 


2 
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to the density within the flow, and seismic velocity is directly pro- 
portional to the density, it is probable that the zones of welding serve 
as seismic layers. Furthermore, the depth and thickness of seismic 
layers beneath a welded zone can be expected to be In error, for the 
thin lower border of less welded or unwelded material (less dense and 
of lower velocity) would not be detected by seismic refraction. Should 
the flow consist of more than one cooling unit (be underlain by a pre- 
vious welded flow(s) of greater age), then the hierarchy of welded and 
unwelded zones could become complicated. The multiple low velocity 
layers which could be entrained in such a section would greatly decrease 
the accuracy of a thickness determined by seismic refraction. 

There are at least eight different seismic velocity layers repre- 
sented in the Valley. In a few cases an uppermost thin (1.2 meters) 
layer of very low velocity (average 0,16 km/sec) material was observed. 
Surely this covering Is composed of uncompacted secondary deposits of 
a (wind and water) weathered nature. Although probably present through- 
out the Valley, this layer went undetected along most profiles due to 
its thinness, the depth of the shot and the spacing of the first few 
geophones . 

in most cases the top seismic layer which was detected has an 
average velocity of O .38 km/sec. In each case where 0.16 km/sec material 
was detected, it was underlain by a layer of O .38 km/sec material. This 
second layer is probably air-fall. Again it is a relatively thin section 
(averaging 5.7 meters thick) in most cases. Along a. few profiles it 
was not detected (the cause is probably the same as that cited for the 
common failure to detect the first layer). 
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Usually the third layer had an average velocity of 0.61 km/sec. 

(However, in a few cases the. next layer was of O .98 km/sec or 1.94 km/ 
sec material; and, on the bench on the western side of Baked Mountain, 

30 meters of 0.30 km/sec material is underlain by a section of 3-2 km/ 
sec material.) The 0.61 km/sec layer is relatively thick and probably 
constitutes an unwelded zone within the pyroclastic flow. Generally, 
this zone is underlain by O .98 km/sec or 1.94 km/sec material, or both; 
the exceptions are a few places in the southern branch where velocities 
of 2.5 to 3.0 km/sec are observed for the next layer, and along one pro- 
file in Novarupta Basin where 73 meters of 0.65 km/sec material overlies 
a section of 4.4 km/sec material. 

When detected, the O .98 km/sec layer is usually overlying 1.94 km/ 
sec material. In the southern branch there are a few instances where 
the velocity of the next layer ranged from 2.5 to 3-1 km/sec, and near 
the terminus of the flow the bottom layer has a velocity of 3-8 km/sec. 

The conspicuous columnar jointing within the flow near the terminus is 
indicative of partial welding. In this area, the bulk of the flow possesses 
a seismic velocity near 1.0 km/sec. It is thus logical to assume that 
the layer with an average velocity of 0.98 km/sec represents a zone of 
partial welding within the pyroclastic flow. The average thickness ob- 
served for this layer is 41.4 meters. 

The presence of dense welding is probably indicated by velocities 
in the range of 1.9*1 km/sec. Often this seismic layer was the deepest 
horizon detected by the refraction surveys; however, in the few instances 
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when deeper sections were observed, it is possible to determine that this 
7 . 0 ne can range from 33 meters to 73 meters thick. 

The uniqueness of the upper five layers ( 1 — 0 . 1 6 km/sec; 2--0.38 km/ 
sec; 3 — 0 .61 km/sec; 4--0.98 km/sec; and 5"- 1.94 km/sec) is demonstrated 
in Section 5b of Sbar and Matumoto's Profile 2 in the south branch of 
the Valley; elsewhere, layers 1 through 4 and 2 through 5 are encountered 
along individual sections (refer to Appendix D) . The failure to detect 
the uppermost layers along sane profiles is no doubt due to their thin- 
ness relative to the surveying parameters. The absence of 1.94 km/sec 
material indicates that conditions within that section of the flow were 
not conducive to the formation of a densely welded zone. The absence of 
both the fourth and fifth layers shows that no welding occurred in that 
portion of the flaw. The absence of layers three and four above a layer 
of 1.94 km/sec material indicates that this part of the flow was the 
site of extreme temperature or pressure conditions. 

The sixth seismic layer encompasses quite a range of velocities and 
is probably the glacial drift mantling the Naknek bedrock. Across the 
southern branch of the Valley, the velocity of the deepest horizon de- 
tected by Gedney et al . (1970) ranged from 2.62 to 3-0 km/sec. Similar 
velocities were encountered elsewhere in the southern branch, in the 
middle Valley, and in the lower halves of Broken Mountain Valley and 
Novarupta Basin, Since no deeper horizons were ever observed beneath 
this material, its thickness can not be determined. When present, this 
glacial drift constitutes layer six. 

Curtis (personal communication, 1971) claims that seismic velocities 



of 2,8 km/sec are not out-of-line for densely welded tuffs. If this con- 
jecture is expanded to include all of layer 6 (average velocity 2.88 km/ 
sec), it is evident that the pyroclastic flow might be substantially 
thicker than previously assumed. If layer 6 is actually the densely 
welded zone in an earlier cooling unit, it would probably be enveloped 
in lower density, less-welded material (which would constitute low velo- 
city layers and would not be detected by seismic refraction surveying) 
and the total thickness of the deposits may be great indeed. 

The bedrock in the region of the Valley is the Naknek sedimentary 
strata. Near the end of the flow, it is exposed along the river gorges 
beneath less than 50 meters of py rod as t i cs . In this area, Matumoto 
and Ward ( 196 '/, p. 121 ) detected a base horizon of 3-8 km/sec at a 
depth of 46 meters. Similar seismic velocities were observed below zones 
of dense welding in the upper half of Broken Mountain Valley and near 
the edge of the bench along the western flank of Baked Mountain in the 
southern branch of the Valley. Since it is likely that' the glacial 
accumulations were slight in these two places, it seems probable that 
seismic velocities near 3.7 kni/sec are indicative of the Naknek sedi- 
mentary horizon. The density of the Naknek sediments ( 2.6 gm/cc) is in 
accord with such a seismic velocity (Grant and West, 1965 , p. 200 ). 

Only twice were higher velocities observed. Near Novarupta, Sbar and 
Matumoto (1370 recorded velocities of 4.4 and 4.6 km/sec. Nearby 
igneous rocks have densities near 2.5 gm/cc, which could possess seismic 
- compress l ona 1 velocities In this range. It is proposed that in these 
two instances, Sbar and Matumoto penetrated to an igneous stratum of 
bedrock. Thus, the bottommost horizon is the Naknek bedrock which composes 
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layer 7a; or, near the volcanoes, the Igneous horizon acts as layer 7b. 

A few general statements can be made about the thickness of the 
pyroclastic deposits. The possible inclusion of low velocity layers 
means that the seismic information can serve only as minimum estimates 
of the flow thickness. This restriction should be kept In mind through- 
out the following discussion. 

The minimum thickness of the flow increases from approximately 50 
meters at the terminus to about 70 meters in the middle and southern 
branches, to 100 meters in Broken Mountain Valley and Novarupta Basin, 
and perhaps even more in the southeast branch. Generally the air-fall 
pyroclastics constitute less than 10 meters of the pyroclastic flow. 

On the bench on the west of Baked Mountain the air-fall layer is some 
25 meters thick. These unusually thick deposits are probably in part, 
due to slumping of the ashy material from the adjacent steep slopes. 

Unfortunately, there is no case where a profile containing all 
three zones of tuff extends to horizons 6 or 7* Therefore, we do not 
have thicknesses for all 3 zones within any one section. The available 
profiles, however, permit some conclusions concerning the zoning within 
the flow. Where underlain by layer k f layer 3 averages about 15 meters 
thick; layer A ranges from 2k to 62 meters thick in this case; the 
absence of layer 3 makes no discernible difference in the recorded 
thicknesses for layer k\ where layer 5 is in immediate contact with 
layer 3, the latter layer averages nearly 25 meters thick; there is 
no marked difference in the thicknesses detected fo'r layer 5 whether or 
not it is overlain by layers 3 or A (it varies from 33 to 73 meters in 
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both cases); and where not underlain by either layers 4 or 5, layer 3 was 
observed to reach thicknesses of up to 73 meters. 

Layer 4 was not observed in Broken Mountain Valley, although the 
flow is thick, and possesses a zone of dense welding. Assuming that an 
ancestral Novarupta (refer to Chapter V) stood higher than Greasy Pass, 
one could assume that Broken Mountain Valley would receive a large por- 
tion of the flow, and at a higher temperature than the more distant 
branches of the Valley. The high emplacement temperature could be re- 
sponsible for the absence of layer 4. The flow is thinner in the southern 
branch and farther north. The seismic evidence for the southeast branch 
is inconclusive since layer 5 was the deepest horizon detected, but here, 
too, layer 4 is absent. Curtis (1968) showed that most of the teph ra 
eruption was directed northeast away from Novarupta; this would be into 
Broken Mountain Valley and the southeast branch. If the expulsion of 
the tuff was similarly directed, the accumulations would be thicker and 
hotter in this branch as well as in Broken Mountain Valley. This supposi- 
tion agrees well with the meager seismic evidence. According to this 
hypothesis, the portion of the tuff which flowed down Novarupta Basin 
and thence into the southern branch was somewhat cooler, and cooled fur- 
ther with travel. The presence of layer 4 in these areas and farther 

J 

north supports this idea. 

Kienle's profile D-D’ shows no layer 3* The high temperature near 
Novarupta apparently indurated most of the flow. Closest to the dome 
the flow consists of 16 meters of partial welding and 73 meters of dense 
welding; whereas about haifway out of the basin, the ratio has reversed 



108 


to 43 (Deters of layer 4 over 33 meters of layer 5. In the southern branch, 
layer 3 is often present as well as layers 4 and 5. And, near the termi- 
nus of the Valley where the flow was the thinnest and coolest, there Is 
no evidence of dense welding. 

The direction of the eruption and distances from the eruptive vent 
affected the thickness and temperature of the flow, and thus the character 
and extent of the welding. The thicker and/or hotter the flow, the more 
extreme and extensive the welding. 

Comparison of magnetic and gravimetric data along adjacent traverses 
(refer to Figures 4.10, 4.11 and 4,12) clearly demonstrates correspondence 
--usually the two data forms appear as mirror-images of one another with 
a gravimetric low corresponding to a magnetic high and vice versa. One 
explanation for this, feature would be to assume that the py rocl as tics 
deposits and the glacial drift above the bedrock (Naknek sediments are 
of high density and low susceptibility) are of uniform thickness along 
any one profile, with lateral groupings of dense material of low suscep- 
tibility. Such a situation is contrary to density and susceptibility 
measurements to date, i .e . , the denser pyroclastic material has a tendency 
to be of higher s uscept ib i 1 i ty . 

Generally, the pyroclastics In the Valley tend to be of relatively 
low density and high susceptibility in compa r i son wi th their sur roundings . 
Increasing the thickness of the underlying deposits therefore increases 
the magnetic anomaly and decreases the gravimetric anomaly. in fact, 
variable thickness is probably the primary cause of the anomalous gravi- 
metric and magnetic trends across the various branches of the Valley. 
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Figure ^.10, Correlation of adjacent geophysical profiles in the 
southern branch of the Valley. Refer to Figures 
1 . 3 and k . 1 . 
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Figure 4.11a. Correlation of adjacent geophysical profiles in the lower 
valley. Refer to Figures 1.3. and 4.1. 

Figure 4.11b. Correlation of adjacent geophysical profiles in Broken 
Mountain and Knife Creek Valleys. Refer to Figures 1.3 
and 4.1. 
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Figure ^4 . 12a. Correlation of adjacent geophysical profiles in 
Novarupta Basin. 

Figure ^.12b. Correlation of adjacent geophysical profiles in 
Broken Mountain Valley. 
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Since the top of the flow is quite level, the geophysical anomalies Im- 
ply that base has considerably relief. In all cases, the greatest ano- 
malies occur in the vicinity of the present streams. This evidence sup- 
ports the suggestion of previous investigators (Kienle, 1970, p. 6657; 
Curtis, 1968, p. 186) that the present streams overlie the pre-existing 
stream channels. it would seem then, that the gravimetric and magnetic 
anomalies can serve to delineate the configuration of the pre -1912 valley 
floor. The undulating nature of the anomalies in the southern branch 
may indicate faulting previous to the final emplacement of the flow. 

Cursory examination of the seismic profiles corresponding to the 
gravimetric and magnetic profi les shows greater tuff thicknesses and 
higher degree of welding in conj unct i on wi th the magnetic highs. This 
complementary evidence Is further proof that in profile the base of the 
flow Is not level and that the buried deeply Incised stream channels 
underlie the present streams. 

Kienle (I 969 ) determined two possible models capable of producing 
the observed gravimetric anomalies (refer to Figure ^.13). The density 
contrast between the flow and bedrock for his model P now appears to 
have been too great ( 1.59 gm/cc), as a density of about 1.6 gm/cc is nor- 
mal for the tuff, and Naknek sediments have an average density of 2.6 
gm/cc. The flow is probably separated from the Naknek bedrock by glacial 
drift, which would, however, be of lower density than the sediments. 
Kienle's model F with a density contrast of 0.75 gm/cc probably more nearly 
represents the actual situation. The tuff thicknesses suggested by model 
F correspond well with Curtis 1 (1968) estimates based on geomorphol og i c 
extrapolation (refer to Table 1 . 3 ) - 
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Figure ^.13. Summary of Kienle's 09&9) structure sections 
derived from gravity data. 


in most branches of the Valley, the seismically indicated flow thick- 
nesses are less than Curtis and Kienle calculated. Perhaps the seismic 
sections are too thin for they have not allowed for undetectable low ve- 
locity layers of substantial thicknesses; perhaps seismic layer 6 is also 
part of the pyroclastic flow. On the other hand, the seismic data show 
thicker pyroclastic sections in Broken Mountain Valley than Curtis and 
Kienle predict. Curtis had no base control for this branch of the Valley 
and apparently he did not appreciate the depth of the base in this area. 
Although Kienle’s data were poor in Novarupta Basin and Broken Mountain 
Valley, his preliminary thickness estimates for these branches seem ano- 
malously low. Higher emplacement temperatures for the flow near Novarupta 
(e.g., these two branches) would result in the flow being more densely 
welded and therefore more dense. The tuff densities were found to increase 
with proximity to the dome, with the highest (1.96 gm/sec; Kienle, personal 
communication, 1971) from Novarupta Basin (refer to Appendix C) . Thus, 
in these two branches the contrast between tuff and bedrock is lower 
than Kienle’s model F; i.e,, the flow is thicker than suggested by model 
F. 

I 4.3 Magnetic Models of the Flow 

Some of the parameters remained constant In all of the computer 
models. A total field of 53,800 gammas, wi t:h declination of 2l.5°E end 
inclination of 71*2°, was used. All the magnetic bodies were assumed 
infinite along an average, valley axis of 30 °W. 

Four widths were used in modeling the anomalies: 250, 500, 750 and 

10 00 meters. Most of the observed anomalies have a lateral extent con]- 
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parable to one of these values. The depths and thicknesses of the model- 
ing bodies are based on the seismic, layers. The probe of the magneto- 
meter was approximately 2 meters above the surface of the flow. The 
air-fall deposits generally account for an upper section of about 3 meters, 
below which the body of the tuff amounts to 25 meters or more. A 10 
meter wedge below 50 meters could represent a thermal ly-al tered zone of 
glacial drift and/or Naknek bedrock. Additional wedges of 50 and 100 
meters (representing greater flow thickness, internal zoning and/or a 
layer of glacial drift) permit calculation for total thicknesses of up 
to 300 meters (or, of up to 900 meters for 750 meter wide bodies). 

As is clear from the derivation in Appendix E, magnetic suscepti- 
bility is a constant, altering only the magnitude of an anomaly, not 
its form, if all other factors are static, it is possible to calculate 
the amplitude of the anomaly over an identical body of different sus- 
ceptibility by merely multiplying by the ratio of the newly-assigned sus- 
ceptibility to the previously-assigned susceptibility (refer to Table 4.3). 
Of course, the converse scheme permits the deduction of the suscepti- 
bility required to produce an anomaly of specific magnitude if all the 

other parameters of the body are constant. 

-6 

Since a susceptibility of 1250 x 10 emu/cc aptly represents air- 
fal 1 and tuff, and may well approach the apparent susceptibility of the 
glacial drift, It was the best choice for the initial model studies. 
Appendix E.3 lists the total magnetic anomalies at increments of 20 
meters across various wedges of 1250 x 10 ^ emu/cc material. The cen- 
tral anomalies of these wedges is summarized in Table 4.4. 



Table 4.3 


Rock Types in the Valley of Ten Thousand Smokes, 
Their Modeling Susceptibilities, and Conversion 
Factors (for Calculating the Magnitude of the 
Magnetic Anomaly Due to Susceptibilities 
Other than 1250 ^ 10“° emu/cc 
f rom Append i * E . 3) 


Mater! a 1 

Naknek sediments 
fumarolic altered ash 
rhyol i te 
ai r-fall 
tuff 

banded volcanics 
glacial drift 
andes i te 
altered Naknek 
al tered drift 

*Sus cept i b i 1 i ty i s i n 


Suscept i b I 1 i ty* 

50 

50 

250 

1250 

1250 

1250 

( 1250 ?) 

2500 

5000 

(5000?) 

its of 1 0 -6 emu/ cc . 


Conversion Factor 

0 . 04 
0.04 
0.20 
1 .00 
1 .00 
1 .00 
1 .00 
2.00 
4.00 
4.00 
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Table 4 . 4 

Summary of Central Total Anomalies (in Gammas) for 
Computerized Magnetic Bodies of 
Susceptibility - 1250 x 10” fc emu/cc 
at a Depth of 2 Meters 
(refer to Appendix E.3) 


Thickness (m) 
Width (m) f 

3 

28 

53 

63 

103 

203 

298 

250 

253 

328 

364 

378 

423 

488 

519 

500 

2 50 

23k 

316 

324 

355 

414 

451 

750 

250 

297 

312 

318 

340 

389 

426 

1000 

249 

292 

303 

308 

325 

365 

398 
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Unfortunately, It is not possible to derive the unique thickness, 
susceptibility and depth responsible for an anomaly of specified width 
and amplitude. Therefore, these parameters must be studied separately 
within the reasonable limits of the other two parameters. 

Initially, the reduced magnetic anomalies were examined by divid- 
ing the profiles into as few laterally separated anomalies as seem pro- 
bable, and taking their magnitude relative to the regional magnetic 
field only (refer to Figure 4.14). Table 4.5 lists these values. As 
is clear from Appendix E.3, the magnetic anomalies at the high magnetic 
latitude of Katmai and along the specified traverse orientations die 
off rapidly beyond the wedges. It is therefore possible to neglect the 
contribution of neighboring wedges without markedly influencing the 
quality of the results. First, it was assumed that each anomaly was due 
to a wedge of 1250 x 10 ^ emu/cc material at a depth of 2 meters and of 
the modeling width nearest the width of the observed anomaly. From this, 
the nearest modeling thickness capable of producing the observed anomaly 
was determined. In several instances the observed anomaly was of in- 
sufficient magnitude to be the resultant of even a 3 meter thickness of 
-6 

1250 x 10 emu/cc material at a depth of 2 maters. In these cases, the 
top of the causitive body is assumed to be at a depth of 5 meters. In- 
creasing the depth of the wedge requires a sizeably greater thickness 
to produce the same magnetic anomaly. 

It is clear from Table 4.4 that the susceptibility of the uppermost 
few meters is critical; e.g., the anomaly caused by a 3 meter thick by 
250 meter wide wedge at a depth of 2 meters is app roxi mate ly equal to 
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Figure b,]k. Smoothed residual magnetic anomalies and their correspond i ng 
simple and local anomalies along the valley traverses (see 
text) . 
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that from a 295 meter thick by 250 meter wide wedge at 5 meters, and a 
3 meter thick by 750 meter wide wedge at 2 meters yields an anomaly 
similar to that due a 595 meter thick by 750 mater wide wedge at 5 meters 
dep th . 

If fumarolic leaching was especially prevalent in the upper por- 
tion of the flow, one would expect the material to be of very low sus- 
ceptibility. In such a case, the top few meters would contribute a 
negligible amount to the magnitude of the anomaly; this situation could 
be responsible for the unusually low anomalies. It is unreasonable 
however to expect that such leaching should be restricted to the precise 

region of the lower anomalies. If the leaching is widespread, the im- 

-6 

plication is that ridiculously thick sections of 1250 x 10 emu/cc 
material would be responsible for all the anomalies. The initial 
determination of the thickness of the flow (given the approximate width, 
depth and probable susceptibility) yields some erratic results across 
various branches of the Valley. These models do suggest, however, that 
the flow may be some 300 meters thick. 

Another method of analysis is to assume a specific thickness, depth, 
and width for the causitive body. Two estimates of the thickness of the 
flow in each branch were chosen from the conclusions of previous investi- 
gators (refer to Table *4.5) the depth to the top of the flaw was taken 
as 2 meters and the modeling width nearest the width of the observed 
anomaly was used. Given the above factors and the magnitude of the 
observed anomaly, It is possible to calculate the magnetic susceptibility 
required of the body to produce the anomaly. The calculated susceptibilities 



Table 4.5 

Simple Magnetic Anomalies and Three Computer Models 




OBSERVED PARAMETERS 
Anomaly 

Meores t 


An pm. 

Anomaly Width 

Magn i tude Mode ling 

T raverse 

£ 

(m) 

(y> 

Width (m) 

A- A' 

1 

600 

400 

750 

A- A 1 

2 

1800 

200 

10CO 

B-3 1 

1 . 

650 

400 

750 

B-B*. 

2 

. jCD 

300 

250 

8-3’ 

3 

1350 

400 

1000 

3- 3 1 

3a 

350 

550 

250 

3-3 ’ 1 ' 

1 

850 

150 

750 

0-B" ' 

1 

goo 

300 

1C00 

C- C ' 

1 

1100 

250 

lo 0.0 

C-C 1 

2 

250 

500 

250 

c-c 

3 

500 

250 

500 

C-C’ 

A 

500 

2pO 

500 

0-D’ 

\ 

£00 

300 

500 

D-D' 

1 

550 

425 

500 

D-D' 

3 

250 

200 

250 

D-D’ 

A 

350 

250 

250 

£-£' 

* 

“ ~ " 


all 

negatives 

Model 1 : 

g i ven 

suscept i b 1 1 1 ty 

- 1250 X 

10*^ cmu/cc; 


width 

» nearest modeling width 

, and the mag 


n 1 1 u de 

of the observed anomaly 

, and specify' 


ing the depth to the 

top (D) , 1 

whot is the 


required thickness, 

1 . e . , what 

is the 


( D+ Du) 

nearest model 

? 

ing depth 

to the base 


MODEL I MODEL I I MODEL I 1 I 


D 

(D+DD) 

(D+DD) 

K , 

(D+DD) 

, X _ 6 

(m) 

<m) 

(m) 

(emu/cc x 10 °) 

(m) 

(emu/cc x 10 

2 

250 

55 

1620 

155 

1 366 

5 

>300 

55 

826 

155 

725 

2 

250 

55 

1602 

105 

1470 

2 

30 

155 

■ 84! 

250 

746 

2 

300 

155 

3 449 

250 

1312 

2 

>300 

155 

1822 

250 

1616 

5 

250 

55 ‘ 

600 

105 

552 

2 

55 

55 

1233 

105 

1154 

5 

>300 

55 

1032 

155 

506 

2 

2 50 

55 

1717 

155 

1401 

5 

>300 

55 

S2C 

155 

812 

5 

>300 

55 

390 

155 

812 

2 

30 

55 

1 187 

155 

■ 974 

2 

250 

55 

1681 

155 

1 380 

5 

205 

55 

686 

155 

56 1 

5 

300 

55 

858 

155 

7!0 





— 

— 

— 


Models M and III: given width = nearest modeling 

width, the magnitude of the ob- 
served anomaly, and the depth 
to the top of the body = 2 m. , 
and specifying the depth to the 
base (D+DD), what is the required 
suscept I b i 1 i ty (K) 7 
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are all within the range of the measured susceptibilities, but it is 
difficult to explain such lateral alteration in the magnetic properties 
of the flow. Evidently, there may be down -valley and cross -valley 
compositional changes in the flow. 

A casual inspection of the profiles of the residual field shows 
that the regional trends are significant along several profiles. Kienle 
(I969) was able to remove the regional gravity trends because he had 
bedrock stations at the ends of his traverses. Such a technique is 
not possible in the present magnetometer study. Neither are the mag- 
netometer traverses sufficiently dense to permit delineation and removal 
of the regional field by mathematical techniques. Instead, the regional 
trends have to be arbitrarily, visually removed in order to obtain the 
local magnetic anomalies. Figure A .14 exhibits the regional trends and 
the resulting local anomalies. These anomaly parameters are listed in 
Tab 1 e k .6 , 

These anomalies were analyzed as above; first, to determine the 

-6 

necessary thickness of 1250 x 10 emu/cc material, and next to determine 
the necessary susceptibility in two cases. The results we re no more 
satisfying than those from the examination of the residual anomalies in 
Table k. 5 . The few instances in which seismic data is available near 
magnetic profiles permits a check on the plausibility of the magnetic 
models. Table k.l presents these correlated data. Comparisons of the 
computerized magnetic models with the seismic information reinforces 
many of the earlier indications of Tab les h . 5 and 4 . 6 . 

If the seismic thicknesses are assumed to be reasonable estimates 



Table 4.6 . 


Multiple Local Anomalies end Three Computer Models 


OBSERVED PARAMETERS 1 





An oma 1 y 

Neares t 



(D+CD) 

(m) 


Ar.om. 

Traverse §• 

Anomaly Width 

Cm) 

Megn i tude 

(t) 

Model log 
Width (m) 

D 

(m) 

(d+od) 

(m) 

(emu/cc x 10 

A-ft* 

1 

620 

235 

SCO 

2 

30 

55 

1 125 

A- A' 

2 

730 

230 

750 

5 

500 

55 

222 

A-A' 
A- A 1 

3 

4 

3 CO 
600 

170 

285 

250 

500 

5 

1 

105 

30 

55 

55 

1125 


5 

200 

55 

250 

5 

30 

55 

189 

n_A 1 

1 

4C0 

455 

500 

2 

300 

55 

1 SCO 


2 

300 

585 

250 

2 

>300 ■ 

53 

2050 

3-£ l 

3 

250 

515 

250 

2 

300 

155 

1 

a- £ 1 

4 

400 

455 

500 

2 

300 

155 

1 473 

B-B 1 
3 1 

5 

6 

350 

550 

740 

340 

250 

500 

2 

2 

>300 

105 

155 

155 

2038 

1104 

3-8' 

7 

350 

170 

250 

5 

105 

155 

477 

5- 3 ’ 1 

1 

eso 

340 

750 

2 

105 

55 

1362 

b-D’ ' 

’ 1 

?2 0 

455 

1000 

2 

400 

55 

1879 

C-C 

c-c 

r- r * 

! 

2 

2a 

750 

750 

200 

265 

230 

585 

750 

750 

250 

2 

5 

2 

30 

500 

>300 

55 

55 

55 

1142 

322 

2010 

r> r r 

3 

450 

270 

500 

2 

30 

55 

1069 

r r * 

4 

530 

340 

500 

2 

105 

55 

1345 

c-c 

5 

4 00 

170 

500 

5 

205 

55 

672 

0-0 1 
D-D' 
n- n 1 

1 

2 

2a 

650 

650 

350 

400 

170 

570 

750 

750 

250 

' 2 
5 
2 

250 

300 

>300 

55 

55 

55 

1602 

681 

1258 

p-p< 

3 

300 

200 

250 

5 

250 

55 

686 

D-D' 

D-D' 

4 

5 

750 

550 

170 

115 

750 

500 

5 

5 

300 

105 

55 

55 

68 1 
488 


) 

200 

115 

250 

5 

55 

105 

340 


2 

250 

170 

250 

5 

105 

105 

503 

r r } 

3 

1 50 

585 

250 

2 

>300 

105 

1 729 

f - r 1 

4 

400 

340 

Soo 

2 

105 

105 

1138 

E-E' 

E-E’ 

5 

6 

250 

280 

115 

230 

250 

250 

5 

5 

55 

250 

105 

105 

340 

679 

{The 

models are 

based on the 

same requirements as gl 

ven In 

Table 

4.5.) 



MODEL 1M 


( 0 +DD) 

K -6 

(m) 

(emu/ce x 10 

155 

226 

155 

736 

• 155 

477 

155 

326 

155 

i 54 

105 

1 605 

105 

1723 

250 

1279 

250 

1316 

250 

1833 

250 

532 

250 

422 

105 

1250 

10 S 

1750 

155 

375 

155 

786 

155 

1640 

155 

878 

155 

1104 

155 

552 

155 

1 366 

155 

581 

155 

1562 

155 

548 

155 

581 

155 

373 

155 

322 

155 

477 

155 

1640 

155 

. 1104 

155 

322 

155 

645 



Table 4,7 


Correlation of Four Computerized Magnetic Models with Observed Local 
Magnetic Anomalies and Nearby Seismic Refraction Profiles 





Model i ng 




Th i ckness 

of 



Mode 1 1 ng 

Model 1 

Mode! 1 i 

Model 11! 

Mode 1 1 V 




Width 

Magn i tude 



Se 

i sm i c 




Thickness Nearest 


Magn i tude 






Meares t 

of Local 



Layers (m)* 



Total Seismic 

(D+00) 

of 

K 

(D+00) 



Anom. 

Observed 

Anomaly 








Th i ckness 


Anomaly 

c 


Locat i on 

T reverse 

U 

Width (m) 

(y) 

1 

2 

3 

4 

5 

6 

7a 

(m) 

(m) 

(y) 

(emu/cc x 10 ) 

(m) 

Lower Valley 

A* A' 

2 

750 

230 

C ) 

2 

7 

37 

7 

— 

7 

— 

<5 

— 

--- 

... 

Southeast Branch 

3-5 ' 

3 

250 

515 

( ) 

1 

18 

--- 

n 

{ ) 

{ ) 

— 

300 

... 

... 

105 

Sou theas t 3 ranch 

B-S * ' 

5 

250 

740 

( ) 


22 

— - 

7 

( ) 

( ) 

— 

>300 

... 

— 

300 

Broken Mtn. Valley 

B-3 ’ 

1 

500 

455 

1 

2 

25 

— 

? 

( ) 

( ) 

(105) 

300 

355 

1605 

130 

5 rcker, K* n. Va 1 1 cy 

6-3 ’ 

2 

250 

535 

1 

2 

25 

— 

7 

( ) 

( ) 

(105) 

>300 

423 

i 72® 

155 

Broken M.tn. Valley 

S-B’ ' 

1 

730 

340 

3 

5 

25 

— 

53 

7 

( ) 

105 

105 

340 

12=0 

55 

Broken Mtn. Valley 

B-B" ’ 

1 

1000 

455 

2 

7 

22 

— 

75 

— 

? 

105 

400 

325 

1 750 

205 

Scut hern Branch 

C-C' 

I 

• 750 

235 

( ) 

1 

47 

— 

... 

? 

( ) 

55 

30 

312 

i 142 

... 

Southern Branch 

C-C’ 

2 

750 

230 

( ) 

4 

43 

— 

— 

7 

( ) 

55 

<5 

312 

222 

— 

Southern Branch 

C-C 

2a 

250 

585 

( 5 

2 

14 

25 

? 

(’) 

( ) 

— 

>300 

— 

— 

155 

Southern Branch 

C-C‘ 

3 

500 

270 

( ) 

5 

31 

— 

33 

7 

( ) 

65 

30 

324 

1041 

— 

Southern Branch 

C-C 

4 

500 

340 

( ) 

2 

— 

105 

— 

? 

( ) 

105 

105 

355 

1 12$ 

— 

Southern Branch 

C-C ' 

5 

500 

170 

( ) 

25 

— 

— 

— 

? 

( ) 

30 

<5 

297 

723 

— 


( ) 
Model 

Model 


Model 


net detec ted 

detected, but no thickness information 
su s peered 

£ 

I: given depth to top = 2 u, susceptibility ■ 1250 x 10 emu/cc, width Q nearest modeling width, and the magnitude of the observed anomaly, what 

is the required nearest modeling depth -d the base, (D+DD)? 

li: given depth to top - 2 m, susceptibility a 1250 x 10 "& emu/cc, width = nearest modeling width, and (D+DD) = modeling thickness nearest the total 
seismic thickness, what is the resultant magnitude of the anomaly? 


Ml: 


given death to top = 2 m, (D+DD) = modeling thickness nearest the total seismic thickness, width, a nearest modeling width, and the magnitude 
of the observed anomaly, what is the required susceptibility it? 


given a body of the nearest modeling width at a depth of 1 m extending to 30 m and of susceptibility ™ 1250 x 10"^ emu/cc, directly 
by a body cf the same width but of susceptibility - 2500 x 10"° emu/cc, what is the nearest modeling depth to the base, (D+00), of 
second body in order that the anomaly be of the observed magnitude? 


unde r 1 a i 
the 


n 


Model (V: 
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of the true thickness of the flow along the traverses , then one must con- 
clude that the susceptibility of the flow Is variable. Clearly, the al- 
terations in the flow thickness are insufficient to be the sole cause of 
the local magnetic changes, 

Concentrat i ons of material of substantially lower susceptibility 
-6 

than 1250 x 10 emu/cc are to be found in the Lower Valley, near the 
borders of the flow in all branches of the Valley, and in Novarupta Basin. 

tt should be understood that the modeling susceptibility (1250 x 10 ^ 
emu/cc) for the flow was based on air-fall and tuff samples from the 
upper regions of the flow only. No samples we re obtained from the Lower 
Valley where the flow has changed in appearance and is more homogeneous 
and more rhyolitic. 

It has been .found that the susceptibility of rhyolitic material is 
relatively low. Obviously, were an appreciable amount of the flow pri- 
marily rhyolitic, the resulting magnetic anomalies would be far smaller 
than those calculated for a flow comprised entirely of the modeling 
susceptibility. Thus, the more rhyolitic composition of the flow near 
its terminus is responsible for its lower susceptibility and resulting 
lower magnetic anoma 1 i es .. Extensive pockets of rhyolitic material near 
their source, Novarupta, could account for the indicated low susceptibi- 
lities in Novarupta Basin. 

A different explanation for the longitudinal pockets of low suscep- 
tibility material along the borders of the Valley is probable. Much of 
the ash along the borders is thought to be weathered debris which has 
slumped from the adjacent steep slopes. Perhaps the prolonged weather- 
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i rig of these portions of the flow which were originally deposited on the 
steep slopes has leached the ash sufficiently to decrease its suscepti- 
bility almost by half, as is suggested by the magnetic modeling studies. 

As was previously stated, it is doubtful that fumarolic leaching 
could be extensive enough to account for the observed anomalies. Fuma- 
rolic activity was never observed to be continuous over areas as large 

as those outlined by the magnetic anomalies. On the other hand, In many 

-6 

cases a sei smi cal ly-determi ned thickness of 1250 x 10 emu/cc material 
is insufficient to produce the observed local magnetic anomaly; greater 
thickness and/or higher susceptibility is required. 

Frequently the highest magnetic anomalies were encountered near 
streams. If greater thickness is assumed responsible, then one would 
conclude that the previous streams must underlie the present gorge. 

Spurr (1900) reported that the previous streams were incised about 30 
meters into the old valley floor of glacial drift. But even an extra 
50 meters of the modeling material below a minimal flow thickness of 
50 meters cannot produce the observed anomaly. A wedge 250 meters wi de 
at a depth of 55 meters could increase the anomaly by 95 gammas; far 
less than the difference of about 200 gammas frequently found near the 
streams. The models in Tables A.5, A. 6 and A. 7 sho*w that increased 
thicknesses of over 100 meters are required to account for the "stream" 
anomalies, if the susceptibility remains at 1250 x 10 ^ emu/cc through- 
out the flow. Therefore, one must conclude that at least a portion of 
the flow must be of greater susceptibility than the modeling suscepti- 
bility. 

There has been some speculation that magnetite concentrations along 
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the present stream beds are responsible for the high anomalies encountered 
nearby. Certainly even a thin wedge, with a susceptibility as high as 
has been measured for magnetite, could easily account for the magnitude 

of the observed anomaly, we re the accumulations not too deep (even taking 
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into consideration the void of 1250 x 10 emu/cc material above the 
stream). However, the gorges are not nearly as wide as the observed 
anomalies. Perhaps magnetic accumulations along the buried, possibly 
wider previous channels could be in part responsible for the highs asso“ 
dated with the present streams. 

It was frequently noted that fumarolic activity was particularly 
concentrated along the streams. Assuming that the present waterways are 
indeed above the previous channels, it is clear that they outline the ma- 
jor site of accumulation for both ground and surface waters. Thus the 
neighborhood of the streams could provide abundant steam to be expelled 
by the fumaroles and to aid in their transport of metallic constituents 
from the flow. The stream area would also be expected to undergo a 
greater amount of settling and fracturing since it is the thickest por- 
tion of the flow. Such conditions are favorable for the establishment 
of fumarolic vents. Welding would also be more likely to occur in the 
deeper sections of the flow; the degassing would further add to the 
fumarolic emanations. The chance of concentration and retention of mag- 
netic minerals along fumarole vents has been discussed. Since such accu- 
mulations are restricted to the immediate vicinity of the vents, and are 
infrequently preserved, it is doubtful that fumarolic magnetite could form 
a substantial zone the width of the observed magnetic anomalies. 
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The susceptibility of one of the samples of welded tuff was nearly 
-6 

2500 x 10 emu/cc. If the densely welded zone within the flow (s of 
this susceptibility, then the thickness of the flow required to produce 
the local anomalies is closer to that indicated by the seismic data* 
(Model IV of Table -4.7 lists the depth to the base of the flow assuming 
that the densely welded section is overlain by 30 meters of 1250 x 10 
emu/cc material.) it therefore seems possible that the welding process 
can double the apparent susceptibility of the flow. 

Another method of increasing the magnitude of the model anomalies 
Is to consider that there may be a substantial thickness of high suscep- 
tibility glacial fill beneath the flow. Fifty meters of 2500 x 10 ^ 
emu/cc material below a typical flow thickness would not sufficiently 
increase the model anomaly: if it is 250 meters wide and 55 meters deep, 

the increase equals almost 120 gammas; if 100 meters deep, the increase 
is only about 70 gammas. Even if the upper 10 meters of the drift had 
been thermally-altered during emplacement of the flow such that ft's sus- 
ceptibility had been enhanced to 5000 x 10 ^ emu/cc, the increase in 
the model anomaly would only be about another 50 gammas for a 50 meter 
thick flow underlain by glacial drift. In addition, there is no reason 
to expect either of these features to be restricted to the stream areas. 
A previous andesitic lava flow of similar susceptibility (2500 x 10 
emu/cc) and thickness might have been restricted to the stream channels, 
but there is no geologic evidence for such a stratum. 

In practice, the magnetic model studies tend to imply that the 
glacial drift below the flow is either relatively thin or of negligible 
apparent susceptib i 1 j ty . To model the larger anomalies requi res the 
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flow itself to be of higher susceptibility than 1250 x 10 emu/cc. 
Assuming higher susceptibilities for the welded zones is an effective 
means of • i nc re as i ng the value of the model anomalies to the observed 
magn i tudes . 

Were Fenner to examine the geophysical data, he would no doubt 
construe it as support for his hypothesis of the eruption. The un- 
dulating bedrock profiles would indicate the collapse of the old 
valley floor following the expulsion of the tuff flow from fissure 
feeders. The multiple anomalies across each branch would outline the 
portions of the flow originating from separate feeders. The increased 
magnetic anomalies in these areas would be due to compositional varia- 
tions outward from each vent and alterations within the glacial drift 
and bedrock through which the tuff was erupted. Since the valley 
floor would have been thinnest and therefore weakest along the stream 
channels, then the greatest alteration would be expected here, account- 
ing for the highest magnetic anomalies in this area. 

Assuming that a typical valley cross-section consists of 50 meters 

of tuff over 50 meters of drift, the total magnetic anomaly for a 250 
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meter width (of 1250 x 10 emu/cc) is A23 gammas. If the character 
of the emplacement of the flew had been so encompassing as to realign 
the magnetization vectors throughout the drift, Its apparent suscep- 
tibility might be about doubled, which could increase the anomaly by 
60 gammas. Assuming that the drift underwent extreme, alteration 
capable of doubling its susceptibility again, the anomaly could be 
increased by yet another 120 gammas. Had the Naknek bedrock been 
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similarly altered and its susceptibility enhanced to 5000 x 10 emu/cc 
throughout a 50 rreter section, the anomaly would be greater by almost 
300 gammas. Of course such extensive alteration is extremely unlikely, 

4,4 Concl us E on s 

The pyroclastic flow in the Valley of Ten Thousand Smokes is not 
a simple, uniform unit. Not only does it possess the vertical zonation 
typical of flows which have undergone welding deformation, but it is 
variable laterally as well. As would be expected, cross-sections of 
the flow become increasingly complex as one nears Novarupta. 

The presence of several individual anomalies along each profile 
suggests that some of the parameters of the flow are also variable 
laterally. Correlation of these anomalies among adjacent profiles 
within any one branch of the Valley shows that the caus i t i ve variations 
are more or less continuous, and are parallel to the axis of that branch 
of the V alley. This apparently linear nature of the anomalous magnetic 
bodies permits a specific type of two-dimensional analysis. Adaptation 
of this model analysis technique to the computer allows the calculation 
of several possible geologic models of the flow. Comparison of the 
anomalies produced by the models, with the observed anomalies, enables 
limitations of the probable flow parameters. 

The excellent "mi rror- i mage" correspondence of gravimetric and 
magnetic anomalies along adjacent traverses suggests that the primary 
cause of both anomalies is variable thickness of the low-density, high- - ' 
susceptibility pyroclastic material above the dense, low-susceptibility 
Naknek sedimentary bedrock. Since the top of the flow is essentially 
flat, th i s p ropes i t: 1 on req u i res t h a t the b as e o f the f 1 ow have con s i de rab 1 e 
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relief. In most cases, the deepest sections are indicated beneath the 
streams, implying that the present streams are above the pre-1912 
channels. Generally, the indicated base relief can be reconci led with 
a plausible pre-1912 valley configuration; however, in the southern 
branch, the undulating nature of the anomalies may outline faulting 
of the bedrock prior to the final emplacement of the pyroclastics in 
1912 . 

The sparse seismic data available to date reinforces the conclusion 
based on the other geophysical evidence that the base of the flow has 
high relief. Several intermediate layers within the flow are exposed 
by the refraction profiling. These seismic layers are probably caused 
by zones of differing degrees of welding. As would be expected, were 
Novarupta the eruptive source, the zones of dense welding become more 
prominent as one nears Novarupta, where the flow would have been the 
hottest; and there and elsewhere, the degree of welding is greater in 
the thicker sections of the flow. 

The inability of magnetic models (using thickness based on the 
gravi rnetri c , geomorpho! og i cal and seismic evidence, and a susceptibili- 
ty based on the average of the present studies) to produce the observed 
anomalies suggests that variable thickness is not solely responsible 
for the multiple magnetic anomalies encountered along each profile. It 
is necessary that the susceptibility also vary across the flow. Model 
studies suggest that the high anomalies may be due to increased sus- 
ceptibility in the welded section. The lower than average susceptibilities 
suggested for the borders of the flow are thought to result from weathered 
ash which has slumped from the adjacent steep valley walls. Fumarolic 
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leaching could also decrease the susceptibility of the flow, but its 
effects are thought to be restricted in depth and width. The longi- 
tudinal (from south to north) compositional change of the flow from 
more andesitic to more rhyolitic Is the probable cause of the lower 
susceptibilities indicated in the lower Valley. Pockets of this more 
rhyolitic ash near the eruptive vent, Novarupta, are thought to be 
responsible for the low susceptibilities suggested by the magnetic models 
for Novarupta Basin. 

The preparation of a model magnetic profile and Valley cross-section 
based on the ground magnetometer survey traverses alone Is unjustified. 
The complicated inter-relationships of changing thickness, composition 
and history (both weathering and fumarolic) throughout the flow precludes 
any such model profile. if used in conjunction wi th other geophysical 
data, however, the magnetometer data can help place limits on the various 
parameters of the flow. 

The magnetometer survey data are in agreement with the flow thick- 
nesses suggested by Curtis (19 68) , Kienle (1369, Model F) and Kienle’s 
recent seismic refraction profiles. In fact, the magnetometer data may 
indicate even greater thicknesses for some branches of the Valley. 
According to Curtis (1968, p. 207), such thicknesses show that the total 
volume of the eruption(s) in 1912 was great enough to equal the collapse 
of Katrnai crater, the subsidence surrounding Novarupta, and a substantial 
amount of other regional subsidence. 



CHAPTER V 


RECONNAISSANCE MAGNETOMETER SURVEY IN THE 
VICINITY OF NOVARUPTA’ VOLCANO 

5.1 Description of Novarupta 

Six mi les west and 780 meters below Katmai Crater is the volcanic 
dome, Novarupta. Present hypotheses consider Novarupta as the primary 
vent for the 1912 pyroclastic eruption. Supposedly this vent enlarged 
throughout the eruption, passed through a violent stage, and finally 
extruded a mass of viscous lava as a dome (Fenner, 1950b, p. 708 ). 

Other sources of the flow deposits are assumed to have been a series of 
eruptive fissures at the head of the Valley of Ten Thousand Smokes, 
with vents concentrated along a southeast-trending zone, approximately 
parallel to that marked by the line of Falling and Cerberus Mountains, 
and the domes at the base of Mount Mageik (Williams, 195*», p. 58), 

The dorrs of Novarupta is circular, approximately A00 meters in dia- 
meter and 91.5 meters high. It is surrounded by a high wall or "corona” 
of ejected pumice and glass blocks. The height of this crater rim rises 
continuously from a low point of 61 meters above the general level out- 
side the corona on the western side until it merges with the gouged-out 
face of Stumbling Mountain on the northeast, where Naknek sediments are 
exposed. Mild gas emission still occurs along the crest of the rim, 
and many of these fuma roles are far too hot to touch. 

The strata of the crater rim dip away from the dome at moderate 
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angles, but on the inward side are cut off sharply, presenting a steep 
face to the dome. The intersection of this face with the surface of 
the dome generally makes a V-shaped trough or moat, but locally the moat 
has a flat bottom. 

Novarupta dome itself consists of glassy though slightly porous 
lava. in forming the dome, the viscous lava congealed to a glass while 
it was still being thrust upward, resulting in a great portion of it 
being shattered. Much of this shattered material has remained on the 
slopes of the dome, although some blocks have dropped into the encircling 
moat . 

Flow banding is exhibited by the dome. On all sides the banding 
can be seen dipping toward the center of the dome (Curtis, I 968 , p. 192). 
Fenner attributes this orientation to the overturning of the outer lay- 
ers as the central mass was pushed up (Fenner, 1950, p. 7^9). The flow 
banding itself is due in part to alternations of glassy and cellular 
material but some Is due to alternations of light-gray rhyolitic glass 
and dark-brown andesitic scoria (Fenner, 1923, p. 56). The dark bands 
of andesitic composition are particularly common on the northeast side 
of the dome, while the central part is much closer to pure rhyolite 
(Curtis, I 968 , pp. 192, 194). At the summit of Novarupta dome is a 
trough (Fenner, 1925b, p, 219) ♦ 

5-2 Patterns of Subsidence Around Novarupta 

in general, the pattern of subsidence in the upper Valley is con- 
centric to Novarupta dome. Most of the faults in the vicinity of 
Novarupta seem to display simple vertical throws; others are step- like 
(Fenner, 1925b, p. 202). All the faults appear to record settling 
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toward Novarupta. Grtiben structures separate Stumbling Mountain from 
Broken Mountain. On a smaller scale, fissures intersect the edge of 
the dome itself. Thick depos i ts of pyroclastics conceal any offsets 
in Novarupta Basin. The prominent subsidence features are outlined in 
Fi gure 5.1. 

Curtis (I968) estimates that the circular scarp at a radius of 
approximately 2 km centered just northeast of Novarupta, has a downthrow 
towards the dome of at least 150 meters. The major structural features 
on the border of this scarp are the northern face of Falling Mountain, 
the faulted lower northern slope of Mt. Trident, Fenner Ridge, the 
faulted southeastern slope of Broken Mountain, and Greasy Pass. 

The sharp scarp face of Falling Mountain bears witness that a huge 
rock slide probab ly .preceded or accompanied the pyroclastic eruption. 

The irregularities of the landslide debris are not visible, as they are 
smoothed over by the pyroclastic deposits. Curtis (1988, p. I90) sup- 
poses that the northern part of Falling Mountain collapsed into a void 
as subsidence of Novarupta began. Fenner (1926, pp . 197" 198) reported 
a line of crater vents located at the base of the falling cliff which 
he attributed to a deep fissure. Many small fumaroles issued from the 
bare rock face of the mountain in the early years (Griggs., 1922, p. 242). 
Fenner (1920, p. 588) discovered that the fumarolic emanations were 
actively altering the igneous rock composing Falling Mountain, resulting 
in a loss of cohesive strength, and thus contributing to the numerous 
landslides which character i zed this mountain from the beginning. 

The faults which shattered the lower slopes of Mt , Trident were 
the scene of strong fumarolic activity shortly after the eruption 
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Figure 5*1* Flap of structural features in the vicinity of Nova- 
rupta Volcano. In all instances, the faults are 
downthrown towards the dome. Radials correspond to 
the major axes of the patterns of tephra distribu- 
tion given by Curtis (1968). Refer to Figures 5.2 
and 5.3b. Elevation contour interval is 500 feet. 
Based on aerial photos from Curtis ( 1 96 8 , Plate 5) 
and U.S.G.S., Katmai Quadrangle, aerial photos I 36 
and 137* 
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(Fenner, 1925b, p. 199). In 19&9, quiet steaming still occurred along 
the deep fissures along Fenner Ridge and the slump blocks on the 
southern slopes of Broken Mountain. Throws of almost 30 meters have 
been recorded for the larger faults on Broken Mountain (Fenner, 1925b, 
p. 217). And, Fenner (1950b, p, 716) has estimated the total downthrown 
displacement of these faults to be about 90 meters, but the total 
bedrock displacement is probably much greater. The undisturbed appear- 
ance of the pumice beds overlying these faults prove that the present 
observable offsets record only gradual readjustment after the deposition 
of the pyrocl as tics . The initial fracturing could have been violent 
and more extensive than the present exposures indicate (Fenner, 1925b, 
p. 218). Even without allowing for burial of its base beneath the 
pyroclastic flow, Greasy Pass is over 90 meters above the top of the 
f 1 cm in Novarupta Basin. 

The major circular scarp just described is very prominent in aerial 
photos of the head of the Valley. There are three oval depressions with 
in the scarp which are also quite conspicuous in aerial photos. The 
outline of each of these resembles that of the crater rim of Novarupta 
dome, which is the most obvious of the oval features. The major axis 
of the ovals is about 2 km and they all intersect to the northeast of 
the dome near the center for the outer circular scarp. 

The crater of Novarupta is the clearest of the oval features; a 
second oval depression Is offset to the north, barely intersecting 
the lava dome, and most visible where it cuts the crater rim of 
Novarupta and breaks the summit of Stumbling Mountain. To a viewer 
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standing on Novarupta dome, the cuts in the rim would appear as “radial" 
fissures. They are the site of continuing fumarolic activity. Tri- 
dent Basin, located a little southeast of the dome, also has this same 
general • out! i ne . The basin heads at Fenner Ridge on the east and is 
separated from Novarupta Basin by a slight rise to the west. Apparently 
surface runoff was disrupted by the sinking of Trident Basin; several 
small lakes are now situated in this basin. Reportedly, a stream from 
this drainage area has cut a gorge through the rise and into Novarupta 
Basin. This last feature implies that the stream was able to cut through 
the pyroclast i cs while the basin was sinking. 

Each of these small collapse ovals could outline the top of the 
Novarupta vent during a major eruption stage. Then the location of 
an oval would be directly related to a specific conduit orientation. 

The sequence of orientations is implied by the "freshness" of the collapse 
The crater r?m (oval III) is clearly the most recent. Trident Basin 
(oval l) , although recording substantial subsidence, has no distinct 
fissures like those of oval II which cut the corona and Stumbling Moun- 
tain. Thus it appears that Trident Basin (oval l) predates oval II, 

which predates the crater rim (oval Ml). 

If the collapse ovals do record previous orientations of the Nova- 
rupta conduit, they should correlate with the isopachus tephra trends 
as. measured and contoured by Curtis 0968) (refer to Figure 5.2). 

Rad 1 a 1 s from- Novarupta along the main axis of the tephra contours are 
superimposed on the subsidence features in Figure 5.1. The first layer 

Its distribution is elongate over oval l. 


contoured by Curtis is C. 



139 



Figure 5.2. Isopachs of tephra layers C, D, F 

and G. (After Curtis, 1968, Figures 
8, 11, 12 and 13). Refer to 
Figure 5.1. 
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Radials for layers 0 and F actually correspond to the latest oval, No. 
ill. Layer G, from the last major eruption of Novarupta (Curtis, 1 9 68 , 
p. 198), is strongly directed northward over oval II. Thus there is 
some evidence that Novarupta' s conduit underwent changes of position 
throughout the eruption. Some of these previous positions have been 
preserved in oval features intersecting the lava dome. 

Beyond the major circular scarp are several arcuate faults on 
Broken Mountain and on the lower slopes of Mt. Trident, and the 
parallel fumarolic lineations near the mouth of Novarupta Basin. Q.uiet 
steaming characterizes many of the concentric and arcuate features. 

Surely these denote subsurface fractures related to the family of 
concentric fractures about Novarupta. 

The oval patterns of subsidence close to the dome probably record 

t 

the sequence and direction of the changes in the orientation of the 
eruptive vent of Novarupta. The larger circular scarp was probably 
formed as the overall area collapsed in response to the removal of 
vast amounts of magma from the underlying reservoir. The oval features 
and surrounding concentric and arcuate fractures therefore probably 
outline an underlying intrusive body related to the erupted pyroclastics 
and the extrusive dome of Novarupta. The presence of this cooling mass 
is also indicated by the continued fumarolic activity, which is restricted 
to the head of the Valley. 

5.3 Magnetometer Survey Around Novarupta 

It was hoped that magnetometry could be used to delineate the ig- 
neous subsurface at the head of the Valley. To test the applicability 
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of this method of detailed subsurface mapping within the ground magneto- 
meter regime of the Valley, several reconnas i ssance ground magnetometer 
traverses were conducted around Novarupta. Most of these traverses were 
oriented radially to the dome. In addition, a few spot readings were 
taken in the immediate vicinity. All. work was accompli shed wi th the 
vertical field magnetometer. The locations of these stations and tra- 
verses are shown in Figure 5»3a. 

A contour map of the anomalous vertical field in the vicinity of 
Novarupta ts presented in Figure 5.3b. A 500 gamma contour interval 
shows well the concentric pattern of the anomaly about a high centered 
just northeast of the dome itself. Two interesting features of this 
anomaly are its gentle horizontal gradients and insignificant distor- 
tion by local topography. These features indicate that the source is 
both deep and thick. As for the shape of the anomaly: near the dome, 
the concentric oval contours correspond beautifully with the outline of 
the crater rim; farther out, the contours are elongate more down valley, 
and are similar in shape to the outer collapse oval of Figure 5-1. I n ~ 
spection of total field anomalies along 1970 traverses N-N" and N'"-n 
reveals that the contours extending off the left of Figure 5. 3b curve 
back around to the south, then southeast, essentially parallel ing- the 
500 gamma contour in that quadrant of the figure. 

The elongation of the contours and the outer oval down-valley may 
indicate structural control. It is along this same bearing that Curtis 
defines the axis of the Naknek anticline, which is slightly offset from 
the volcanic line (see Figure 1.1). Williams also felt that the eruptive 
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fissures were aligned along this trend. The elongation may overlie the 
bedrock contact of Naknek sediments to the north from volcanics to the 
south. Or, the elongation may merely reflect the control of confining 
mountains bordering both sides of Novarupta Basin, 

A more extensive survey will be required to follow the magnetic 
anomalies westward beyond Falling Mountain and past Mageik Basin. Such 
exploration could reveal a relationship of these collapse areas and the 
hypothesized bedrock contact. A more encompassing survey might also 
yield more information about possible contemporaneous settling of the 
Baked-B roken Mountain complex. 

The magnetic iow just northeast of the dome in the center of the 
otherwise inwardly increasing anomaly is based on two spot readings 
only; one on the summit of Novarupta, the other on the gouged-out western 
face of Stumbling Mountain. Perhaps the abrupt terrain affected the 
recorded field. However, Anma (1972) observed the same low as he crossed 
the dome in his aeromagnetic survey. This inner low may indicate that 
some central portion of the Novarupta intrusive is still near the 
Curie temperature of the magnetic minerals, since the magnetic suscep- 
tibility of volcanic rocks is drastically reduced as this temperature 
Is approached. The continued fumarolic activity in this area only is 
proof that considerable heat is still present in this subsurface. But 
since rocks are above their Curie temperatures long before they are 
molten, and the affected volume of rock in this instance is probably 
small and also shallow, one would not expect detection of this zone 
of anomalous temperatures by the seismic method of Berg and Kubota (I 967 ) 



A more plausible explanation for this central low may be based on 
the composition of the dome: the central zone of Novarupta is of al- 

most pure rhyolite* whereas the margin has been contaminated by streaks 

of more basic material. Susceptibility studies show that the banded 
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lava has a susceptibility of 1085 x 10 emu/cc, compared to only 30 ^ 
x 10 ^ emu/cc for the rhyolite. The fact that the larger conical 
anomaly centered on Novarupta seems to be caused by higher susceptibility 
material supports Forbes 1 hypothesis that magmatic differentiation was 
responsible for the variable ejecta of the 1912 eruption. Accordingly, 
as the rhyolitic magma was erupted by Novarupta, the remainder of the 
magma in the reservoir would have become increasingly basic; consequently 
it had the potential of attaining relatively higher susceptibilities. 

The banding observed at the margin of the dome suggests that some of this 
h i gher suscept i b i 1 i ty material was dragged up along the vent contacts 
during the expulsion of the rhyolitic plug at the end of the eruptive 
stage of Novarupta. The small concentrat i on of low-susceptibility rhyo- 
lite forming the central dome of Novarupta as contrasted to the mass of 
higher suscept ibi 1 i ty material remaining in the subsurface can account 
for the magnetic low observed directly over the dome. 

Anma (1972) has suggested that this central low is caused by adja- 
cent dipoles underlying Novarupta and Stumbling Mountain (which would 
correspond to collapse ovals l! and III). However, the complete en- 
circlement of high values about the rim of Novarupta 1 s crater is con- 
trary to his model and tends to discredit this hypothesis. 

The local high anomaly directly west of the dome corresponds to a 
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zone of previously intense fumarolic activity. The local high at the 
northern end of Slippery Pass is marked by a large area of quiet steam- 
ing and bright, multicolored clay, bounded on the east by a sizeable 
fissure. A magnetic low is also indicated along the ’'valley" extending 
northeast of the dome from Novarupta Basin to Broken Mountain. 

The individual profiles in the Novarupta region are jagged with 
minor peaks and troughs imposed on the overall dome-shaped anomaly. 

There seems to be some correlative value of anomalies between adjacent 
radial traverses (refer to Figure 5.**). Such correlation indicates 
near-surface features concentric about Novarupta, Fractures in the 
py rod as t i cs , and possibly fumarolic alterations along these, could 
account for the narrow, shallow anomalies. 

Broader breaks in slope along the magnetic pprofiles, like those 
clearly exhibited along traverses G ,M -g and l-l 1 , could indicate bedrock 
structure. Major subsidence faults could be responsible for these off- 
sets . 

No evidence for fractures radial to Novarupta were revealed by this 
survey. However, the nature of magnetic profiling precludes the de- 
tection of anomalies paralleling traverses. 

5.4 Anomalous Geophysical Parameters In Novarupta Bas-in 

Adjacent total magnetic anomaly profiles across the lower half of 
Novarupta Basin (E-E 1 and N’tN") possess a similarity (Figures 4.8 and 
4.12) which suggests that multiple, linear magnetic bodies paralleling 
the axis of the basin are their sources. 

Model studies (refer to Table 4.6) show that the material causing 
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Figure 5.4. 


Vertical magnetic anomaly profiles radial to Novarupta 
Dome. Refer to Figure 5.3a for locations. 
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several of the local anomalies must be of quite low susceptibility, even 
if the flow is only 100 meters thick as has been indicated by seismic 
measurements (Table 1.3). It was proposed (Section 4.3) that pockets 
of predominately rhyolitic material might be responsible for the small 
local anomalies. 

Another unusual property of these magnetic profiles is the abrupt 
change in the regional trend about half-way across the basin (refer to 
Figures 4.8, 4.12 and 4.13). Perhaps this feature reflects the encroach- 
ment of a wedge of volcanlcs related to the volcanic peaks of the Aleu- 
tian Range into the Naknek sedimentary province. 

Yet another conspicuous quality of the magnetic profiles in the 
lower half of Novarupta Basin is their negative values relative to the 
overall regional geomagnetic field. It has been suggested that this 
magnetical ly- low zone may represent the low of the dipole field due to 
Novarupta, perhaps including contributions from the other neighboring 
igneous masses, too; but the magnitude is too great to substantiate 
this proposition, and the alignment of the anomaly is not that expected 
of a dipole at Novarupta dome. Another remotely possible explanation 
for this zone would be the presence of a reversely-magnetized body 
beneath the flow in this area. 

Anomalous lows in other parameters have been recorded in . this same 
portion of Novarupta Basin. Zeis and Allen (1923) measured relatively 
low fumarolic temperatures in this zone in 1919. Their data has been 
contoured and is shown in Figure 5.5* There are too few data on the 
composition of the early emanations from the fumaroles in this particular 
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Figure 5.5a. : Location of fumaroles studied by Allen and Zies In 1919 
(Zies, 1929 » pp. 2-3)- Contour interval is 1000 feet. 

Figure 5.5b. Isothermal map of fumarole vapor temperatures as measured 
by Allen and Zies in 1919 (1923, PP- 104-106). Contour 
interval i s 1 00°F. 
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area of the Basin to permit any predictions as to their most probable 
primary source. Perhaps the gasses vented here were of a secondary 
nature, in contrast to the fumes high in CC^ and the sulphur gasses 
which characterized the high temperature vents at the head of the 
Valley, and were indicative of connections with a deep seated source 
(Lovering, 1957, p. 1588), 

The area of the relatively low fumarolic temperatures observed 
in 1919 lies beyond the prominent collapse features in Novarupta Basin 
(beyond the larger collapse oval shown in Figure 5.1). if this area 
is located beyond the bedrock faults related to the deflation of an- 
cestral Novarupta, it would not have been as accessible to the heat 
and vapors from the deep seated source as were the fumaroles within 
the disturbed areas closer to Novarupta. The same isolation would hold 
true were this area between widely spaced bedrock fractures; another 
major fault caused by the subsidence of the Baked Mountain-Broken Moun- 
tain-Novarupta complex is thought to be aligned with the western face 
of Baked Mountain, which would be to the west of this zone. 

In this same portion of Novarupta Basin, Klenle (19&9, PP . 126- 
127) took gravity measurements. He found the regional field in the basin 
to be quite different from that encountered in the other branches of 
the Valley. The regional gradient which plunges almost 7 mg a 1 s from 
Baked Mountain to Falling Mountain is almost twice the regional gradient 
observed elsewhere in the Valley, and the field is some 6 mgals more 
negative. 

In all probabi 1 i ty , the distinctly negative nature of both the 
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magnetic and gravity fields in Novarupta Basin can be explained in terms 
of regional tectonics. The alignment of these features might be struc- 
turally controlled by a deep seated fault with the northern side down- 
thrown by some 3 km relative to the southern side (Kienle, 1369, p. 23). 
The Naknek section also thickens as one approaches the range from the 
north. Magma reservoirs have been located at relatively shallow depths 
below the volcanic peaks (Kuboto and Berg, 1967). All three of these 
structural possibilities are capable of reducing the observed gravity 
and magnetic fields because they increase the amount of low-density, 
low-susceptibility material between the stations on the surface and the 
dense, high-susceptibility volcanic basement. These arguments are parti- 
cularly applicable if the lower portion of Novarupta Basin is a relatively 
undisturbed block within the area-wide subsidence towards Novarupta. 

Unfortunately the seismic profiles are not complete in this area; 
however, seismic profile D-D' (refer to Appendix D) shews the flow to 
thin from 95 meters to 80 meters from the base of Novarupta towards the 
center of the basin. Since the top of the flow is relatively level 
along this profile, the data suggest that the bedrock dips towards Nova- 
rupta. Such an interpretation is in keeping with the often recorded 
pattern of subsidence towards Novarupta dome. The flow thickness 
suggested by Kienle ( 1969 , Model F) for gravity profile KV-4 seems 
unusually low, only 50 meters (refer to Table 1.3). The lack of bedrock 
stations for reference may have led to an erroneous value; or, since 
the tuff is considerably denser here near its source than farther down 
valley (refer to Appendix C) it would make a lower density contrast 
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with the bedrock than assumed in Model F (a lower density contrast would 
have the effect of increasing the gravimetri cal ly-indi cated thickness). 
Kienle's 1970 seismic refraction profile D-D' shows the flow to be about 
80 meters thick near this area. Profile 3 of Sbar and Matumoto is very 
erratic but it does suggest that the flow is from 40 to 70 meters thick 
along the same line as Kienle‘s profile (refer to Appendix D and Figure 

1.3). 

If the bedrock here is actually higher relative to the general 
subsidence closer to Novarupta, then the nearly level surface of the 
flow would require that the pyroclastics are thinner in this anomalous 
area than on either side. Such an interpretation would also be in keeping 
with the relatively small magnitudes of the smoothed local magnetic 
anomal i es . 

5.5 Conclusions 

Preliminary magnetometry in the vicinity of Novarupta Volcano 
indicates the presence of a broad, dome-like structure underlying the 
head of the Valley. The nature of the major anomaly indicates that 
the causitive body is both deep and thick. Generally, the magnetic 
anomalies decrease radially from the dome. The contoured anomalies 
form concentric ovals centered just northeast of the dome itself. 

The direction of elongation of these ovals parallels the Aleutian Range, 
and is probably related to the overall structure of the area. The 
inconsistent central low is caused by the small concentration of low 
susceptibility rhyolite forming the central zone of Novarupta dome. 

Minor anomalies superimposed on the main anomaly imply bedrock faulting 
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as well as fracturing of the overlying pyroclastics and possible fuma- 
rol ic al teration. 

Oval topographic depressions intersecting Novarupta dome may pre- 
serve a record of conduit orientations. Alteration in orientation 
of the vent at Novarupta could account for the varied patterns of dis- 
tribution displayed by successive tephra layers. Partial agreement 
of extensions from Novarupta along the direction of elongation of the 
tephra contours with the major axis of the collapse ovals supports 
such an hypothesis. 

The radial orientation of the traverses precluded detection of 
radial features. The survey can therefore give no evidence of a sug- 
gested relict conduit connecting Mounts Katmai and Trident with Nova- 
rupta. 

The impressive agreement of the shapes of the magnetic contours 
and the subsidence features indicates that both outline the intrusive 
body beneath. A more comprehensive magnetometer study at the head of 
the Valley could produce a detailed map of the Intrusive body which 
is expressed at the surface by Novarupta and its concentric fault 
system. Aerial surveying would reduce the complication caused by near 
surface anomalies. 

The unusually negative values recorded for both gravity and magnetic 
fields in the lower half of Novarupta Basin are related to their proxi- 
mity to the volcanic axis of the Alaska Peninsula. Increased Naknek 
thicknesses and shallow magma reservoirs are probably the principal 
causes of the regional trends adjacent to the volcanic peaks. The 
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anomalously low fumarollc temperatures in this area in 1919 may mean 
that this region did not have the connections v^i th the deep-seated 
source responsible for the hotter emanations observed closer to Nova- 
rupta in a st ructural ly-more disturbed area. It is proposed that the 
lower half of Novarupta Basin may be a relatively upthrust block be- 
tween two major faults: one caused by the settling encompassing the 

entire Baked Mountain-Broken Moun tai n-Novarupta complex, and the other 
by the more localized collapse around Novarupta itself. 

The amount of subsidence recorded in faults about Novarupta can 
account for an ancestral Novarupta of sufficient altitude to allow 
deposition of pyroclastics in all areas where they are to be found. 
Spurr mapped a high, conical structure situated at the head of the 
Valley in I898, It is proposed that this mountain was ancestral 
Novarupta. The present Novarupta volcano is the last remnant of 
this mountain, and the arcuate and circular collapse features preserve 
the pattern of its deflation. 



CHAPTER VI 


SUMMARY OF THE GROUND MAGNETOMETER SURVEY 

Magnetic susceptibilities were determined for many samples 
collected from the Valley of Ten Thousand Smokes. It was found that 
the pyroclastic flow, the air-fall pyroclastlcs , and the banded lava 
and pumice all have susceptibilities near 1250 x 10 emii/cc. The 

-6 

average susceptibility of the andesitic lava is close to 2500 x 10 

emu/cc, but the susceptibility of the rhyolitic material is an order 

of magnitude smaller. Naknek sediments normally have susceptibilities 

less than 50 x 10 ^ emu/cc, but if subjected to extreme heat and 

volcanic emanations, the altered sediments can attain susceptibilities 

of 5000 x 10 emu/cc. It is difficult to predict the susceptibility 

of the glacial drift, as none was sampled. Probably the drift is 

primarily composed of randomly-oriented, andes i tic debris . Its 

-6 

susceptibility is most likely less than 1250 x 10 emu/cc. 

A small scale total field magnetometer survey was conducted over 
a zone of relict fumaroles near the terminus of Broken Mountain Valley 
Including both crater vents and filled fissure vents. It was found 
that the large, narrow total magnetic anomalies sh ow spectacular 
agreement with the surficial geology. The close association of the 
anomalies with the fumarolic markings is a strong indication that 
pockets of magnetic minerals have been preserved along the vents at 
shallow depths. Sharp anomalies were often encountered near extinct 
fumaroles in all branches of the Valley. It is thus proposed that the 
preservation of such accumulations of magnetic minerals along fumarolic 
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vents Is common throughout the Valley. 

In general, the vertical magnetic anomalies in the vicinity of 
Novarupta Volcano increase as one approaches a point just northeast 
of the dome. The contoured vertical anomaly Is in fact composed of 
nearly concentric ovals. Evidently, the dome is simply the surficial 
expression of a sizeable, conical-shaped intrusive, probably of andesitic 
composition. The inconsistant low in the center of the Novarupta 
magnetic anomaly is probably caused by the small concentration of 
rhyolitic lava forming the central zone of the dome. Corresponding off- 
sets in the magnetic anomalies among adjacent traverses radial to 
Novarupta suggest that there is concentric faulting about the dome. 

This faulting apparently includes fractures in the bedrock as well as 
in the pyroclastic flow, and would be related to the general pattern 
of subsidence around Novarupta. 

As one approaches the volcanic peaks at the head of the Valley, 
the magnetic anomalies tend to be increasingly more negative. Increased 
thickness of the Naknek sediments as well as shallow magma chambers are 
no doubt the principal causes of this regional trend. 

Cross-valley magnetometer traverses were made In each branch of 
the Valley. The pattern of individual anomalies along adjacent, 
nearly-parallel traverses is very similar. Evidently, the magnetic 
bodies responsible for the anomalies are linear, and trend parallel 
to the axis of each branch of the Valley. 

The mirror-image correspondence of magnetic and gravimetric 
anomalies along nearby cross-valley profiles suggests that the primary 
cause of both anomalies is variable thickness of high-susceptibility, 
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low-density pyroclastics above the low-susceptibility, high-density 
Naknek sedimentary bedrock. Since the top of the flow is relatively 
flat, this proposition requires that the base of the flow must have 
considerable relief. In most cases, the implied base relief can be 
reconciled with a plausible pre-1912 valley configuration. Since 
the deepest sections are usually indicated beneath the present 
streams, it follows that the streams have been established above 
the p re- 1912 channels. The wavy pattern of the anomalies in the 
southern branch of the Valley probably outline major bedrock fault- 
ing which probably occurred prior to the final emplacement of the 
pyroclastic flow. 

Using thicknesses based on the gravimetric, seismic and geo- 
morphological evidence, and a susceptibility based on the average 
for the pyroclastics sampled for this study, it is impossible to 
model the observed anomalies. Clearly variable thickness is not 
the single cause of the multiple magnetic anomalies encountered 
along the cross-valley traverses, it is necessary that the suscepti- 
bility of the flow also varies laterally. The susceptibility studies 
do suggest that the denser welded tuff can possess much higher 
susceptibilities than were observed for the unconsolidated airfall 
pyroclastics. Welding is also more prominent in the thicker sections 
of the flow. If the welded portion of the flow in the thicker 
sections of the flow is of higher susceptibility than the average 
measured for the pyroclastic flow to date, then reasonable thicknesses 
can be used to model the flow. To model the small anomalies over the 
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sections of the flow near the valley walls often requires lower than 
average susceptibilities. It is thought that the slumping of weathered 
ash from the adjacent steep slopes could be responsible for the lower 
susceptibilities. The highly rhyolitic composition of the flow in 
Novarupta Basin and in the lower valley has been cited as the cause 
of the low susceptibilities suggested for these areas by the magnetic 
model studies. 

in conclusion, if used in conjunction with other geological and 
geophysical data, the magnetometer data can. set limits on estimates 
of the thickness and composition of the pyroclastic flow. 
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APPENDIX A 

REDUCTION OF MAGNETOMETER SURVEY DATA 
A.l The Reduction Technique 

Reduction of ground magnetic survey data to magnetic anomaly 
values along metrically-described traverses involves data selection, 
interpolation, conversion and correlation. 

Ground survey data include the magnetometer reading(s), surface 
station spacing, time of occupation, altimeter reading, and field 
notes. A unique magnetometer reading must be selected for every 
station where multiple readings were made. Because time and eleva- 
tion were only noted intermittently, values of these parameters had 
to be interpolated and assigned to each station. These interpolated 
values were then utilized in correcting for diurnal variations and 
for calculating the horizontal station spacing, respectively. 

a) Reduction of Total Field Data 

The Lamor frequency, L, as recorded by the Rubidium vapor 
magnetometer at the magnetometer base, can be converted to gammas of 
total field strength, R, by applying: 



A. 667 x 10 


where 4.667 is the calibration constant of the instrument. 

A base datum of 53,000 gammas was chosen to reduce each Rubidium 
vapor magnetometer reading to a pos i ti ve fluctuation of the geomagnetic 
field. This diurnal correction, B, corresponding to each reading 
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of L , is g i ven by : 

B * 5300 - R. 

Values of diurnal correction, B, as calculated from the geomagnetic 
variation as measured at the magnetometer base are presented in 
Figure A. 1 . 

The Rubidium vapor magnetometer was inoperative during the 
surveys conducted on July 13 and July 14, 1969. From a composite 
of all other diurnal variation curves (field versus time) a rate 
of daily geomagnetic fluctuation was assumed. Commonly, the field 
decreased at the rate of 9.8 gammas per hour in the mornings until 
1300 local time; then it increased at the rate of 10.9 gammas per 
hour until evening. Since the base station was occupied with the 
Elsec magnetometer at the start and conclusion of these survey days, 
it was a simple task to impose these common fluctuation rates to a 
datum dictated by the Elsec base station data. This was a procedure 
followed In applying a diurnal correction for magnetometer traverses 
A-A‘ and B-B' . 

In reference to the above, it should be noted that for simulta- 
neous readings at the magnetometer base, the Elsec magnetometer 
recorded a total field 771 garranas greater than that recorded by the 
Varian magnetometer (see Table A.l). This difference is assumed to 
result from the difference in the location of the proton-precession 
magnetometer over the wooden post and the Rubidium vapor magnetometer 
sensing head burled a few feet away and an uncalibrated crystal 
oscillator in the Rubidium magnetometer. In practice, this discrepancy 
Is of 1 little consequence as it is the variation which is of interest 
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Table A . 1 


Simultaneous Magnetometer Readings at the 
Magnetometer Base Station 


Date 

Time 

Elsec Reading 
(gammas) 

Varian Reading 
(gammas) 

Difference 
of columns 3&4 

7-7-69 

9:30 

53701 

52930 

771 

7-7-69 

15:20 

53690 

52919 

771 

7-9-69 

16:42 

53726 

52949 

774 

7-9-69 

17:26 

53728 

52908 

771 

7-10-69 

19:35 

53731 

52952 

782 

7-11-69 

1 1:30 

53678 

52957 

770 


(the average deviation of the difference from 771 gammas is 2.5 gammas} 
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and the 771 gamma difference was constant over the whole survey period. 

It does, however, require that a base datum of 53,771 gammas be used 
for July 13 and 14, so that the survey readings on these days will be 
corrected to the same reference as the other traverses; |.e., in these 
cases, the diurnal correction is found by: 

B - F + 53771 “ E 

where 1 F * is the amount of fluctuation incurred at the "common" rate 

1) since the time of the initial Elsec base reading before 1300, or 

2) leading to the final Elsec base reading after 1300. The time and 
value of the Elsec base readings determined the locus of the F values 
by specifying a point on each rate line. To effect the same datum for 
all surveys, i.e., to cause all survey data to be corrected relative 
to a singular magnitude for the geomagnetic field at the magnetometer 
base station, a constant must be added to the survey data which have 

no corresponding Rubidium vapor base data. This constant is determined 
by the difference between the Elsec base reading, E, and 53771 gammas. 

The regional geomagnetic field in the Valley as interpolated 
from maps of the earth's field in Alaska (Deel, 19W is approximately 
53,800 gammas with inclination 71.2 degrees and declination 22.5 degrees. 
This mean field for the general area was used to reduce survey values 
to anomalies in the total field. 

The proton precession (Elsec) magnetometer was commonly read more 
than once at each station to eliminate electronic errors and short term 
fluctuations in the earth's field. The method used to select the single 
value used in further computations was (in order of preference): 



a) a repeated reading; 

b) the nearest even integer to the arithmetic mean of clustered 
values; and 

c) the nearest even integer to the arithmetic mean of those 
multiple measurements differing by less than ten instrumental units. 
(Statistically, using only even integers reduces the chance of biasing 
the overal 1 data . ) 

Proton precession, P as displayed on the Elsec instrument, can 
be converted to gammas of magnetic field strength, S, using the 
instrument calibration constant, so that: 

2.4051 x 10 9 

b " P 

The diurnal geomagnetic variations were removed from each of 
the base magnetometer readings with respect to time of survey station 
occupation. If B' represents the diurnal correction interpolated to 
the time of occupation of a station of Elsec reading, S, then the 
value S* of the magnetic field at that station (with diurnal effects 
removed) is given by: 

S* = S + B‘ . 

The value of the total anomaly (with respect to a regional field 
of 53,800 gammas) at this station would then be found by: 

TA - S 1 - 53800 . 

Or in a complete form, the reduction of survey data to total 
anomalies could be accomplished by: 

TA = R‘ - S + 800 , 

where fV represents the total field at base as recorded by the Rubidium 
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vapor magnetometer, interpolated to the time of occupation of a survey 
station with field reading S. 

A sample of the computer program designed to reduce total-field 
magnetometer survey data is given in Appendix A2. 
b) Reduction of Vertical Field Data 

As discussed in Chapter 1, the fluxgate magnetometer displays 
differences in the vertical field directly in gammas. Obviously no 
conversion factor is required. 

Removal of diurnal effects from vertical magnetometer survey 
data is treated differently from that for the Elsec field data. 

Repeated base readings with the fluxgate instrument often showed a 
change in the vertical field opposite in sign to the total diurnal 
variation, presumably due to instrument drift. Field stations were 
not systematically reoccupied because the object of the fluxgate 
surveys was to obtain a rapid, general view of the magnetic situation, 
it is therefore impossible to determine a detailed correction for the 
combined effects of instrumental and diurnal geomagnetic variations. 

The vertical field traverse data have been adjusted for diurnal 
drift (instrumental and geomagnetic) by linear interpolation of the 
difference in base readings with respect to times of station occupa- 
tion. The vertical survey readings were further adjusted to give a 
unique value per location In the case of a station common to two or 
more traverses. 

Using 53,771 gammas as the total field at the base station, with 
an inclination of 71.2 degrees, the absolute vertical field at the 
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base would be 50899 gammas. However, to obtain the vertical anomaly, 
VA, one must remove the regional vertical field of 50927 gammas 
(= 53800 x sin 71.2). 

One must also remember that the fluxgate Instrument was set to 
•H 50 gammas at the magnetometer base station at the onset of each 
survey day. Thus, the vertical survey reading, V, which has been 
corrected for diurnal variations, can be reduced to the vertical 
anomaly by: 


VA = V - 177 . 
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APPENDIX A. 2 


C SAMPLE PROGRAM FOR THE REDUCTION OF GROUND MAGNETOMETER DATA 
C TO RESIDUAL ANOMALIES IN TOTAL FlELO 

C. RRFRAREO IN 1070 BY MARLA CAVE TRXU.E 

DIMENSION T( 2001 »Ut?O01 , t < 200 I, XX I 200) , YY ( 2001 ,VG(200 J » . 
i$h( 2*0 >,su 2001 ,r ( tso) ,rh ( sm ,p.m(50I ,exmo> ,txi ltn 
c STATE TOTAL DATA POINTS- N= TOT AL STATIONS ALONG TRAVERSE 
r M= T CJTAL DIURNAL FIELD MEASUREMENTS 

c NK=TOTAL BASE ALTIMETER READINGS 

1 RfAOI 1 «?, END=90 1 N,M,NW 

2 fUMMAT l~M 14,2X1) 

I \ ° U T FIELD DATA- ST AT T ON* t SEPARATED FROM PREVIOUS STATION. 

BY D=C ABLE LENGTHS t YY= ESTIMATED 
EFCT, AT TIME- TsHHURS f. XX=MTNUTFS, 
ALTJMMER KLAI1ING- l = f e l t , magnetometer 
READINGS SH = AT HIGH PROBE POSITION 
AND St. = A T LOW POSITION. 

00 4 t-t.N 

3 FORMAT 1?<F2. 0,2X1 , F4 . 0 , 2 X , F4 . 0 , 2 X , F? . 0 , ?X , 2 1 f 1,2X1 I 

4 READ (1,3) T(IJ ,XXII) ,F(II ,nU) ,YYU) ,SH( II ,SU I I 

C INPUT base MAGNETOMETER RFAOINGS=R. AT TIME .RH*HOURS,C RM- 
c Ml Null S 

DO 6 1=1, M 

5 FORMAT ( P6.1 .21 2X,F2.0> 1 
u A l ( I . ’• ) R .KM, R M 

C INPUT BASF ALTIMETER M t AO I NG S= EX FEET AT TX=HOURS t 
C VG*MI NUTES 

00 A T'l.NM 

. r format n-4.a,2(2x,F2.o> > 

3 READ <1,71 E X , T X , VG 

C TO CONVERT DATA TO METERS, DECIMAL HOURS, GAMMAS 
O' 1 10 I = t , V 

on 1 = t . 0 4*pm ) ♦( YYT I )*0. 304H01) 

T t I > = TU IMXXtl) /60.) 

E(II=tlI 1*0. 304301 
S h ( t 1 = 2.4031 EO/SHt I ) 

!0 SL t n=2.405lF9/SH( SI 

DO ! 2 I = I * M ‘ 

ftHf I )*RH[ I 1 ♦ ( 0 M ( 11/60. ) 

] ? Rt I >=R I I ) /4. 667E-3 
DO 14 1=1, NM 
rx t I 1 =TX1 I 1 + I VG< I > /60. ) 

14 F.x m = EXl I 1*0.304301 
C TO INTERPOLATE TIME FOR EACH STATION 
DO 20 I = 1 , M 
xx ( l ) *0. 


20 YY I I 1=0. 

00 25 1=1, N 

IE ITU) .FO.O.l GO TO 26 
YY ( I 1 = D ( I ) 

xxm = Tm 

25 CONTINUE / 

CALL YLNTPL < YY, XX , 0 , T , N) 

C SUBROUTINE YLNTPL DOES LINEAR INTERPOLATION BETWEEN KNOWN VALUES 
C TO REDUCE ALTIMETER PAT A TO BASE DATUM = 2550 FfET 
DU 10 1*1, NM 

30 r X [ I 1 *771 . 21 -E X< 1 1 
DU 3 ? 1 = 1 , N 

32 vr.r i > = o. 

CALI. rLNIPL ( TX, ex,T,VG,NJ 
00 3 5 I = 1 , N 
35 E('n=Ei n+vam 

C TO 1NTT Mi,<n Ml CirvATIMN FOR EACH STATION 
DO 40 I = 1 , N 
XXI I 1 = 0. 

40 YYl I 1=0. 

on 4 S ! = t,N 

IE I E ( II.EO.VGdn GO TO 45 
YY [ I 1 =EI J } 

XX I 1 I - 1 t 1 ) 

45 CONTINUE 

C TO CONVERT SURFACE SPACING INTO HORIZONTAL DISTANCE FROM 
C THE INITIAL STATION 

XX ( l 1 =0. 

om=o. 

DO 50 1=2, N 

XX(ll«SO*T (A3$( ( Dill o*2)-C CEL I- 11-e< 111**211) 

50 ix n *di i-i i mxtn 

C TO CORRECT MAGNETOMETER DATA TOR DIURNAL VARIATIONS 
00 60 1=1, M 
60 R ( 1 >=5.3000. -Rt I 1 
DO 62 1=1, N 
62 VG< 11=0. 

CALL YLNTPL ( RH, R , T , VG, N > 

DO 65 T = 1 ,N 
SHI I ) =SH( 1 ) + VG ( T 1 
65 St. m*SL II I + VG( II 

C TO RFOIJCF DATA TO ANOMALIES IN REGIONAL FIELO*53BOO GAMMAS 
1)0 .30 I = 1 »N 
SH( 1 > = SHt t 1-53300. 

30 SL t I 1 =SL t I 1-53300. 
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C TO CALCULATE the VERTICAL GRADIENT FOR PROBE SEPARATION OF 
C 1.524 METERS 

DO <?0 1=1, N 

or 1 VC ! {)•=( Sl!tl-SM(t ) > / 1 • 5? 4 
THE DATA REDUCTION IS NOW COMPLETE 

tn? ST*T!ON=f f THE 4 NO MAE Y IN THE TOTAL FIELD * SH{ It GAMMAS 
The VERTICAL GRADIENT OF THIS = VO! II G/METER 
THE STATION IS 0(11 METERS FROM STATION 
THE STATION ELEVATION IS E ( I 1 METERS 

the t i m f dp nrnjPATtON is tu> hours 

the OUTPUT CAN COVSrST OF PRINTED DATA, PUNCHED CARDS, OR 
PLOTTED PROFILES. 

V< APPROPRIATE FORMAT SHOULD HE CHOSEN PY THE INTERPRETER. 

"A CMl MIT 
END 


SUBROUTINE YLNTPL CX,Y.XE.Y£,L> 

dimension xm,Ytn ,xei n ,yei n 

D8 1 = 1 
DO J=? 

100 [ F (XF(ll-XUm05,101,IO4 

id Yen )=vi j» 

OP TO 1 0 Ci 

104 J=J+1 

GO TO 100 

105 on SIS K=l,L 

IE < X! J-K) >110,114,112 

no writer, nn x,j,i 

111 FORMAT! ’ ’ILLEGAL VALUE ’tl?,' BEFORE ',1?,' p 0R *,!?» 

on to u s 

11? IF (Xfc m .fcO.XI J-K) » GU TO 113 

YEm*YI J-K)M <YU»-Y(J-K> )/(X( J>-XC J-K) ) ) * < XE (I) -X [ J-K 1 ) 
GO TO 106 

in YC( ! )'V(J-K 1 
GO TO 106 

114 IF { (J-K)-ll 110,112,115 

1 15 CONTI NUF 

106 1-1=1 

10T IF (I .r.T.L) GO TO 109 
GO TO 
1 ON IM t URN 
END 


ON 

CO 
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APPENDIX B. la 


Survey Information for Magnetometer Traverses 
in the Valley of Ten Thousand Smokes 




Date 


Beari ng 


Survey 

T raverse 

Type* 

July '69 

Duration 

(degrees E of N) 

Va 1 ley 

Party* 

A-A' 

TG 

14 

13:40-18:46 

36 

Lowe r 

TSB 

B-B ‘ (1) 

TG 

13 

10:15-13:25 

22 

Broken 
& Knife 

TSB 

B-B 1 (2) 

TG 

14 

9:55-12:30 

22 

Knife 

TSB 

B-B' ' 

T 

7 

13:48-15:12 

57 

B roken 

TSB 

B-b 

V 

10 

18:17-19:05 

57 

Broken 

F 

B 1 ' -b‘ 1 

T 

\ 7 

12:50-13:27 

116 

Broken 

TSB 

b'-b* 1 

V 

10 

17:50-18:01 

1 16 

Broken 

F 

b-b ' 

V 

10 

18:02-18:16 

114 

Broken 

F 

B-B' * ' 

T 

7 

9 : 35 - 11:26 

79 

Broken 

TSB 

B-B 1 ' ' 

V 

10 

16:27-17:41 

79 

Broken 

F 

b« * « -b 1 ' 

T 

7 

13:31-13:45 

188 

Broken 

TSB 

C-C* (3) 

V 

5 

14:12-17:10 

93 

Upper 

TF 

C-C (4) 

V 

6 

11:20-15:25 

93 

Upper 

TF 

D-D 1 

TG 

10 

13:00-17:50 

94 

Upper 

TSB 

E-E 1 

TG 

8 

10:48-14:44 

182 

Novarupta 

TSB 

F-F 1 

V 

14 

afternoon 

128 

Novarupta 

F 

F 1 -F' ' 

V 

14 

afternoon 

84 

Novarupta 

F 

F' 1 - F 

V 

14 

afternoon 

175 

Novarupta . 

F ■ 

G'-G 

V 

23 

17:00-17:30 

76 

Novarupta 

SB 

G l -G* 1 

V 

23 

17:30-17:36 

151 

Novarupta 

SB 

G ' 1 1 ~ 9 

V 

23 

17:40-18:10 

65 

Novarupta 

SB 

H-H * 

V 

23 

16:30-16:48 

131 

Tri dent 

SB 

I-! * 

V 

23 

15:27-16:09 

127 

Trident 

SB 

1 « _ 1 « t 

V 

23 

16:09-16:14 

‘ 37 

Tri dent 

SB 

1 " '-1 " 

V 

23 

16:14-16:19 

127 

Tri dent 

SB 

l-J 

V 

23 

14:44-15:27 

66 

T r i dent 

SB 

J‘-J 

V 

19 

18:56-19:15 

195 

Fenner 
Ri dge 

TB 

J * ' -J * 

V 

19 

19:15-19:20 

154 

Fenner 
Ri dge 

TB 

J-K M ‘ 

V 

19 

18:49-18:56 

49 

Fenne r 
Ri dge 

TB 

K‘ 1 '-K* ‘ 

V 

19 

17:27-18:49 

57 

Kn i fe 

TB 

K*-K' 1 

V 

19 

15:29-17:27 

152 

Kn i fe 

TB 

K-K' 

V 

19 

12:21-15:29 

112 

Knife 

TB 

L-L ' 

V 

17 

af te rnoon 

191 

Novarupta '5) 

F 

Fumarol e 
Grid 

T 

12 

10:22-13:48 

22 

Broken 

TSL 

Fuma role 
Grid 

T 

16 

afternoon 

22 

Broken 

TGM 
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> From 
i From 
i From 
i From 
1 From 

V = 

T = 

G — 

Baked Mtn. to Knife Creek 
Knife Creek to Mt. Griggs 
Baked Mtn. to mid-valley 
mid-valley to the Buttress Range 
Stumbling Mtn. to Broken Mtn. 

vertical field 
total field 

vertical gradient of total 

** T = Trible 
S = Stone 
B = Bingham 
F - Forbes 
L = Lofgren 
G = Gedney 
M = Matheson 




Table B. 1 b 


Survey Information 
in the VI cini ty of 

for Spot Readings of the Vertical Field 
Novarupta Volcano (measured by Forbes) 

at ion 

Date 

July '69 


Time 

Location 

1 

23 


16:52 

rim of Novarupta Crater 

2 

23 


16:54 

crater rim 

3 

23 


16:56 

crater rim 

4 

13 


? 

summit of Novarupta 

5 

17 


afternoon 

crater rim 

6 

17 


afternoon 

crater rim 

7 

13 


7 

west slope of Stumbling Mtn. 

8 

17 


afternoon 

crater rim 

9 

? 


? 

crater rim 

10 

17 


afternoon 

summit of Stumbling Mtn. 

H 

13 


? 

crater rim 

12 

14 


afternoon 

west slope of Broken Mtn. 

13 

17 


afternoon 

west slope of Broken Mtn. 

14 

17 


afternoon 

south slope of Greasy Pass 

15 

13 


? 

Greasy Pass 

15 

17 


afternoon 

Greasy Pass 



APPENDIX 0.2 

Magnetoneter Survey Data 


STATION 
N'<J M P r ? 


1 

2 

3 

4 

5 
ft 
7 
a 
o 

10 

n 

n 

14 

15 
15 
17 
10 
1° 
20 
21 
2? 
2? 
74 
2 4 
,?6 
27 
25 
-*o 

30 

31 

32 

33 
74 
35 
?6 
37 
33 
39 


MAGNETOMETER TRAVERSE A -A' 


magnetometer TRAVERSE A -A' (CONT.l 


STATION 
nl STANCE 
t**rTERS> 


STATION 
ELEVATION 
t METERS) 


TOTAL 
anomaly 
< G A U M AS 1 


0.0 

55. 5 

1 05 . ° 
144.0 

2 20.1 

33S . o 

3 53.6 

4 I 4. 7 

570.5 
6 26.5 

532.6 

613.7 

694.7 

750.5 

006.6 
0 62. 6 
9 1 8 . 7 

974.7 
1 0 3 4 . 4 
1 056.9 

1143.0 
1 196.0 

1355.1 

1311.2 

1767. 3 

1423.4 

1479.5 

1674.6 

1 591.7 
1647.3 

1703. 8 

1 759.9 

1015.7 

1071.7 

1927.5 

1 943.5 
20 10.3 
2095.0 

2150.8 


407.5 
3° 4. 3 
3HC.6 

369.9 

7 6 9.9 

?72 . 4 
371 .4 
301 .7 

774.9 
37 7.9 
377.8 

377.8 
7 0 0.8 
3 06.9 
3 06.8 
3 83.7 

332.7 
7 80.5 
7 0.0.5 

330. 4 

700.3 
3 83.3 
7 03.2 

303.1 
3 0 3.0 

782.9 

332.9 

704.7 

785.7 

335.7 
382.6 

362.5 
3 7 6.3 

379.3 

305.3 

382.2 

308.3 

792.1 

388.1 


-106. 
-13 3. 
-2 39. 
400. 
17?. 
781 . 
4?7. 
299. 
405. 
419. 
375. 
792. 
474. 
-17?. 
-40 5 . 
32. 
191 . 
-47. 
6 66 . 
301. 
343. 
148. 
157. 
359. 
183. 
194. 
127. 
24 4. 

213. 
-109. 
74. 
236, 
56. 
2 32. 
153. 
-? 7 . 
-9 2. 
177. 
13 3. 


VERT HER 

STATION 

STATION 

station 

op total 

NUMBER 

DISTANCE 

ELEVATION 

{ G/MET ER I 


IMETERS) 

{mptERS) 

-4.6 

40 

7206. 1 

■378.9 

-17.4 

41 

2262.1 

381.9 

-??. 7 

4? 

2317.9 

787.9 

5.4 

43 

2374.0 

387.8 

62.7 

44 

2430.0 

?04.7 

68.4 

45 

24 86. I 

387.0 

1.6 

4 6 

2542.1 

389.2 

-3.1 

47 

7598.2 

342.2 

8.5 

48 

2654.2 

390. T 

22.5 

49 

2710.0 

7^6 . 7 

34.0 

50 

2766. 1 

398.2 

-49 . 2 

51 

7827.2 

399.7 

7.0 

52 

2077.9 

405.0 

-45.4 

53 

2934.0 

408.8 

-84.4 

54 

2983.7 

470. 9 

-70.0 

55 

3044. 7 

423.0 


-9.2 

40.4 
172.0 

0.0 
-2.3 
16.9 
-5 6.1 
0.0 
-19.2 
-19.9 
-18.4 

20.5 
-13.8 

-117.9 

64.3 

40.5 
29.0 

- 10.0 

2.3 

-47.3 

-50.8 

-7.7 

-8.4 


TOTAL 

ANOMALY 

(GAMMAS1 


VERT OEft 

OP TOTAL 

f c./metlri 


-109. 

-0.3 

219. 

-17.6 

7 3 5. 

-9.7 

4 7. 

-16.8 

19 7. 

-10.7 

14 5. - 

-17.5 

118. 

-3.0 

701 . 

- 7 6 . P. 

-37. 

-35.7 

-54 1 . 

-43.7 

30. 

11.4 

-215. 

0.0 

-218. 

- 19. ft 

-647. 

-13.4 

847. 

-785. 1 

-420. 

36.2 



MAGNETOMETER TRAVERSE B 


-R* 


STATIST 

STATION 

STATION 

NU V *ER 

01 STANCE 

ELEVATION 


(METERS) 

( METERS) 

l 

0.0 

740.3 

7 

77.6 

746.1 

? 

65.4 

741.0 

4 

>33.1 

7 3 7.0 

5 

110.2 

776.4 

& 

1 34. 7 

715.5 

7 

160.3 

704,.? 

a 

192.3 

696.0 

a 

206.5 

0 30.3 

to 

2 63. 7 

66 5 . FI 

11 

316.5 

659.9 

17 

375.7 

643.2 

n 

631.3 

653.9 

i* 

>.3 7.4 

654.0 

1 5 

54?. 4 

651.0 

16 

540.5 

648.0 

17 

655.6 

64 6.6 

■ IS 

7U.6 

645.2 

I'J 

367 . 6 

641 .7 

70 

•53 3.6 

63 9.1 

?1 

8.7 >. 6 

64 0.9 

22 

035.7 

640.1 

2 3 

. .8 

639.3 

24 

1047.9 

641.4 

25 

1103. 3 

644.3 

26 

H5S. 0 

655.7 

77 

1714.4 

66 3.3 

28 

12TC.5 

664.7 

29 

1326.4 

66 9.4 

’0 

1337.4 . 

666.9 

01 

1410.0 

677.1 

02 

1431.0 

657.6 

13 

1457.5 

647.5 

34 

1411.1 

677.3 

35 

1508.1 

619.7 

36 

1 534.0 

61 1 .3 

17 

1557.0 

604.0 

33 

t 6<J ). 0 

603.1 

30 

1646. 1 

592.6 


TOTAL 

VERT OER 

ANOMALY 

OF total 

(GAMMAS! 

(G/METER) 

72. 

-6.3 . 

-99 . 

- 14.9 

-72 9. 

-53.0 

- 249 . 

-61.6 

-539, 

-80. 1 

-375. 

-21.7 

-715. 

-10.9 

03. 

-3 0.8 

317 . 

- 16.7 

761 . 

-21.5 

463 . 

-15.4 

439 . 

-32.1 

3 0 1. 

-52.7 

170. 

-64 . 5 

177. 

-7.9 

190. 

-33.3 

433. 

-53.6 

573. 

-27.3 

7,4 2. 

13.7 

96. 

-62.9 

- 210 . 

-57.0 

-0. 

-45.7 

• - 304 . 

-17.9 

-255. 

-17.2 

-35 0. 

-34 . 2 

-379. 

-44.3 

-52 l . 

0.0 

-512. 

— 10.0 

- 461 . 

-2 7.9 

- 355 . 

-29.5 

- 559 . 

-335.7 

-40? . 

25.0 

-25 7. 

-31.2 

50. 

4.7 

362. 

-6.4 

477. 

-41.7 

607. 

-3.2 

*0 6 . 

“74.0 

746. 

-35.3 


MAGNETOMETER TRAVERSE B 


tCONT.) 


-S' 


STATION 

STATION 

ST AT 1 ON 

TOTAL 

VERT OER 

NUMBER 

OT STANCE 

ELEVATION 

ANOMALY 

OF TOTAL 


(METERS) 

(.METERS) 

(GAMMAS) 

(G/V.ET5R) 

40 

1700.9 

587.1 

222. 

-34.9 

41 

1725.6 

582.6 

1 54 . 

-4 1 . 2 

4? 

1756.7 

532.3 

215. 

-32.6 

4 3 

1812.3 

592.7 

-6. 

-55.1 

44 

1063.7 

579.1 

165. 

1.6 

45 

1924. 0 

577.7 

24. 

— 41.8 

46 

1 980.8 

573.8 

-4. 

-55 . 1 

47 

2036.5 

567.1 

309. 

-2 6.3 

43 

2064.5 

646.6 

295. 

31.1 

49 

2170.6 

55 l .2 

345. 

53.9 

50 

2176.1 

559.7 

417. 

14.7 

51 

2232.0 

555.2 

495. 

-3.1 

5? 

2?H7. 9 

560.2 

44 3. 

1.6 

53 

7344.0 

560.0 

46 3 . 

17.1 

54 

2400. 1 

660.0 

4? 2 ♦ 

20.9 

55 

2456.1 

561.6 

646. 

10.9 

56 

7512.2 

563.2 

675. 

10.2 

57 

2568.3 

563.3 

79 2. 

26.7 

60 

2673, 9 

871 .0 

722. 

18.0 

59 

2630.0 

549.5 

68 7. 

13.0 

60 • 

2734.1 

560.4 

51 5 . 

3.9 

61 

2 7)2,2 

569.6 

379. 

-10.8 

62 

- 2848.2 

463.1 

350. 

-3.1 

63 

2904.3 

570.6 

355. 

-70.1 

44 

2960.4 

571.3 

439 . 

34. { 

65 

3016.5 

571 .4 

33 3. 

4.6 

66 

3072.6 

669.9 

365. 

6.2 

67 

3123.6 

671.4 

37 7. 

1.5 

68 

3104. 7 

571.5 

31 1 . 

-0.8 

69 

3240.0 

560.5 

333. 

9.3 

70 

3296.8 

570.0 

219. 

12.3 

71 

3352.9 

670,1 

74. 

-23.7 

72 

3409.0 

570.1 

67. 

-9.9 

73 

3465.1 

568.6 

98. 

-8.4 

74 

3521.2 

670.2 

178. 

13.0 

75 

3577.2 

567.2 

190. 

-3.0 

76 

3633.3 

565.7 

ICO. . 

-2.3 

7 r 

36H9. 4 

V. t . 1 

7 . 

It .4 

78 

3745.5 

567.3 

-105. 

8.3 


N> 






MAGNETOMETER TRAVERSE 

3 -O’ 

ICONT. ) 


MAGNETOMETER 

TRAVERSE D»* 

-b* • » 

STATION 

STATION 

STATION 

TOTAL 

VERT DER 

STATION 

STATION 

STATION 

TOTAL 

YJmreR 

DISTANCE 

ELEVATION 

ANOMALY 

OF TOTAL : 

NUMBER 

distance 

' elevation 

ANOMALY 


imstfrsi 

(METERS) 

( 0 AVMASl 

(G/MRTERI 


(METERS ) 

( METERS) 

(GAMYAS1 

7° 

SADI .6 

597.4 

-420. 

-9.7 

1 

3.0 

733.1 

474 . 

*0 

>K<7. 7 

*67.4 

-753. 

1.5 

2 

9.4 

733.0 

491 . 

81 

3913.3 

567.5 

-91. 

-1.5 

3 

1.2 / 

732.9 

4 64 . 

92 

3964.9 

567.5 

-138. 

-1.5 

4 

L.R 7 

732.8 

462 . 

■ 33 

4(1?6. 0 

567.9 

-125. 

-1.5 

5 

2.4 

737.7 . 

452. 

84 

4981. 7 

573 .7 

343. 

-4.6 

6 

3.0 

732.6 

46?. 

36 

413?,. I 

SH7. 5 

397. 

-49.4 

7 

3.6 

732.5 

36? . 

8 6 • 

4190.9 

599.7 

39. 

-23.6 , 

8 

4.2 

732.4 

4 34. 

87 

4244.8 

615.0 

373. 

-29.3 

9 

4.8 

732.3 

46.3 • 






10 

5.4 

732.2 

423. 






J 1 

6.0 

73 7.1 

389. 






12 

6.6 

7 52.0 

34 0 . 






13 

7. 7 

731.o 

305. 






14 

7.8 

731 .B 

304. 






15 

8.4 

731.7 

246. 



MAGNETOMETER 

TRAVERSE b 

- b* 

STATION 

station 

station 

VERTICAL 

NUMBER 

DISTANCE 

ELEVATION 

ANOMALY 


1 METERS ) 

(METERS) 

(GAMMAS) 

l 

0.0 

6A0.7 

93. 

2 

2 A. 0 

680. 8 

-4 7. 

3 

56.1 

681.0 

33. 

4 

84.1 

682.5 

-167. 

5 

112. 1 

634.0 

-152. 

6 

140.2 

635.5 

-257. 

7 

169.2 

637.0 

-197. 

3 

196.2 

688.5 

-427. 

9 

224.2 

691.5 

-77. 

10 

257.3 

694. 5 

13. 

1 V 

2R0.3 

697.5 

73. 

12 

303.3 

700.7 

63. 


\JS» 




MAGNETOMETER TRAVERSE 

B -S’ 1 


STATION 

STATION 

STATION 

vertical 

TOTAL 

NUM RER 

DISTANCE 

FV.FVAT ION 

AimOMAL v 

ANOMALY 


1 METERS) 

(METERS) 

(GAMMAS) 

(GAMMAS) 

1 

0.0 

765.2 

-32,0 

-95.0 

2 

5. I 

764 .0 

'> H . 0 

****** 

2 

32. 1 

75 3.6 

-57.0 

****** 

4 

95. 1 

754.6 

* «*#** 

- 5 n • 0 

S 


75 9.2 

-97.0 

& A ft <r ijf ^ 

6 

35.1 

7 4 7.8 

-47.0 

£ 9 * * <r * 

7 

110.7 

747.0 

****** 

-167.0 

R 

111.7 

742.3 

-57.0 

****** 

o 

T . 1 

796.9 

- 102.0 

****** 

10 

164.6 

731.5 

-287.0 

-276.0 

11 

17 3.5 

729.9 

- 102.0 

****** 

12 

197.9 

728.3 

-107.0 

****** 

) i 

70. ' 

77‘i . 2 

- 1 1 7 . i) 

- 3 w. 0 

14 

242.9 

703.3 

. -267.0 

****** 

IS 

265.1 

691 .4 

213.0 

267.0 

16 

7 ,/> ? . 6 

684.fi 

239 . 0 

****** 

17 

320.1 

■680.2 

198.0 

206.0 

19 • 

?4B . 0 

677.1 

223.0 

****** 

1 ° 

778 . 9 

674.1 

303.0 

265.0 

?<> 

4 9 1.9 

6 7 2.4 

V 3.0 


21 

*31.9 

670.8 

153.0 

140.0 

2 ? 

440.9 

670.5 

203.0 

* -<* 9>(* fr 

7 S 

4 8 !(.0 

6 7 0. ) 

sm.il 

t 6 *> * 0 

24 

516.0 

671.3 

243.0 

$4 £* * 

25 

544.0 

672.2 

173.0 

136.0 

2t> 

5 72 . 1 

67 0.7 

3 4 3.0 

* 4 <l <* * ? 

27 

600. 0 

669.2 

273.0 

90.0 

7? 

628.1 

670.3 

73.0 

* * $ * 

2 « 

64r- . t 

6 n . 4 

283.0 

too.o 

TO 

634.1 

673.0 

2 5 3.0 

****** 

?1 

712.1 

674.5 

243.0 

107.0 

82 

740.2 

674.5 

133.0 

****** 

38 

768.2 

674.5 

133.0 

-32.0 

34 

796.2 

676.4 

-27.0 

****** 

29 

4 74.2 

678.2 

13.0 

-115.0 

8 4 

852.3 

670.5 

13.0 

****** 

97 

990.2 

680.7 

98.0 

-24.0 

3 3 1 

9 16. 3 

680.8 

****** 

-177.0 

30 

°97.4 

682.3 

****** 

-254.0 


MAGNET OMETFR TRAVERS F B -8” 5CQNT . ) 


STATION 

STATION 

STATION 

VERTICAL 

TOTAL 

N1JMRFR 

0! STANCE 

ELEVATION 

ANOMALY 

ANOMALY 


(METERS) 

(METERS) 

(GAMMAS) 

(GAMMAS) 

40 

1049.0 

690.0 

****** 

-332.0 

41 

1102. 7 

702.3 

****** 

-594.0 

42 

1 156.7 

, 717.5 

****** 

-304.0 

43 

1210.6 

/ 733.1 

****** 

230.0 


MAGNETOMETER TRAVERSE B*' -b‘ ’ 


STATION 

NUMBPR 


I 

? 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 
18 
14 


station 

STATION 

distance 

ELEVATION 

(MCIFRSl 

(ME TLPS) 

0.0 

733.1 

27.7 

728.6 

ft 0.9 

710.9 

136.5 

703.0 

192.2 

696.1 

249.1 

700.7 

27 A. 1 

699.8 

304.1 

698.9 

332.1 

698.6 

340.2 

698.4 

380.2 

700.0 

416.2 

701 .5 

444.2 

701.4 

472.3 

701.4 


vertical 

TOTAL 

anomaly 

anomaly 

(GAMMAS) 

(GAMMAS) 

****** 

230.0 

****** 

290.0 

****** 

-380.0 

****** 

-173.0 

****** 

-104.0 

63.0 

-41.0 

73.0 

****** 

-47.0 

-152.0 

-47.0 

****** 

-62.0 

-184.0 

-9 7.0 

****** 

-107.0 

-229.0 

-117.0 

****** 

-17.0 

- 15L.0 




MAGNETOMETER TRAVERSE B -B'*' 


STATION 

station 

STATION 

VERTICAL 

TOTAL 

NU M ? E R 

01 STANCE 

ELEVATION 

ANOMALY 

anomaly 


(MATERS) 

(METERS) 

( GAMMAS ) 

t GAMMAS) 

1 

0.0 

766.6 

-32.0 

-77.0 

2 

26. 8 

76 3.0 

-«2. 0 

****** 

7 

53.6 

759.4 

-87. 0 

****** 

4 . 

56.0 

763.6 


-103,0 

$ 

84. 6 

744.7 

-77.0 

* ft ft ♦ ft * 

A 

107.2 

752.1 

-107.0 

ft ft ft ft ft # 

T 

111-1 

752.7 

ft * ft ft ft 

-1 l?.0 

8 

136.0 

74 fi , 5 

-107.0 

ft << ft P. ft ft 

o 

160.9 

744.9 

3.0 

ftftftft ft ft 

10 

166.5 

744.1 

-2.0 

-5T.0 

11 

193.2 

735.6 

-297.0 


i: 

219.0 

727,6 

-347.0 

ft ft ft ft ft ft 

n 

220. 3 

726.2 

* ***** 

-347.0 

14 

242.8 

719.7 

-4\7.0 

« o ft ♦ ft* 

IS 

246.6 ' 

71° .5 

-477.0 

****** 

16 

273.3 

7 09,9 

-197.0 

6.0 

17 

300.9 . 

704 .9 

243.0 

****** 

i* • 

329. 5 

699.7 

198.0 

22% mO 


3 36.4 

697.5 

183.0 

* ftC* ftft 

?o 

383.7 

659.9 

****** 

1.6*. 0 

21 

336 . R 

694.3 

143.0 

*£ A * ft £ £ 

22 

414.2 

693.0 

223.0 

ft* ft A <f * 

23 

637.3 

682.2 

****** 

?65 .0 

24 

441.5 

69C.9 

293.0 


25 

469.8 

68 8. 7 

523.0 

* ft ft ft ft <i 

26 

405.2 

636.7 

413.0 

3 3 3-0 

27 

530.5 

687.0 

198.0 

* $ $ 

?P 

54 t. 3 

697.6 


313.0 

2° 

6^5.7 

6 87.3 

373.0 

* « **« * 

. 30 

60 1.2 

687.6 

263.0 

****** 

31 

607.4 ■ 

ese.r 

«««««* 

349.0 

32 

636.7 

6B7.9 

383.0 

****** 

33 

663.5 

6* 8- 1 

323.0 

326.0 

34 

686.6 

687.7 

433.0 

****** 

35 

710.5 

637.2 

273.0 

****** 

36 

717.5 

6«5.0 

* * * * <i * 

9 3.0 

37 

732.9 

686.8 

323.0 

****** 

39 

747.5 

686.5 

283.0 

214.0 

30 

7 79.0 

680.0 

303.0 

* *<■ * ** 



MAGNETOMETER TRAVERSE 

ft -0**’ 

(CDNT.l 

STATION 

STATION 

STATION 

VERTICAL 

total 

NUMBER 

DISTANCE 

ELEVATION 

ANOMALY 

ANOMALY 


( METERS! 

(METERS) 

(GAMMAS) 

( GAMMAS) 

40 

807.5 

639.4 

****** 

3TI .0 

41 

809.2 

689.4 

343.0 

****** 

4? 

8 31.6 

689.4 

****** 

257.0 

43 

838.7 

690.9 

363.0 

****** 

4 4 

869.0 

' 64? . 4 

723.0 

****** 

45 

88 7.5 

693.9 

****** 

183.0 

46 

999.4 

693.9 

463.0 

****** 

4 7 

929.3 

695.3 

363.0 

* « * * * * 

4? 

943.6 

696.0 

323.0 

150.0 

49 

974.6 

697.0 

263,0 

****** 

50 

999.7 

695.4 

* ***** 

64,0 

51 

1005.7 

698.0 

223,0 

* ** * * * 

52 

1036.9 

699.0 

223.0 

****** 

53 

1056. 7 

697.5 

253.0 

109.0 

54 

1 093.7 

700.5 

153.0 

****** 

55 

1111.7 

701 .4 

3.0 

-151.0 

57, 

1139.5 

705.6 

-92.0 

****** 

57 

1167.2 

709.8 

-7 2.0 

-274.0 

58 

1 1 84.9 

7(4.2 

-237.0 

****** 

59 

1202.5 

718.6 

-397.0 

****** 

60 

1220.0 

723.0 

-527.0 

* * * * * * 

61 

1230.3 

725.6 

****** 

-550.0 

62 

1 277.6 

727.4 

-437.0 

ft * ft ft 

63 

1255.3 

731.8 

-352.0 

* ft ft ft ft ft 

64 

1269.9 

736.5 

-347.0 

ft tr t 0 ft ft 

65 

1284.7 

739.2 

-347.0 

-385.0 


Vn 




MAGNETOMETER 

TRAVERSE C 

- C* 


magnetometer 

TRAVERSE C 

- C f (CONT 

STATION 

3 t a r r r.N 

STATION 

VERTICAL 

STATION 

STATION 

station 

VERTICAL ' 

N<J*8ER 

01 STANCE 

ELEVATION 

ANOMALY 

number 

distance 

ELEVATION 

ANOMALY 


IMSTEPS1 

{METERS! 

{GAMMAS! 

{METERS! 

(KFTFRS J 

{GAMMAS ! 

1 

0.0 

605.9 

123. 

40 

770.5 

532.5 

298. 

?. 

30.5 

596.2 

173. 

41 

601.0 

532.5 

233. 

3 

38.4 

593.3 

43. 

42 

316.9 

532.2 

-47. 

u 

M.O 

691.9 

-197. 

4? 

531.8 

531.9 

78. 

5 

91.4 

5 39.5 

-37. 

44 

847.3 

531.3 

373. 


107.0 

587.4 

>207. 

45 

671.1 

531.3 

34 3. 

7 

122.2 

535.2 

313. 

46 

8 ? 3 . 1 

5 30.4 

773. 

P 

1 36.2 

583. 1 

413. 

47 

920.6 

529. 1 

53. 

9 

133.0 

5.40.6 

24 3. 

48 

923.9 

520.1 

173. 

i n 

193.5 

576.1 

333. 

49 

954.6 

530. 0 

333. 

U 

213.7 

569. 7 

73. 

60 

964. 1 

530.4 

228. 

12 

2 ?rt .9 

546.6 

33. 

51 

9 8 4.7 

530.4 

23 3. 

1 1 

244. 1 

563. 6 

238. 

57 

1016.5 

530.4 

133. 

!'• 

272.2 

559.4 

233. 

63 

1022.9 

530.4 

173. 

is 

2 74.6 

657.4 

?n. 

54 

1047.6 

4 3 0.4 

273. 

16 

305.1 

• 552-9 

233. 

55 

1076.4 

530.4 

36 3. 

1? 

3 ?4 . 3 

6 4'?. q 

263. 

56 

1080.8 

530.4 

4?3. 

13 

336.8 

544.0 

423. 

57 

1100.8 

530.4 

193 . 

19 

352.3 

546.6 

123. 

58 

11 16. 8 

530.4 

273. 

20 

363.6 

544. ! 

243. 

59 

1140. 6 

529.4 

103. 

21 

3 6 7.9 

543.8 

233. 

40 

1 1 49.0 

5 ? a . 2 

773. 

22 

39 3.7 

542.2 

183. 

M 

1171.7 

528.7 

273. 

23 

670. 2 

540.4 

233. 

62 

1200.3 

5 2 1.6 

30 3. 

74 

4 51.7 

5 39.2 

237. 

63 

1270.4 

5 2 7.3 

?4 3 . 

25 

4 59.6 

5 3 8.9 

317. 

64 

1230.8 

627.0 

2 4 H . 

2o 

490.4 

536.0 

31 3. 

66 

1761.6 

526.4 

423. 

27 

504. 1 

537.7 

333. 

66 

1273.3 

526. 1 

5 83. 

26 

52 1 .? 

537.4 

31 3. 

6 7 

1 2 H9. 9 

525.6 

483. 

2 ° 

5 42.3 

936.4 

223. 

66 

1320.7 

5?4. 3 

37 3. 

30 

564.8 

536.1 

103. 

69 

1336.2 

524.0 

673. 

31 

58?. 1 

535.4 

3. 

70 

. 1346.9 

523.0 

483. 

3? 

614.2 

535.5 

473. 

71 

1377. 7 

520.3 

47 3. 

13 

616.0 

535.2 

62 3. 

77 

1 392.6 

575.5 

733. 

34 

64-3.0 

534.9 

413. 

73 

1405 . 1 

577.0 

623. 

35 

676.0 

534.3 

-67. 

74 

1420.7 

577.0 

3 3. 

3 6 

707. 1 

634.0 

32 3. 

76 

1451.5 

52 7.0 

223. 

37 

709.6 

571.4 

383 . 

76 

14!?2.2 

527.0 

173. 

3? 

740.1 

53 3. 1 

323. 

77 

1485.3 

527.6 

173. 


7i.l .7 

53 2. 8 

273. 

78 

1514. J 

57 /. 3 

-7. 




magnetometer 

TRAVERSE C 

- C» (CONT 

STATION 

STATION 

STATION 

VFRTICAL 

NUMBER 

distance 

ELEVATION ■ 

ANOMALY 


(METERS! 

(METERS! 

(GAMMAS) 

79 

1544.1 

527.0 

323. 

*>n 

1562.7 

527.0 

257. 

51 

1575.2 

5 7 7.3 

-87. 

8 2 

1606.0 

527.9 

-437. 

93 

1621.5 

528.2 

123. 

?4 

1634.0 

52 8.8 

-477. 

95 

t 6 6 5 . 1 

5?9. 7 

28. 

96 

1695.9 

570. 7 

173. 

87 

I 726.4 

530.7 

153. 

83 

1735.5 

530.7 

163. 

89 

1763.5 

530.7 

198. 

90 

1796.8 

530.7 

203. 

91 

179.3.9 

530.7 

203. 

9? 

1329.7 

530.7 

303. 

93 

1356. 2 

530.7 

247. 

94 

1360.5 

530.7 

233. 

95 

1791.6 

530.7 

323. 

96 

1915.7 

530. 7 

223. 

97 

1922.4 

530.7 

238. 

92 

1963.2 

53C.7 

123. 

99 

1970.2 

530.7 

703. 

100 

1932.7 

530. 7 

373. 

101 

2013.8 

530. 7 

443. 

107 

2C’9.l 

530.7 

383. 

193 

2047.7 

5 30.7 

463. 

104 

2073.7 

530.7 

443. 

105 

7109.5 

530.7 

213. 

106 

2119.6 

571.0 

20 3. 

107 

2177.9 

571.3 

1 78. 

' 109 

21 63.7 

531.9 

28. 

109 

2196.4 

572.6 

53. 

no 

2232.4 

533.4 

53. 

ui 

2261.3 

534.0 

153. 

112 

2297.4 

534.9 

3. 

113 

232^.7 

575.6 

3. 

114 

2354.0 

535.8 

-17. 

115 

2334.8 

535. 8 

-62. 

116 

2415.7 

5 76.1 

283. 

117 

2456. 1 

536. 8 

243. 


MAGNETOMETER TRAVERSE C 


C* (CONT.) 


STATION 

STATION 

STATTON 

NUMBER 

distance 

ELEVATION 


( METERS1 

(METERS) 

158 

2486.6 

537.7 

119 

2517.0 

538.6 

! 20 

2547.5 

579.5 

121 

2578.0 

540.4 

122 

2608.5 

54 1.3 

123 

2639.0 

547.2 

174 

2669.4 

547.2 

125 

2700.5 

545.6 

126 

2731.3 

54 R. 0 

177 

2745.3 

5 4 ». . 9 

128 

2762.1 

850.5 

1?9 

2793.2 

552.6 

130 

2807.6 

557.9 

131 

7824.3 

55 7. 8 

172 

2858.1 

558. 1 

1’3 

7836.2 

560.2 

174 

2898.9 

561.1 

136 

2914.2 

567.6 

136 

2943.8 

56 6.6 

137 

7853. 7 

563.5 

1 3R 

2973.7 

870.0 

139 

3004.4 

877.9 

140 

3032.2 

577.9 

1.4 1 

3035.2 

578. 2 

1 42 

3067 . 6 

5 f! 3 . 1 

143 

3093.1 

8 8 8.6 

144 

3108.4 

. 5 39.3 

145 

3126.6 

592.5 

146 

3144.9 

595.0 

147 

3157.4 

596.8 

148 

3188.2 

600.8 

149 

3205.0 

602.9 

150 

3219.0 

6Q4.4 

151 

3250.7 

60P.7 

15 2 

3252.5 

508. 1 

163 

3279.4 

617.6 

154 

3281 .2 

614. 8 

155 

3312.0 

61 8.7 

156 

3727.2 

621.2 


VFRTTCAl 
ANOMALY 
f GAMMAS 1 

1 1 7. 
133. 
263. 
773. 
193. 

7 2 3. 
363. 

2 7 7 . 
273. 
263. 
773. 
323. 
243. 
273. 
723. 
103. 

7. 
-77. 
-147. 
-237. 
-4 77. 
-117. ' 
-312. 
-367. 

— 4 7. 
-142. 
-7. 
-47. 
-217. 
-777, 
-52. 
-227. 
73. 
-32. 
-67. 
-72. 
-77. 
-117. 
-6 7 . 


•VJ 



M4GNET0KETER TRAVES.SE C 


C» (CONT.l 



STATION 

STATION 

.NUMBER 

DISTANCE 
(ME TE«S ) 

157 

334?. 4 

1 58 

3372.9 

159 

3397.0 

150 

3433.1 

161 

3433.0 

162 

3465.8 

153 

3402.8 

164 

3433.3 

165 

34-J2.1 

166 

3524.7 


STATION 

VERTICAL 

ELEVATION 

ANOMALY 

< METERS! 

(GAMMAS ) 

623.6 

-137. 

62 8.5 

-297. 

634.6 

-447. 

636.1 

-447, 

643.7 

-447. 

660. 1 

-417. 

651.7 

-45Z. 

659.9 

-377. 

662.3 

-347. 

676.4 

-347. 


STATION 

NIJMRFR 


1 

? 

3 

4 

5 

6 

7 

8 
9 

10 
I l 
12 
n 

15 

1* 

17 

IB 

19 

20 
21 
22 
28 

24 

25 

24 

27 

28 
?? 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 


MAGNETOMETER TRAVERSE 0 


- 0 * 


station 

STATION 

TOTAL 

VERT DER 

DISTANCE 

ELEVATION 

anomaly 

OF TOTAL 

( mete R S ) 

< METERS) 

(GAMMAS) 

CC/METCRJ 

0.0 

62 5.5 

-32. 

-0.3 

56.6 • 

617.9 

313. 

■ 103.2 

t 10.3 

605.8 

334. 

16.8 

165.8 

597.6 

172. 

4.0 

221.1 

688.1 

364 . 

-42.3 

2 76.1 

576.9 

270 . 

-(9. 3 

331.6 

563.4 

44 3. 

35.3 

3“ 7. 6 

546.0 

470. 

1.8 

443.5 

661.7 

24 7. 

5 . 8 

499.6 

560.2 

302. 

-3.8 

5 54. 7 

549.6 

1 39. 

-13.5 

610.5 

555.6 

91 . 

-11.9 

665.6 

5 5 5.0 

-59. 

-75. 1 

722.6 

552.6 

113. 

1.6 

778.7 

582.3 

-74. 

-29.8 

8 34.8 

842.0 

-152. 

-29.7 

890.9 

551.1 

-79. 

— 4 3 * 6 

847.0 

549.6 

110. 

-4.0 

1 003.0 

552.4 

163. 

13.2 

1069.1 

851.3 

132, 

19.9 

1115.2 

551 .2 

272. 

-177.4 

1171.3 

550.6 

433. 

0.0 

1339.5 

54 5.1 

573. 

-32.9 

1367. 6 

546.0 

191. 

44.5 

1395.6 

54 7.0 

-99. 

19.3 

1451.7 

547.7 

214. 

60.4 

1 507.3 

546. 5 

-14 l. 

-35.1 

1535.6 

552.2 

-83, 

7. 1 

1 563 . 6 

552.1 

-371. 

-92.9 

1619.7 

851.9 

-9 0. 

-1.8 

1675.8 

551.6 

-6 6. 

13.4 

1731.9 

553.2 

-218. 

3.9 

1738.0 

5 5 1.7 

-225. 

-70.1 

1844, l 

550.3 

131 . 

36.4 

1900.2 

552 .6 

-342. 

-85.4 

1966.2 

553.3 

-29! . 

10.9 

2017.3 

554 . 5 

-11. 

12.6 

206R.4 

555.7 

207. 

17.5 

2124.5 

658.7 

39. 

-id. 9 


'•-J 

CO 



MAGNETOMETER' TRAVERSE 0 -D« 


( CONT . > 


MAGNETOMETER TRAVERSE 0 


-0 


(CONT. » 


STATION 

STATION 

station 

TOTAL 

VERT OER 

NO*nE« 

DISTANCE 

(METERS) 

ELEVATION 

(METERS) 

ANOMALY 

(GAMMAS) 

OP TOTAL 
(G/ METER) 

4? 

2130.6 

656.4 

355. 

21.5 

*) 

2236.7 

657.9 

-2. 

-4. 7 

42 

.2 2 72.3 

557.9 

-22 4. 

-39.0 

43 

2320.3 

558.0 

-12. 

l .6 

44 

2343 .9 

558.0 

-124. 

— 7.8 

48 

24C8 ,0 

665.9 

16. 

0.8 

46 

2461. 1 

64 f>.? 

-257. 

-119.6 

47 

2517.2 

556.2 

-113. 

13.3 

48 

2573.2 

557.1 

-4 8. 

5.5 

4Q 

24?9.2 

561.1 

-4. 

-18. t 

59 

2665.3 

562.1 

-18. 

12.6 

5 l 

2741.4 

662.4 

98. 

II . 1 

52 

2797.5 

562.4 

66. 

— ?h . 0 

44 

2847.6 

562.4 

262. 

-28.6 

54 

290° . 7 

552.7 

281. 

-17.5 

54 

7064. S 

568.2 

93. 

-7.9 

46 

3021.5 

571.9 

-46. 

-12.6 

57 

3077.5 

574.9 

-101 . 

-6 . 5 

59 

3133.5 

576.0 

4 6. 

38.7 

40 

3199.4 

562.9 

-706. 

9.4 

60 

32*5.4 

4 37.2 

-313. 

10.9 

*1 

3301.3 

?9 1 .8 

-376. . 

-38.0 

4? 

3367.0 

608.2 

-63. 

11 .8 

6 3 

3412.9 

602.5 

-11. 

37.3 

64 

3440.6 

604.0 

-166. 

-0.8 

*>5 

34t>tJ.9 

606. 1 

-zie. 

— 98.1 

66 


607.1 

-157. 

29.0 

6 7 

3524. 8 

610.1 

— 1 86* 

7 . 8 

63 

3552.6 

614.1 • 

-201. 

-6.2 

60 

35R0.4 

617.8 

-124. 

-16.4 

70 

1 3608.1 

67 2.1 

-II. 

2.4 

71 

3636.0 

62 5.4 

— 2 6. 

-35.3 

72 

3663.3 

631 .8 

* -75. 

. — ?5 • 9 

73 

3691.2 

6? 4. 6 

-61 . 

— 42 . 4 

74 

3 7 1 7.4 

64 4.7 

91 . 

19.7 

75 

3745.3 

647.1 

134. 

-41.1 

76 

3772.9 

652.0 

58. 

-20.5 

7 T 

3 799. 1 

66 7.1 

-97. 

I 8. P 

76 

3«27.t 

663.0 

—370. 

4.7 


STATION ST AT ! DN 

number distance 

( METERS) 

79 385*. 9 

90 3910. 4 

M 3964.2 

92 *018.9 

83 4069.4 


STATION 

total 

ELEVATION 

anomaly 

(METERS) 

( GAMMAS 

667.3 

-451. 

675.3 

-484. 

/ 681.1 

-471. 

TOO. 3 

-498. 

724.7 

785. 


VERT DER 
OP TOTAL 
(G/M&TEP) 

-* 0.4 

- 3.1 

-4.6 

802.6 

-34.0 




MAGNETONS TER TRAVERSE 

E -E * 



MAGNETOMETER TRAVERSE 

E -E ' 

(CONT. 1 

STATION 

STATION 

STATION 

TOTAL 

VERT OER 

station 

STATTON 

STATION 

TOTAL 

VERT DER 

number 

. 01 STANCE 

ELEVATION 

ANOMALY 

OP TOTAL 

NUMBER 

Of STANCE 

ELEVATTON 

ANOMALY 

OF TOTAL 


(WETFRSI 

1 METERS) 

tO.AMMASl 

IG/METER) 


(METERS) 

(METERS) 

(GAMMAS) 

IG/ METER) 

1 

0.0 

747.1 

-997. 

10.7 

40 

l 139.5 

■ 682.8 

-2867. ' 

501.2 

2 

23. 5 

731.9 

-1327. 

-26.9 

41 

1167.5 

6°3.7 

784. 

175.1 

3 

53.6 

724.5 

-1398. 

-20. L 

42 

1 195.5 

683.8 

' 22 0. 

-65.0 

4 

77.1 

7! 5.4 

-1283. 

-18.7 

43 

1 223. 5 

/ 686.4 

-776. 

-30.6 

5 

104.7 

710.2 

-1016. 

-0.8 

44 

1 251 . 5 

688.9 

-30 6. 

-171.2 

6 

132. 1 

704.1 

-896. 

8.4 

45 

1279.4 

691.0 

-745. 

-163. 1 

7 

160.1 

707.0 

-664. 

24.6 

46 

1707.4 

691 .5 

-64 0. 

-151.5 

8 

18 7.9 

69R.3 

—929 . 

-26.6 

47 

1 335.4 

691 .9 

-75. 

14.9 

Q 

21 5. p 

606.2 

-365 . 

-26.6 i 

4 8 

1 363.4 

6 93.4 

-374, 

3.1 

10, 

243. 7 

692 . R 

-951. 

-137.7 

49 

1391,4 

683.4 

-596. 

-40.0 

n 

271.5 

6 8 9.? 

-856. 

-70.0 

50 

1419.3 

696 . 5 

-40 t . 

4.6 

! ?. 

?99 . 5 

6 9 7.3 

-729. 

-65.7 

61 

1447.2 

676.5 

-74 4. 

-74.9 

n 

327.5 

666.1 

-612. 

-18.6 

6? 

1476.2 

696.5 

. -674. 

0.0 

i <• 

365 . 5 

684 . 3 

-572 . 

-11 .6 

53 

1503. 3 

69 6 . 5 

-557. 

7.7 

15 

353.4 

691.2 

-253. 

26.9 

54 

1631.4 

696.5 

-501. 

-2.3 

16 

411.4 

69 0.3 

-504. 

-1.5 

55 

1 557.7 

707.1 

-496. 

-13.5 

17 

4-»0. 5 

679.7 

-618. 

-11.5 

56 

1585. 3 

707.1 

-437. 

-6.2 

IB 

467.5 

679.8 ' 

-891. 

-88.9 

57 

1597.4 

708.1 

-421 . 

-5.4 

10 

49 5.4 

676.4 

-724. 

-155.6 

58 

1625.3 

710.2 

-462. 

-4.6 

?n 

57 3.4 

67 5.1 

-189. 

3556.9 

59 

1683.4 

711.1 

-320. 

16.3 

21 

5 51.4 

672.1 

-28. 

19. T 

60 

1680.7 

716.3 

-387. 

9.5 

2 ? 

570.4 

677.0 

136. 

-9.5 

61 

1708.7 

720.1 

-573. 

-2.3 

23 

607.2 

669.0 

-449. 

-6.2 

62 

1736.5 

723.9 

-855. 

-11.4 

24 

675.2 

669.0 

-75. 

59. H 

63 

1764.1 

728.5 

-826. 

-16.0 

25 

6° 5. 9 

669.0 

-i ms. 

-253.4 

64 

1 702.0 

731.5 

-294. 

21.8 

26 

71Q.3 

670.6 

808.. 

322.4 

65 

1947.8 

737.6 

-33. 

• -0.3 

27 

747 . 3 

66O.0 

96 . 

' -372.1 






29 

775. 3 

668.1 

-545. 

-91.5 






2« 

803. 7 

670.6 

429. 

263.5 






?0 

? 5 1 . 3 

6 70.6 

-4 64. 

-1.6 






31 

B59.4 

670.6 

-42 8. 

21.7 






32 

837.4 

670.3 

-253. 

18.7 






33 

915.5 

670.6 

-18 2. 

10.2 






3* 

947,6 

670.6 

-283. 

-24.9 






35 

971 . 5 

672. 1 

-400. 

-59.6 






36 

999.6 

673.6 

98 9. 

204.5 






37 

1042.4 

678.7 

-472, 

16.5 






. 38 

103 3.5 

679.7 

-389. 

28.0 






39 

mi. 4 

681.8 

-802. 

-71. T 













magnetometer 

TRAVERSE R 1 * 

F 


MAGNETOMETER 

TRAVERSE F 1 

- F 









STATION 

STATION 

STATION 

VERTICAL 

STATION 

STATION 

STATION 

VERTICAL 

NUMBER 

DISTANCE 

ELEVATION 

anomaly 

NUMBER 

D I STANCE 

ELEVATION 

anomaly 


{ METERS1 

tMETERS) 

(GAMMAS) 


{ METERS t 

(METERS) 

(GAMMAS) 









t 

0.0 

71 3.2 

980. 

1 

0.0 

725.4 

-542. 

2 

30.4 

713.5 

9 80. 

2 

30.4 

72 5. 7 

-5 Q 2 • 

3 

61.0 

713.8 

530. 

3 

61.0 

726.0 

-531. 

4 

91.4 

714.1 

879 . 

A 

91,4 

7 2 6.3 

-521 . 

5 

171.9 

714.5 

1329. 

5 

121.9 

727.0 

-701 • 

6 

152.4 

715.1 

979. 

6 

152.4 

727.3 

-300* 

7 

132.9 

715.7 

1726. 

7 

132.9 

72 7.6 

380. 

B 

213.4 

716.3 

548. 

8 

213.4 

727.9 

-30. 

9 

243.8 

7 1 7 . 3 

1428. 

o 

2 43.0 

72 S. 5 

— 69. 

10 

274.3 

719.3 

1603. 

l 0 

2 74.3 

728.3 

- 309. 

1 1 

304.8 

7 2 0.9 

1727. 

1 1 

3 04.0 

72<3.l 

36 1. ! 

12 

335.3 

722.4 

1427. 

12 

335.1 

729.4 

34 1 . 

1 3 

365.3 

72 3.9 

1627. 

13 

365.8 

730.0 

502. 

14 

3 96.0 

725.4 

1176. 

14 

396.2 

730.3 

7. hi. 

15 

42-6.4 

727.0 

926. 

15 

426. T 

7 30.6 

Z?2 . 

16 

454.9 

728.5 

901 . 

16 

457.2 

730.9 

2 7 1 - 

1 7 

437.4 

731.5 

1025. 

17 

48 7. 7 

731.5 

643 . 

1 8 

517.2 

737.6 

3 75. 

15 

513.2 

734.0 

473. 

19 

547.1 

743.7 

850. 

to 

543.3 

736.7 

1244. j 

20 

5 7 7.0 

749.8 

66 5. 

23 

5 73.3 

739. 1 

794 . 

21 

606.9 

755.9 

304. 

21 

609.3 

741.9 

1244. 

22 

636. 7 

762.0 

724. 

22 

6’9. 5 

744.3 

1 444 . 

23 

666.6 

76 3.1 

9 7 4 . 

23 

670.0 

747.1 

94 5. 

24 

696.5 

774.2 

1 073- 

24 

700.4 

74°. 5 

995 . 

2 5 

726.3 

730.3 

1 248. 

25 

7 30*6 

752.2 

645. 





26 

761.1 

754.7 

74 6. 





? T 

7<n .6 

"75 7. 1 

84 6. 





28 

321 .7 

759.6 

74 6. 





29 

3 52.2 

762.0 

747. 





30 

382.4 

765.7 

947. 





31 . 

° l 2 . 9 

749.3 

1047. 





32 

943.1 

773.0 

1097 . 





33 

973.2 

776.6 

1048. 





34 

1003.7 

780.3 

1248. 








MAGNETOMETER 

TRAVERSE P * 

- F" 

STATION 

STATION 

STATION 

VERTICAL 

NUMBER 

DISTANCE 

ELEVATION 

ANOMALY 


(KETF^SI 

mereRSJ 

(GAMMAS) 

l 

0.0 

725.4 

-542. 

2 

30.4 

775.1 

-663. 

3 

M .0 

724.6 

-603. 

4 

9 1.4 

77 3.9 

-44. 

5 

121.0 

7 2 3.6 

-394. 

6 

152. 4 

723.0 

-634. 

7 

19 2.9 

722.4 

-745. 

n 

213.4 

721.8 

-606. 

Q 

243. g 

721.5 

-645. 

10 

2 74.3 

720.9 

-846. 

11 

304. 9 

7 2 C. 2 

534. 

12 

335.3 

719.6 

414. 

13 

V.5.8 

719.3 

684. 

14 

3 06.2 

716. 7 

-1 ) 7. 

15 

425.7 

718.1 

63. 

16 

457.2 

717.6 

103. 

17 

437.7 

71 7.2 

272. 

19 

518. 2 

716.6 

-4 8. 

10 

545.6 

716.0 

632. 

20 

5 70. i 

716.4 

1331. 

21 

509. h 

715.1 

471. 

22 

640.1 

714.5 

471. 

23 

6 70.6 

713.9 

9 91. 

24 

701.0 

713.2 

08 0. 



MAGNETOMETER 

TPAVF»SE g* 

G 

STATION 

STATION 

STATION 

VERTICAL 

NUMBER 

OT STANCE 

ELEVATION 

ANOMALY 


(METERS) 

( METERS) 

(GAMMAS) 

1 

0.0 

716.3 

450. 

2 

30.4 

717.5 

440. 

3 

61.0 

719.0 

379. 

4 

91.4 

720. 2 

4 49. 

5 

121.9 

72) . B 

759. 

6 

152.4 

72 3.0 

818. 

T 

15 2.9 

7?4. $ 

1078. 

8 

213.4 

728.7 

1158. 

9 

2 4 3.8 

727.3 

1228. 

10 

274.3 

723.5 

1 107. 

1 l 

304 . 0 

729.4 

1 C2 7. 

12 

. 320.1 

739.3 

1117. 

13 

338.3 

7 31.5 

1197. 

14 

36 5.5 

736.4 

1477. 

15 

39 5.6 

74 1.6 

1536. 

16 

425.8 

74 5.5 

1 366. 

17 

4 66.0 

751.6 

1016. 

t« 

486.2 

75 6.5 

566 . 

19 

51 6.3 

761.7 

936. 

20 

546.5 

756.6 

1 355. 

21 

576.7 

7 7 1 . B 

1 155. 

22 

606.9 

776.6 

1156. 

•23 

637.0 

781.8 

935. 

24 

66 7.2 

7 86.4 

161 8. 



MAGNETOMETER TRAVERSE G ,, ’~ g 


MAGNETOMETER TRAVERSE H - H* 


STATION 

MJMRER 


STATION 
D! STANCE 
( mf TERS) 


STATION 

ELEVATION 

( metersi 


VERTICAL 

ANOMALY 

(GAMMAS) 


1 

0.0 

2 

30.4 

3 

61.0 

4 

01.4 

5 

121.0 

ft 

1 S2.4 

7 

152.9 

R 

2 13.4 

9 

228. 6 

m 

2 4 1 . P 

1 1 

2 74.3 

12 

3 04 . 8 

13 

375.3 

14 

3&8.8 

15 • 

396.2 

tv 

4 26.7 

IT 

457.2 

I 8 

487.7 

in 

518.2 


716.3 

623. 

717. 2 

612. 

717.8 

732. 

71 8 • 7 

822. 

710.6 

86? . 

720.5 

94? . 

721.5 

021 . 

722.4 

1061. 

722.7 

1 363. 

724.3 

141 3. 

7? 4. 2 

1 ?o 3 • 

724.1 

1354. 

72;.. 0 

1 364. 

727.0 

12 84. 

727.0 

1324. 

728.8 

1345. 

720.7 

1365. 

73C.6 

136 5. 

7 31.5 

1365. 



STATION 

STATION 

NUMBER 

DISTANCE 
1 YPTCRS) 

1 

0.0 

? 

18.2 

3 

30.4 

4 

■' 45.7 

5 

60.7 

6 

90.8 

7 

171.0 

8 

151 .? 

9 

181.4 

10 

m.a 

11 

242.3 

1? 

272.5 

1 3 

3 02.7 

14 

332.8 

15 

3 63.0 


STATION 

VERTICAL 

ELEVATION 

ANOMALY 

< METERS) 

< GAMMAS) 

847.3 

1003. 

8 48. 3 

1803. 

8 40.? 

1 7R3. 

850.4 

1553. 

847.0 

1 8 « 3 . 

843.7 

1 7°2 . 

84 0. 3 

14 7?. 

A 37.0 

141?. 

8 33.6 

127?. 

mo. 6 

123?. 

827.5 

1021- 

8 31.2 

801. 

8 34.5 

731 . 

837.9 

711. 

841.3 

521 . 



MAGNETOMETER 

traverse g* 

- G' * 

STATION 

station 

station 

VERTICAL 

N'jmm 

distance 

elevation 

ANOMALY 


INTERS) 

(METERS) 

(GAMMAS) 

i . 

0.0 

716.3 

450. 

?. 

45.7 

7 1 8 . 4 

440. 

3 

01.4 

720. 2 

600. 

4 

137.2 

722.4 

1061. 




MAGNETOMETER 

TRAVERSE T 

- I f 


MAGNETOMETER 

TRAVERSE I * 

- I” 

scat ion 

STATtON 

station 

VERTICAL 

STATION 

STATION 

STATION 

VERTICAL 

NUM?E° 

DISTANCE 

ELEVATION 

ANOMALY 

NUMBER 

distance 

elevation 

ANOMALY 


(METERS) 

(ME7ERS1 

t GAMMAS ! 


(MFTFRS1 

( METERS) 

(GAMMAS) 

l 

0-0 

56 1.1 

1601. 

t 

0.0 

914.4 

-5°3. ' 

2 

30.4 

559. 5 

1551. 

2 

45.7 

9 14.4 

40 T. 

7 

61.0 

55 7.4 

1401. 

\ ' 7 

91.4 

9 14,4 

508. 

4 

°1.4 

566.7 

1 262. 

4 ' 

177.2 

914.4 

299. 

5 

121 

063. 1 

1 3 22. 





6 

\ 52. 4 

R51.D 

1352. 





r 

137.9 

843.9 

1 72 2. 





r 

213.4 

8 4 4.7 

1173. 





9 

223.6 

8 44.6 

1143. . ) 





10 

2 34.9 

8 4 7.5 

1047. ’• 





1 l 

2 79.9 

R4D.7 

373 . 


MAGNETOMETER 

TRAVERSE I” 

!*• 

12 

24 2.0 

8 3 8.2 

744. ■ - 





13 

257.3 

836.7 

7*4 4 . 





t4 

2 72.5 

R 3 5 . 2 

774. 





15 

30"3 . 0 

R73.6 

754. 

STATION 

STATION 

station 

VERTICAL 

16 

2 ? 7 . 5 

R ? 7 . 1 

724. 

NUMBER 

DISTANCE 

elevation 

ANOMALY 

17 

363.9 

8 7 0.6 

684. 


(METERS) 

( METFRS) 

(GAM«AS> 

1 3 

3 94.4 

637.1 

675. 





1® 

424.9 

333.6 . 

685. 

1 

0.0 

S83.9 

-321, 

20 

4 55.4 

8 3 5.2 

605. 

2 

29.6 

3 9 L . 5 

149. 

?! 

465.9 

336.7 

435. 

3 

59.2 

809 . 2 

74 ft. 

77 

516.3 

638.2 

505. 

4 

9 8.7 

0 06.8 

528. 

?3 

546.5 

841.9 

295. 

5 

118.3 

. 914.4 

29ft. 

?4 

57 6. 7 

845. 9 

2 1 6 . 





25 

606.9 

849. 8 

2 6. 





26 

67 7.0 

R 5 3 . 4 

-144. 



- 


7 7 

66 5.7 

P.63.8 

-274. 





7 3 

6 >4. 3 

877,9 

-744. 





29 

77 3.0 

887.9 

-234. 





30 

752.9 

390.0 

-203. 


MAGNETOMETER 

TRAVERSE J' ’ 

- J 

31 

767.7 

896.1 

-22 3. 





3? 

510. B 

908.3 

-713. 





33 

< , 540.6 

914.4 

-593 . 









station 

STATION 

STATION 

VERTICAL 





NUMBER 

DISTANCE 

ELEVATION 

ANOMAl Y 






(METERS) 

( ME I hRS ) 

< GAMMAS) 



” 


1 

0.0 

R99.8 

761 . 





2 

41.0 

8 49.2 

8 51. 





3 

91.5 

899.? 

661, 





4 

121 .9 

898.9 

1 57 1 . 





5 

i 213.4 

<3 9. 7.6 

1 237. 





6 

30 3.0 

866.7 

714. 



MAGNETOMETER TRAVERSE I - J 


STATION 

STATION 

STATION 

VERTICAL 

NUMBER 

DISTANCE 

ELEVATION 

ANDMAl y 


(METERS) 

IMF Tf RSI 

( GAMMAS 1 

I 

0.0 

P M . 0 

1601. 

2 

30.4 

858.9 

1701. 

3 

6 0.0 

86 7.1 

1.59 0, 

4 

91.4 

856.3 

1400. 

5 

1 21.0 

853.4 

1310. 

6 

157.4 

.9 6 7.5 

1 40 0 . 

T 

182.8 

86 1.9 

1330. 

R 

213.3 

651.0 

134 0. 

o 

7 4 3.8 

8 6 0.1 

5 270. 

10 

274.3 

849.2 

1 190. 

11 

304. ft 

84ft. 6 

1179. 

. 1? 

335.2 

84 7. 7 

1170. 

13 

3 65.7 

H 4 6 . 7 

1030. 

I* 

304.2 

84 5.8 

1018. 

IS 

426.7 

644.7 

948. 

I i> 

441.9 

84 7. 3 

938. 

1 7 

467.2 

84 8. 3 

ftSft. 

t P 

472.4 

843.9 

33ft. 

10 

487. 6 

649.6 

85 8. 

.70 

518.1 

850.4 

82 7. 

21 

533.4 

851.3 

817. 

22 

548.6 

851.9 

577. 

.7 3 

570.1 

8 4?. ft 

617. 

• 74 

600.6 

ft 5 3. 4 • 

557. 

25 

640.0 

8 4 6.2 

7 L 7 . 

26 

6 70. 5 

858.9 

676. 

71 

701-0 

8M .2 

4 76. 

?P 

731.5 

864.4 

676. 

7 *> 

767.0 

6 67.2 

626. 

30 

797.4 

869.0 

51 l . 

31 

8 22.9 

872.6 

515. 

32 

8 53.4 

875.4 

255. 

33 

8 83.9 

8 78. 1 

415. 

34 1 

014.4 

880.9 

45. 

75 

044. ft 

8 R 3 . 6 

ISO. 

36 

1005.8 

886.7 

714. 



magnetometer 

TRAVFRSE J 

-<• ' » 

STATION 

STATION 

station 

vertical 

NUMBER 

or stance 

ELEVATION 

anomalt 


( mf rr-RS ) 

{METERS) 

(GAMMAS) 

1 

0.0 

8 86.7 

714. 

7 

40.0 

903.7 

9 34. 

3 

118.5 

844.6 

924. 

4 

1 79.5 

897.6 

BO 9. 



M ACNE TOHfi TE R 

TRAVERSE 

K- K» 

STATION 

station 

STATION 

VERTICAL 

NUMBER 

(MS 1 '.Wf, K 

ELEVATION 

AN OK A L Y 


(MrTcNSl 

(ML TER SI 

( GAMMAS! 

1 

0.0 

59 3.5 

33 H. 

2 

304 . d 

54 j, 5 

44£ . 

3 

609.6 

5 8 4.3 

3 5 7. 

4 

914.1 

5 94. r 

517. 

5 

1218.4 

5 91.6 

396 . 

6 

1523.4 

614.0 

765. 

7 

1675.5 

62'). 5 

415. 

8 

1 4 1 9 .3 

67 3.3 

1014. 

0 

1445.5 

62 4.4 

714. 

10 

70 71 . 7 

4 7 '■>.•♦ 

554. 

11 

2754.0 

6 12. J 

6 6 3. 

12 

2294.5 

6 1 J . 9 

7 0. 

13 

2 ~j 3 5 . 1 

639.8 

937. 

14 

27 70.9 

642.5 

77 1. 

15 

3075.4 

6 3 1.6 

629. 

16 

3340.2 

6 10.0 

5 7a. 

17 

3685 . 0 

623.1 

2 0. 

18 

3 9 e i , 4 

64 7.1 

442. 


CO 

VT 



macnetume ter TRAvrRsr k 


MAGNETOMETER TRAVERSE <•»•- K* * 


• - K* * 


STATTCN 

NUMBER 


1 

? 

3 

* 

5 

6 
7 
9 
o 

to 
11 
12 
i 3 
1 * 
15 
1 6 
1 T 
19 
19 
?0 

71 

72 
? 3 
?* 

25 

26 
27 

78 

79 
30 
?1 

72 

73 


ST AT I OM 

STATION 

VERTICAL 

STATION 

STATION 

STATION 

VERTICAL 

DISTANCE 

ELEVATION 

ANOMALY 

NUMBER 

oi stance 

ELEVATION 

anomaly 

(METERS! 

[METERS) 

(GAMMAS! 


( ME TER SI 

(METERS! 

(GAMMAS) 

0.0 

6 A 7 . 1 

*62. 

1 

0.0 

897.6 

809. 

30A. 5 

660. 9 

211. 

2 

*7 . * 

903. 7 

72*. 

609.1 

6 6 6.9 

320- 

3 

60. * 

906.8 

72* . 

9 1 A. 1 

676.7 

200. 

* 

121.0 

501.9 

5*5. 

1219.3 

692.5 

22*. 

5 

136.3 

900.9 

*15. 

1 VI9.9 

702.6 

81 . 

6 

151.5 

899.8 

5 25. 

1359.1 

699.8 

59. 

7 

211.2 

887.3 

*3 5. 

1370. A 

699.3 

32. 

fl 

272.7 

887.0 

*15- 

1395.6 

6°9. 3 

12. 

9 

337.2 

897.6 

3*6. 

1 A .10.° 

6 9 9 . n 

-73. 

1 0 

3'>3.2 

001.9 

316. 

1*16. L 

699.8 

-191. 

LI 

*63,9 

897.0 

216. 

1*31.3 

699.8 

-93. 

12 

51 a. 5 

391.5 

206 . 

1 AA6. 6 

690. * 

-187. 

1 3 

575.2 

885. B 

1*6. 

1 AS] . 9 

701.3 

-312. 

' l* 

6 ■* A - 

S ?/,,! 

t * '> . 

1 *77.1 

702.9 

-382. 

15 

69 A. 3 

359.2 

26. 

1*92. 3 

7 0 A . A 

-277. ■ 5 • 

16 

755.3 

862.0 

207. 

1 9 n /. . 

7 09 . f, 

-707 . 

17 

816.0 

855.6 

117. 

1670. 7 

71 3. A 

-106. 

IB 

8 70.6 

861. A 

— U t f> m 

1570.9 

716.9 

-27. 

19 

8*6.* 

E5S.9 

-f? 5 ♦ 

! 6R-3.2 

716.3 

-120. 

20 

91 *. A 

856.7 

-136. 

1 t. 6 

7 1 7. 11 

- 1 ;’o. 

71 

rw.? 

8 5 7.1 

- 5 57. 

1675.8 

719.3 

. -8*. 

72 

997.9 

880. * 

9 3. 

1630.7 

720.9 

-18*. 

23 

10*2.7 

3*1.9 

12 P.. 

1703.7 

713.7 

-170. 

2* 

1103.7 

338.5 

23. 

1 7 1- 2 . 7 
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APPENDIX C 


Accumulated Density Determinations 
for Rocks from the Katmal Area 


Class if icatlon 

Reference 

Rock Type 

Location 


Ident 5 ty 

Naknek sediments 

Klenle, 1969, p. 63 

Naknek sandstone 

Grosvenor Lake 


KG-71 

Naknek sediments 

Klenle, 1969, p. 63 

Naknek marly shale. 

Grosvenor Lake 


KG-33 

Naknek sediments 

Klenle, 1969, p. 63 

Naknek shale 

Grosvenor Lake 


KG-39 

Naknek sediments 

Kienle, 1969, ,p. 63 

Chisik conglomerate 

Grosvenor Lake 


KG -69 

Naknek sediments 

Klenle, 1 969 , p. 63 

Naknek sandstone 

Grosvenor Lake 

/ 

/ 

KG-72 

rhyol i t i c glass 

Curtis, 1968, p. 192 

rhyol 1 te g 1 ass 

Novarupta dome 



mixed lava 

Klenle, 1571 , pers.com. 

porous banded lava 

Novarupta central dome 


KN-3-1-&5 

mixed lava 

Klenle, 1971 , pers.com. 

porous banded Java 

Novarupta east margin 


KN-4-1-65 

mixed lava 

Kienle, 1971, pers.com. 

porous banded lava 

Novarupta east margin 


KN-4-2-65 

andcs Stic lava 

Kienle, 1971 , pers . com. 

andes 1 te 

Fa 1 1 i ng Mtn . 


KFm-1-70 

andesitic lava 

Kienle, 1971 , pers.com. 

andes i te 

Falling Mtn. 


KFM-la -70 

andesitic lava 

Klenle, 1571 , pers.com. 

andes i te 

Fissure Dome 


KFD-l-70 

andcs I tic lava 

Kienle, 197', pers.com. 

andes i te 

Fissure Dome 


KFD- la-70 

andes i t i c iava 

Kienle, 1269, P- 63 

porous andes 1 te 

Mt. Trident 


KT-28 

endes i tic lava 
teph ra 

Kienle, 1 969 , p. 63 
Griggs, 1922, p. 293 

basaltic andes i te 
ash 

Knife Peak 


KKP2-1 

teph ra 

Lovering, 1957, p. 1593 

fumarolic altered ash 

Lower Valley, Fumarole #1, 

zone 1 

1 T L 53 

teph ra 

Lovering, 1957, p. 1593 

fumaroltc altsred'ash 

Lower Valley, Fumarole #1 , 

zone 2 ■ 

2TL53 

tephra 

Lovering, 1957, P. 1593 

fumarolic altered ash 

Lower Valley, Fumarole #1, 

zone 3 

3TL53 

tephra 

Lovering, 1957. P- 1593 

f umaro tic altered ash 

Lower Valley, Fumarole #1 , 

zone 4 

47L53 

tephra 

Lovering, 1957, p. 1593 

fumarolic altered ash 

Lower Valley, Fumarole If 1, 

zone 5 

5TL53 

tephra 

Lovering, 1957, P. 1 59 3 

funarollc altered ash 

Lower Valley, Fumarole ff 1, 

zone 6 

6TL53 

tephra 

tephra 

Lovering, 1957, P* 1593 
Curtis, 1 963 , p. 207 

normal ash 
tephra Layer D 

Lower Valley, Fumarole H\ , 
upper southeast branch 

zone 7 

7T153 

tuff 

tuff 

Curtis, 1963, p. 133 
Curtis, 1968, p. 183 

tuff 

upper tuff unit 

Middle Valley, Knife Creek 
Middle Val ley 

Gorge 


tuff 

Curtis, I960, p. I 83 

tuff 

southern branch. Fissure Lake 


tuff 

tuff 

tuff 

Curtis, 1 968 , p. 184 
Curtis, 1968, p. 183 
Curtis, 1968, p. 194 

dense tuff 
dark gray tuff 
center of d i ke- 1 f.ke 

southeast branch, Explosion pit 
upper southeast branch 
upper southeast branch 



body of tuff 


Densi ty . 
(gm/cc) 

2.60 

2.60 

2 . 1*8 

2.75 

2.75 

2.25 

to 

2.30 

2.22 

1.70 

1.81 

2.51 

2.47 

2.49 
2.55 
2.44 
2.47 

1 .03 

0. 75 

1 . 10 
1.18 
1.05 

0. 84 

1. C0 
0.93 
1.09 

1 . 15 
1.35 

1.50 
1.73 

to 

1.85 

1.60 

to 

1.72 

1.75 



Class i f teat ion Reference 


Rock Type 


Location 


Ident 1 ty 


tuff 

Kienle, 1371, 

pers . com. 

tuff 

tuff 

Kienle, 1371, 

pers.com. 

tuff 

tuff 

Kienle, 1971, 

pers . com. 

tuff 

tuff 

Kienle, 1371, 

pers . com. 

tuff 

tuff 

Kienle, 1971, 

pers . com. 

tuff 

tuff 

Kienle, 1971, 

pers .com. 

tuff 

tuff 

Kienle, 1971, 

pers , com. 

tuff 

tuff 

Kienle, 1971 , 

pers . com. 

tuff 

tuff 

Kienie, 1971 , 

pers . com. 

tuff 

tuff 

Kienle, 1971 , 

per 5 . com. 

tuff 

tuff 

Kienle, 1371, 

pers.com. 

tuff 

tuff 

Kienle, 1371, 

pers.com. 

tuff 

tuff 

Kienle, 1971, 

pers.com. 

tuff 

tuff 

Kienle, 1971 , 

pers.com. 

tuff 

tuff 

Kienle, 1971, 

pers .com. 

tuff 

tuff 

Kienle, 197' , 

pers . com. 

'tuff 

tuff 

Kienle, 1971, 

pers . com. 

tuff 

tuff 

Kienle, 1971. 

pers . com. 

tuff 

tuff 

Kienie, 19/0, 

p. 6 64ii 

ash 


Kiddle Valley, Knife Creek Canyon 
Middle Valley, Knife Creek Canyon 
Middle Valley, Knife Creek Canyon 
upper south branch, Corner Lake 
upper south branch, Fissure lake 
upper south branch, Fissure lake 
upper south branch, Fissure Lake ? 
upper south branch, Fissure Lake 
upper south branch, Fissure Lake 
upper southeast branch, Peckish Springs 
upper southeast branch. Peckish Springs 

upper southeast branch, Peckish Springs 

upper southeast branch, Peckish Springs 

upper southeast branch, base of Knife 

Creek Glacier 

upper southeast branch, base of Knife 
Crock Glacier 
Novarupta Basin 
Nova rupta Basin 
Novarupta Basin 
Lower Valley 


KPFG-1-70 

KPFG-la-70 

K?FG- 1 b-70 

KPFC-1-70 

KPFll-1-70 

KPFll-2-70 

kpfli-3-70 

KPFU-4-70 

kpfu-5-70 

KPFP- ! -70 

KPFP-2-70 

KPFP-4-70 

KPFP-5-70 

KPFK3-13-70 

KPFK3*lb-70 

KPFN-2-70 

kpfn-3-70 

KPFN-4-70 
5 sor. plus ■ 

throughout 

section 


Dens 1 ty 
(gm/cc) 


1 .62 
1-53 
1.53 
1.52 
1.61 
i -63 
1.61 
1.66 
1 .67 
1.57 
1 .62 
1.73 
1.85 
1.73 

1.61) 

1.33 
1 .31* 

1 .36 
1.03 
average 
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APPEND! X D 

Accumulated Seismic Refraction Data 


Parameters of Seismic Lay 




Layer 

Layer 

i nvest igator 

Prof i le 

Location 1 

2 

Hatumoto 


Lower 

2.4 

and Ward( I9&7) 

Valley 

0.3 

Kienle (I97I) 

E-E’ 

lower 

2.0 



Valley 

0. 38 

Klenle (1971) 

H-H' 

Middle 

1.0 



Valley 

? 

Klenle (1971) 

F-F 1 

Southeast 

1.0 



Branch 

0. 30 

Sbar and 

Prof 1 1 c 1 

Southeast 


Motomoto( 1970 Sec t ion 7 

Branch 


Sbar and 

Profile 1 

Broken Mtn. 


Hatumoto 11 

Section 1 

Valley 


Sbar and 

Profile 1 

Broken Mtn- 


Hatumoto " 

Section 2 

Va 1 ley 


Sbar and 

Profile t 

Broken Mtn. 


Hatumoto 11 

Section 3 

Valley 


Sbar and 

Profile 1 

Broken Mtn. 


Hatumoto 11 

Section 4 

Valley 


Sbar and 

Profile 1 

Broken Mtn. 0. 5 

2.0 

Hatumoto " 

Section 5 

Valley 0.1 

0.3 

Cedney (19?0) 

B-B 1 

Broken Htn. 

4.4 


Profile 5a 

Valley 

0.64 

Gedney " 

B-B' 

Broken Mtn. 

4.7 


Profile 5b 

Valley 

0. 38 

Klenle (1970) 

C-C' 

Broken Mtn. 2. 5 

5.0 


Profile la 

Valley' 0 . 29 

0.60 

Klenle " 

C-C’ 

Broken Mtn. 1.5 

7.0 


Profile lb 

Valley 0.2? 

0.43 

Klenle 11 

C-C’ 

Broken Mtn. 2.0 

7.0 


Profile 2a 

Valley 0.17 

0.34 

Klenle " 

C-C* 

Broken Htn. 1 .0 

1.5 


Prof lie 2b 

Volley 0.17 

0.37 

Klenle " 

C-C 

Broken Mtn . | , 0 

2.0 


Prof i le 3a 

Valley 0.13 

0.37 

Klenle M 

C-C 

Broken Mtn. ) .5 

6.0 


Prof i lc 3b 

Valley 0.20 

0.31 

Sbar and 

Profile 2 

Southern 


Matumoto(l37>)Section 7 

Branch 


Sbar and 

Prof i le 2 

Southern 


Hatumoto " 

Section $ 

Branch 


Sbar and 

Prof i le 2 

Southern 0.5 

3.1 

Hatumoto “ 

Section 5b 

Branch 0 . 13 

0.27 

Sbar and 

Profile 2 

Southern 0.7 

4.7 

Hatumoto " 

Section 5a 

Branch 0.13 

0. 40 

Sbar and 

Prof i le 2 

Southern 1.2 

1.6 

hatumoto " 

Section 4b 

Branch 0.13 

0. 45 

Sbar and 

Prof i le 2 

Southern 0.6 

8.3 

Hatumoto “ 

Section 4a 

Branch 0.14 

0.25 


thickness (meters) 
ers v a Co c iW iKx/ o BeS 


Layer 

Layer 

Layer 

Layer 

Layer 

Total 

3 

4 

5 

6 

7 

Thickness 

7.3 

36. 



7 

45.7 

0.5 

1.00 



3.8 



44.0 

? 



46.0 


.99 

2. 29 




28.0 


34.0 

7 


63.0 

0.66 


1.79 

2.53 



18.0 


? 

/ 


19.0 

0.67 


1.89 




22.0 


? 



22.0 

0 . 61 


2.2 


- 


25.0 


7 



25.0 

0.62 


1.8 




22.0 


7 



22.0 

0. 56 


1.9 




25.0 


7 



25.0 

0.61 


1.8 




22.0 


7 



22.0 

0.52 


2.1 




24.0 


7 



26.5 

0.44 


2.0 




31.0 


7 



35-4 

0 . 62 


2.25 




22.0 


7 



26.7 

0. 65 


1.75 




25.0 


58.0 

7 


90.5 

0.71 


2.10 

3.36 



23.0 


72.0 

7 


103.5 

0. 71 


2.10 

3.2$ 



22.0 


73.0 


7 

104.0 

0.62 


2.09 


3.71 


31.0 


46.0 


7 

79.5 

0.62 


2. 09 


2. ?1 


32.0 


42.0 


7 

77.0 

0.62 


2.44 


2.68 


32.0 


42.0 


7 

81.5 

0 . 43 


2.44 


3.68 



35.0 

? 



35-0 


0.8? 

2.0 





50.0 


7 


50.0 


1.1 


2,5 



16.6 

42.1 

7 



62.3 

0.59 

1.2 

1.8 




45.8 


40.8 


7 

92.0 

0. 74 


1.8 


3.8 

7.3 

56.3 


? 


66.4 

0.62 

0. 87 


2.1 



25.9 



7 


34.8 

0.65 



2.5 







- 



Layer 

Layer 

Investigator 

Profile 


Location 

1 

2 

Sbar and 


Prof i le 

2 

Southern 


30.0 

Matumoto( 

97 1 )Sec t ion 

3 

Branch 


0.30 

Sbar and 


Prof i le 

2 

Southern 


25.0 

Matumoto 

It 

Sect ion 

2b 

8 ranch 


0.32 

Sbar and 


Prof i 1 e 

2 

Sou thern 


25.0 

Matumoto 

1 t 

Section 

2a 

Branch 


? 

Gcdrtey {1970} 

A-A' 


Southern 


1.4 



Prof i le 

4a 

Branch 


0.25 

Gedney 

tt 

A-A’ 


Southern 


23.1 



Prof i 1 e 

4b 

B ranch 


0.35 

Gedney 

H 

A-A* 


Southern 


3-7 



Prof i le 

3a 

Branch 


0.33 

Gedney 

1 1 

A-A* 


Sou them 


1.6 



Prof i lc 

3b 

B ranch 


0.38 

Gedney 

M 

A-A' 


Southern 


1.9 



Prof i le 

2a 

B ranch 


0.54 

Gedney 

1 1 

A-A' 


Southern 


4.5 



Profile 

2b 

Branch 


0.50 

Gedney 

II 

A-A' 


Sou them 


2.1 



Prof i le 

la 

Branch 


0.25 

Gedney 

1 1 

A-A' 


Southern 


5.3 



Prof i 1 e 

lb 

Br anch 


0.33 

Sbar and 


Prof i 1 e 

3 

Novarupta 



Matumot o{l97l) Sect ion 

2 

Bas i n 



Sbar and 


Prof i le 

3 

Novarupta 



Matumoto 

1 1 

Section 

3a 

Basin 



Sbar and 


Profile 

3 

Novarupta 



Matumoto 

1 1 

Sec t i On 

3b 

Basin 



Sbar and 


Profile 

3 

Novarupta 

0.9 

1.2 

Matumoto 

M 

Section 

4o 

Bas i n 

0.11 

0.42 

Sbar and 


Prof i le 

3 

Novarupta 

1.1 

1.4 

Matumoto 

1 1 

5ection 

4b 

Bas 5 n 

0.12 

0 . 59 

Sbar and 


Prof! 1 e 

3 

Novarupta 



Matumoto 

1 1 

Section 

5a 

Bas i n 



Sbar and 


Prof i 1 e 

3 

Novarupta 



Matumoto 

1 1 

Sect ion 

5b 

Bas i n 



Sbar and 


Prof i le 

3 

Novarupta 



Matumoto 

II 

Section 

6a 

Bas I n 



Sbar and 


Prof i le 

3 

Novarupta 



Matumoto 


Section 

6b 

Bas i n 



Sbar and 


Prof i le 

3 

Novarupta 



Matumoto 

n 

Section 

7 

Bas i n 



Klenle (1970)0-0' 


Novarupta 


5.0 



Section 

a 

Bas 1 n 


0.36 

Klenle ( 1970) D- D * 


Novarupta 


2.5 



Section 

b 

Bas i n 


0.49 


Layer 

Layer 

Layer 

Layer 

Layer 

Total 

3 

4 

5 

6 

7 

Thickness (1 




? 


30.0 




3.2 






? 


25.0 




2.5 




? 




25.0 


1.0 





47.0 



? 


48.4 

0.60 



3.0 





7 



23.1 



2.1 




43-1 



7 


46.8 

0. 70 



2.87 



13.8 

24. 1 

7 



39.5 

0.43 

0.8? 

2.15 




24.3 


45.6 

? 


71.8 

0.65 


1.42 

2.0 



30.6 


32.7 

7 


67.8 

0.60 


1.8? 

2.98 




35.6 


7 


37.7 


0.80 


2.62 




59’. 7 


7 


65.0 


0.95 


2. 75 



20.0 


37.0 

7 


57.0 

? 


1.8. 

3.1 





42.0 


7 

42.0 



? 


4.6 


17.0 

62.0 

7 



79.0 

? 

1.2 

2.3. 




4.9 

? 




7.0 

0.61 

0.80 





73.3 




7 

75.8 

0.65 




4.4 


12.0 


? 



12.0 

? 


1.5 


- 


10.0 

? 




10.0 

? 

1.1 





15.0 

1 




15.0 

? 

.96 





15.0 

35.0 

7 



50.0 

? 

1.1 

2.3 




15.0 

? 




15.0 

7 

1 . 1 






16.0 

73.0 

7 


94.0 


0. 8? 

1.55 

2.84 




43.O 

33.0 

7 


78.5 


0.8? 

1.55 

2.84 
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APPENDIX E.l 

DEVELOPMENT OF FUNDAMENTAL EQUATIONS 


This mathematical derivation is based on the methods presented 
by Hintler et al., (1962) and Sharma (1966). 

Consider a right-handed coordinate system: x,y,z, such that 
z is positive downwards. Let the unit vectors in the x,y,z, directions 

AAA 

be i,j,k, respectively. 

If V = AxAyAz represents a volume element at (x,y,z) with magnetic 

moment y, then its magnetization, M, is given by: 

Tt _ JT = 1 

M “ V AxAyAz M ’ 

s 

where: 

M — Mi + M j + M k . 
x y J z 


The vector, r, from position (A,B,C) to the element V at (x,y,z) 
is given by 

T* = (x-A) i + (y-B)j + (z-C)k , 

or letting: 

(x-A) = X, 

(y-B) = Y, and 

(z-C) -Z, 

r* = Xi + Yj + Zk, 

which yields: 

/ 2 2 7 
r «/ X + Y + Z . 

The magnetic potential P at (A,B,C) due to V is 


P 


M r 


r 3 
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which, upon substitution and solving, becomes 


P * 


f XH + YM + ZM 
x y z 


(X 2 + Y 2 + Z 2 ) 




AxAyAz 


Consider a body of cross section AxAz which is infinite In the y- 
direction. The potential at (A,B,C) Is given by 


r IMX + HZ 
P = AxA zf (— 


M Y 


,1 


3 / ? Ay 

J-»Wx 2 + v 2 +z 2 ) d /2 (X 2 + y 2 +z 2 ) /2 


solving : 


P = 3AxAz 


M X+M 1 
x z 

X 2 +Z 2 


The negative gradient of the potential in any direction gives 
a corresponding component of magnetic field strength. Thus the vertical 
field strength, V is found by 


V - 


3 P 

3Z * 


V - 3AxAz 


3XZ M ~M (X 2 - Z 2 ) 
x z 


(X 2 + Z 2 ) 

Likewise, the horizontal field strength in the x-direction is given 

by 


H 


-3P 
3x * 


H = 3AxAz 


3XZ M +M (X 2 - Z 2 ) 
z x 


(x 2 +z 2 ) 2 


-3P 


Clearly, = 0 and there is no component of field strength in the 

y-direction. This result follows intuitively from symmetry considerations 
Similarly, one sees that locations (A,B,C) and (x,y,z) are fully 
represented by (A,C) and (x,z), respectively. Also, the magnetization 
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In the y-direction Is inconsequential. 

Assuming that the body causing the anomaly extends from to 


in the z-direction, \l can be expressed as 


/ 


V -8 Ax I |DXM 
C 1 


— -MX 2 

x l(x 2 +z 2 ) I z 


1 


(x 2 +z 2 ) 


+ M 


.2 \ 


(x 2 +z 2 ) 2 


dz; 


which reduces to 


V * 3AX 


Similarly for H: 


H = 3AX 


-XM -ZM 
x z 

Z 2 +X 2 


ZM -XM 
x z 

Z 2 +X 2 


If the magnetization M, has only a component induced by the 


earth's field, then 


M - k F , 

whhere k is the magnetic susceptibility of the causative body. Clearly, 
the components of the magnetization are given by 

M. - kF cos X sin S, and 
x 

M = kF sin X , 
z 

where, X is the regional inclination of the earth's field, F is the 

magnitude of the regional total field, and S is the clockwise angle 

from magnetic north to the positive direction of the strike, +y, 

which is determined by correspondence with x being positive with increasing 



a . 


C\ : — 


x 



Figure E.la. Cross-sectional view of the modeling wedge. The initial 
traverse station corresponds with the origin. 

Figure E.lb. Overhead view of the modeling wedge with corresponding 
traverse designations. 

Figure E.lc. Components of the geomagnetic field* F. 

Figure E.ld. Components of the anomalous total field, T, in the case 
where V and H are small relative to F. 
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station distance along the traverse. Refer to Figure E.2. 

A A 

If it is assumed that the anomalous total field (T - Vi + Hk) is 

i-\ A „ A 

in essentially the same direction as the earth's field (F « F^i+F^j+F^k) 
ie, that V and H are small relative to F, then the magnitude of the 
anomalous field is approximately equal to its projection onto F: 


|f) /f '2+f"2+f“2 

1 x y z 

or referring to the geometrical relationships illustrated in Figures 
E.3 and 


T = H cos I sin S + V sin I . 


Ultimately, V, H, and T can be expressed in terms of their fundamental 
components : 


V = -3kFAx 


(x-A) cos I sin S + (z-c) sinl 


H - 3kFAx 


T = SkFAx 


(x-A) 2 + (z-C) 2 

(z~c) cos I sin S - (x-A) sin I 
(x-A) 2 + (z-C) 2 


(z-c) (cos 2 I sin 2 S-sin 2 I) ~ 3 (x-A) sinl cos I s i n S 


(x-A) 2 + (z-C) 2 


1 


Where x is the midpoint of increment x and z varies from c^ to c^; 
(A,C) is the station location. 

The vertical gradient of the total field anomaly, G, can be found 
by: 

G = 3T /Sz ; 
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which reduces to: 


G = SkFAx 


(cos 2 I sin 2 S - sin 2 I ) [(x-A) 2 - (z~C) 2 ] 
[ (x-A) 2 + (z-C) 2 ] 2 


3 (x-A) (z-C) s i nX cosl sin S 
[(x-A) 2 + (z-C) 2 ] 2 



These expressions for the various anomalies can be greatly simplified 
by letting: 

L - 3k F Ax 
R^= /{ x- A) 2 + (c^'-C)" 2 
R^= /(x-A) 2 +(c 2 “C) 2 
Hl= (x-A) 

M2= (x-A) 

N 1= (C^C) 

N2= (C 2 -C) 

H = H2-M1 
N = N3-N 1 
P = cosl sin S 
Q = sinl 

With the above relationships one obtains: 

H = L (NP-MQ) 

V =-L (MP+NQ) 

T = L[N(P 2 -Q 2 )“2MPQ] 

G - L [M2 2 -Ml 2 )- (N2 2 -N1 2 )] (P 2 -d 2 )+ A (M2N2-M1N1 )PQ . 
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In practice, the anomalies at position (A,C) are determined by solv- 
ing the previous equations for the contribution from each elemental strip, 
Dx X ( c2 - c 1 ) , and summing. This technique can be extended to cover cases 
of several bodies with varying magnetic properties, dimension, shape, 
position and strike, and any number and location of station (provided 
they are exterior to the magnetic bodies). 

A program has been written for the IBM 3&0/A0 to determine V,T, and 
G in this manner; i.e., at evenly spaced stations along straight line 
traverses over level bodies of paral lei orientation and speci fied depth, 
width and thickness and susceptibility, uniformly magnetized in the pre- 
sent regional field. A sample of the program utilized in determining 
the model anomalies is presented in Appendix E.2, 

There are a few limitations of this program which deserve attention, 
although most are clear from its derivation. 

It was found that altering the number of elemental strips within a 
body of arbitrary width does not markedly effect the resultant computed 
anomaly. But, as a rule of thumb, the width of the elemental strips 
should be at least half the station spacing. 

The relative position of each station above the anomalous body can 
be very important, especially i f the body is narrow, and more particularly 
if it is also shallow. Obviously, as the body narrows, so does its ano- 
maly; and as it nears the surface, Its effect is more pronounced. So, 
station location and spacing become critical for realistic representation 
of a shallow and/or narrow body. Carelessness in this detail of the 
modeling can cause the profile of the anomaly to be misshaped and in- 
correct; or, the anomaly missed entirely. Thus, the minimum model width 
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should dictate the maximum station spacing. Also, since the elemental 
strips are in effect, narrow bodies, It is wise to have the stations po- 
sitioned directly over the center of any elemental strip, and definitely 
not alternately at the center and then the edge of elemental strips. 

Station spacing is also critical in computing the anomaly due to a 
thin body in order to avoid producing erroneous profiles. Basically, 
the station spacing should be less than the thickness of the anomalous 
body to be modeled. 

Computation of this type of two-dimensional model is common in 
analysis of ocean-floor magnetic anomalies. Of course, such analysis 
is much simplified since essentially al 1 pa rame te rs are defined except 
the magnetization. The modeling program presented here could be used 
for such analysis wi thout modification. 



APPENDIX E.2 


Sample Computer Program of Two Dimensional Magnetic Analysts 


f . . 

C PRF'M’gD TN 1 « 57 C RY MARLA CAVS TPlBlF 

C***«e#*f-«-********«>>*«'<'*** 4 ‘<»****#*i>«‘«*#e**#^*aee<>*i***e******(<««** 

r 

C jv-wn-jCL IDENTITY 

r. xsjNirtAL rr -jr rooy along traverse* metfrs 
c OO c <>TM TO TOO OP OOOYtM 
C m-^vr'-i T I CAL EXTENT OF 800Y,M 
r VI*‘ ,r J T 7 ONTAL FXTFNT. M 
r. l)X= INCREMENTAL WIDTH, M 
C .HisT-'TU IX IN AnOY 

C s v = S T ?. !< I OP 3 0 : 3Y, CM prom TP.UPNORTHt DEGREES 
C X* = SI'SCE^T I ^!L T TY, ( 131-6 EMU/CC 
C 

C VH>= VERTICAL ANOMALY, gammas 
r T(0* TOTAL ANHMAL Y ! N GAMMAS 
C GIT ) = V c R T . GRAG. OP TOTAL ANOMALY ( GAMMAS 
C 

C SS*S r P!<P OF TRAVFRSF, MFC, CW OF TN 
C Sf-STATIOV SPACING* M 
C K! >*S T AT fm,N OISTANCF 
r {“"Station elevation n a t u m f m j 
r ti = total stations 

c TN=1X STRIPS /SP ALONG X-AXIS 
c f = <r:'. P’NAt TOTAL PTFLP, GAMMAS 
C F I = J NCI. I NAT I ON OF V , II LG 
C F 0 = OECL t NAT I ON OF F 
C 

C s = CVJ ANGLE FROM f •> TO SM* RADIANS 
C R * ACUU ANCLE BETWEEN SS AND X 
C I. = 1 NCR F mc jj t JNnex 

(»<i*ei**tst^*f<i***f«t*s«**«*tt**#<i <t e *e <<♦$«**** »***♦* 

C. 

c 

0 T u T!snN v( ion , moi i ,GHot i , Ai ion 

1 R p A ,1 <I* 7 »PNri = 105 ) JM,XK ,n*OD,W, SS,SM»FI ,FO,F, II ,XN,OX,C 

2 FORMAT t I 3 , 2 X,f *.C, 3 (?X*F 4 . 0 ) ,M 2 X,F 4 .T > , 2 X,P 5 . 0 , 2 X, 13 , 2 X, 
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17 X« C S.O* 2 X*F** 0 * 2 X*FS.O) 

F I :f 1 / F7. 79SP 

OXs^xeiOO. 

c^o^ico. 

Ci=o*ioo. 


c?=(o*oni*ioo. 

R=A3SI ISM-90, i-ssi 

3=3/57.2958 

S=SM-FU 

IF ( SM.LT.FOJ S* 3 60.— S 
S = S/57.2‘»5a 
p=rr>si f n *si\u s> 

0=STNin i 

SP=mx/cos(m )*xn 

M n-G. 
no 17 1=7,1! 

17 At n=M i-ii + sp 
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X= A I II/?)-IW/(Cf)SIRl*?*J »-{0X/tCOS(R 1*2-1 ) 

IJ = 1 1* I F I X( XN1 
on in 1 = 1 , tJ 
X I =FLOAT I I ) 

XP=X! *DX 
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x*xf lox/irrisis 1*2. > ) 

CO TO 20 

19 r.ONTimip 

20 XI. ■ 2 .*XK*0X*F/( I0.**l l) 

j»j = if i xivn Fin xioxj 

no A? I = 1 , 1 T 

vm=o. 

Tt r i®o. 

32 om=o. 
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XX=f ( X- A ( 1 1 1 *CtlS I R ( 1MI_*DX)-IDX/2»I 

n=ci-c 

Z2=r.7-c 
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37 XR 2 = ( XX*P2 I *( Z2**2> . 

33 XM 1 * x X / X R 1 
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1 XN1 = / 1 /Xltl 
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47 CONTINUE 
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?1 r-DR'-’AT <’ ',FS.y f *>, METERS, OF SUSCEPTIBILITY %F6.0,‘ X 10 

1-S I-MU/CC, STRIKING ’.FS.U 
W-ITF I*,*?! DX 
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APPENDIX E.3 

Total Magnetic Anomalies (in Gammas) at 20 Meter Intervals over Wedges of 
1250 x 10* emu/cc Material (from Computerized Two-Dimensional Models) 
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